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POLARIZATION OF THE CONTINUOUS X-RAYS FROM 
SINGLE ELECTRON IMPACTS! 


By BaLeBaiL DASANNACHARYA 
RYERSON PHYSICAL LABORATORY 
UNIVERSITY OF CHICAGO 


(Received November 8, 1930) 


ABSTRACT 


Experiments are described on the polarization of continuous x-rays from targets 
of aluminum sheets varying in the thickness from 6 to 250X10-§cm. The velocities 
of the exciting electrons corresponded to voltages from 27 to 57 kilovolts. It was 
found that the polarization increased exponentially with decreasing thicknesses, sug- 
gesting perhaps a nearly complete polarization for a thickness of about 6X10-*cm. 
The polarization diminished with an increase in the velocity of the exciting electrons. 

The theory of the polarization of continuous x-rays developed by Sugiura seems 
to be in conformity with the experimental results obtained in the experiments de- 
scribed below. 


HE old problem of the radiation due to the stoppage of electrons, 

attacked by Stokes and J. J. Thomson immediately after Réntgen’s 
discovery of the x-rays, was extended by Sommerfeld? by the methods of the 
classical theory. Recently Sommerfeld’ extended his previous work by apply- 
ing the mathematical methods of wave-mechanics to cases where the elec- 
trons were stopped in a single process (very thin anticathodes). 

This new theory was well supported by the measurements of Kulen- 
kampff‘ on the angular distribution of the radiation produced in an aluminum 
target of thickness 6.5X10-' cm. The intensity in the forward direction of 
motion of electrons was from these experiments extrapolated to be zero. 
Duane® working with mercury vapour found that the intensities J; in the 
forward direction and J, at right angles were in the ratio of 1 to about 5.5. 
Defining the percent polarization as (J;,—J;/I,;+J:)X100, this might be 
interpreted as a polarization of 70 percent. Sommerfeld consequently made a 
generalized assumption according to which only those electron beams emerg- 
ing from the target with zero velocity should give rise to complete polar- 
ization, i.e., the radiation corresponding to the short wave-length limit of the 
continuous spectrum should be completely polarized. 

1 A short account was given at the Chicago meeting of the American Physical Society, 
Nov. 27, 1929. Phys. Rev. 35, 129 (1930). 

2 A. Sommerfeld, Phys. Zeits. 10, 969 (1909). 

* A. Sommerfeld, Proc. Nat. Acad. Sci. 15, 393 (1929). 


‘HH. Kulenkampff Ann. d. Physik 87, 597 (1928). 
5 W. Duane Proc. Nat. Acad. Sci. 14, 450 (1928). 
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1676 BALEBAIL DASANNACHARYA 


P. Kirkpatrick® calculated from his experiments with a thick tungsten 
anticathode a polarization of 11 percent for the short wave-length limit with 
a thick tungsten target. P. A. Ross’ using a method of balanced double filters 
obtained nearly complete polarization at the short wave-length limit with a 
thick tungsten target. Wagner and Ott’ reflected x-rays coming normal to 
the direction of the electron stream at the surface of a crystal and observed 
the intensity of the radiation thus reflected in directions along and at right 
angles to that of the electron beam. They found a maximum polarization 
of 47.5 for the short wave-length limit. Thus results so far recorded fail to 
show any general agreement. 

It appeared therefore desirable to make some direct measurements on 
the polarization of continuous x-rays from thin targets, like aluminum foils of 
various thicknesses, with electrons accelerated by different potential drops. 


APPARATUS 


An electron stream (see Fig. 1) from a hot spiral cathode fitted with a 
focussing cup passes through two openings A and B, both of which form a 
single element and are earthed. The distance between A and B is 5 cm. The 





alls 














Fig. 1. Diagram of apparatus. 


opening A is in the molybdenum disk facing the electron stream. B is in an 
iron sheet welded into the cylinder which carries the molybdenum. A centi- 
meter behind B the electron beam meets the aluminum foils stretched across 
two rings G and H. The beam after passing through the aluminum passes 
without striking the sides of the lower ring ZZ into a hollow cylinder Z and 
through a miliammeter to the earth. Circular openings D and E in lead sheets 
are so placed and are of such a size as to prevent x-rays formed at B from 
falling on the graphite block F. The effectiveness of the shielding was tested 
by observing the effect on the readings of the electrometer when the alu- 
minum foils between G and JJ were rotated away from the path of the elec- 
tron beam. Further, when an aluminum foil is in the path of the electron 
beam, it easily gets punctured at the place where the beam meets the foil. 
The size of this puncture, about 4 mm in diameter, and the place where it 
occurs, gives good indication of the path of the beam in the foil system GH. 

§ P. Kirkpatrick, Phys. Rev. 22, 226 (1923). 

7 P. A. Ross, J. Opt. Soc. Am. 16, 375 and 433; Phys. Rev. 28, 425 (1926). 

* E. Wagner and P. Ott, Ann. d. Physik 85, 425 (1928). 




















POLARIZATION OF X-RAYS 1677 


Pinhole photographs were also made for further confirmation. The openings 
in the diaphragms in front of the ionization chamber on a Bragg spectrograph 
were one square centimeter. 

In order to avoid charge developing on the inside walls of the x-ray tube, 
and deflecting the electron beam, the inside surface from A to the end of 
the tube at L was platinized with platinum paint. Before the paint was 
heated it was nearly completely removed from that part of the wall through 
which the x-rays passed to the graphite F, in order to prevent any filtering 
action of the x-rays by the platinum layer. 

The ionization current was measured in a Compton electrometer worked at 
a deflection of the reflected spot of about 30,000 mm per volt, with the scale 
at a distance of 2 meters from the mirror. The current from a transformer 
was rectified through a kenetron before being admitted into the x-ray tube. 


EXPERIMENTAL RESULTS 


The intensities 7; and J; of the x-ray beam emerging from the graphite 
along the direction of the electron beam and at right angles to the same were 
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measured. Correction due to stray radiation when there was no aluminum 
foil in the path of the electron beam at C was found negligible. The thickness 
of the aluminum foils or plates varied from 6.4 to 250X10-' cm. For the 
thinnest aluminum foil my obligations are due to Dr. Kulenkampff. These 
foils were stretched between the rings G and H/, fastened with a paint made 
of fine aluminum powder and a little varnish dried over hot air. 

The experiments of June 23rd (see Table I) showed that the polarization 
increased with the diminution of thickness as well as by a diminution of the 
accelerating voltage of the electrons. The value for 57 KV seems rather too 
low and that for 45 KV rather too high. The effect of voltage and thickness 
seems however to be unmistakable. The experiments of October 18 were 
carried out with one of the thinnest foils used previously, nearly 48.0 10 
cm, and three others. The effects of voltage and thickness (see Fig. 2) are 
here more clearly marked. 

A rough extrapolation for smaller thicknesses indicates that complete 
polarization may result for a thickness of about 6X10-*cm. This thickness 
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TABLE I. 
Date Accelerating poten- Thickness of the Percentage polarization 
tial in kilovolts. aluminum plates, in (2e—I1) / (Ut +dh) X 100 
10-5cm 
a — Fe Nancae : 250 _ -_ | —— * _ 
ae ) 144 pad ao 
96 23.7 
June 23,1929 | 48 | 28.5 
45... 48 | 26.3 
57 48 | 16.3 (?) 
48 . 3 
( 27 19.2 37.1 
12.8 43.5 
6.4 53.0 
ey cass f 12.8 . 39.0 
6.4 . 46.2 
Oct. 18,1929 | 
48.0 > ae 
| 42 19.2 29.8 
12.8 33.8 
6.4 46.7 
Pe sea 12.8 31.5 
, ae | _* oo _ 33.0 





happens to be about the same as that (4X10~* cm) for which the electron 
keeps its initial direction of motion, as calculated from theories of the diffu- 
sion of electrons through matter® for an average accelerating potential of 
34 KV. According to the diffusion theory, however, the length of path along 
which the electron keeps its initial direction should increase with increase in 
speed, which would mean greater polarization for greater electron speeds. 
Since this is contrary to the results of the present experiments, it would ap- 
pear that the approximate agreement between the length of the undeviated 
path and the thickness for complete polarization is probably fortuitous. 

According to Sugiura,!® who has recently developed a theory of the pro- 
duction of continuous x-rays somewhat on the lines followed by Sommerfeld, 
one should expect that the polarization at right angles to the direction of 
motion of the electron beam even for vanishing thickness of the target to be 
incomplete. The polarization should be greater for rays corresponding to 
electrons of smaller speed, a result which is in keeping with the findings of 
our present experiment. 

It isinteresting to note that Kulenkampff" has recorded a value of 45 per- 
cent for the polarization with electron beams of 36 KV presumably for the 
thickness of aluminum the same as the least one in our present experiment. 


® Handb. d. Exptl. Physik, Kathodenstrahlen by Lenard and Becker pages 333, 379. 
Akad. Verl. Buchh., Leipzig, 1927. 

10 Sugiura Sci. Papers I. P. C. R. 11, 251-290 (1929). 

1H. Kulenkampff, Phys. Zeits. 30, 513 (1929). 
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POLARIZATION OF X-RAYS 1679 


Our value for the same is 47.5, within errors of experiment the same as the 
value obtained by Kulenkampff. 

Duane” did not find any great difference between the polarizations with 
jets of mercury vapour as anticathode, the density of the vapour being 
varied from one to four. On an average the polarization decreased from 
51.4 to 48.6, a result not contradicting our own findings, though Duane is 
inclined to believe that the difference perhaps is not very much greater than 
the experimental errors. Secondly, it may be that in Duane’s experiment we 
are nearly at the maximum of polarization; or it may also be that in the case 
of mercury the change of polarization with thickness is small. Thirdly, a 
comparison of Duane’s value (taking it to be nearly the maximum for mer- 
cury) with our own for Al extrapolated to smaller thicknesses than actually 
used, would lead us to suppose that an increase in the atomic number of the 
anticathode would result in a decrease in the polarization. 

This last result as well as results of Duane just referred to would be hard 
to reconcile with the almost perfect polarization obtained by Ross with his 


‘double filter method, and using thick tungsten as anticathode. Perhaps the 


balancing of Ross’s filters, which may have been perfect, or the error in the 
same inappreciable, when the differences are large, would cease to be so when 
the difference to be obtained is about of the order of the inaccuracy. As 
against Ross,’ Wagner and Ott® using platinum as anticathode could get 
only 47.5 as the polarization for the short wave-length limit. 


A RouGuH PICTURE OF THE PRODUCTION OF POLARIZATION IN CONTINUOUS 
X-Rays 

The polarization may depend on the following: 

(1). Initial direction of the impinging electron. 

(2). Its final direction, and probably also on, 

(3). The atomic number of the material of the anticathode, and 

(4). Electron configuration round the nucleus. 

If the energy of the impinging electron is completely utilized in the pro- 
duction of the continuous x-ray, as would be the case for radiation at the 
short wave-length limit of the spectrum, and consequently (2) is negligible 
and if (3) and (4) should happen to be negligible too, the polarization may 
be complete. On the other hand, as we come to frequencies smaller than the 
limiting one, we may think of them as being produced by electrons leaving 
the atom with appreciable speeds in directions which however are arbitrary. 
Consequently (2) becomes increasingly appreciable so that the electric vector 
would deviate progressively from the initial direction. If we are to account 
for more complete polarization for smaller velocities, we have to assume that 
the probability is greater that an electron of a smaller speed would emerge 
from the atmosphere of an atom with comparatively greater loss of energy. 

I wish to thank Professor A. H. Compton for suggesting the problem, 
for his continued encouragement and for his kindness in permitting me to 
work in his laboratory as a visiting foreign guest of the University. 

122, W. Duane, Proc. Nat. Acad. Sci. 15, 805 (1929). 

13 P, A. Ross, J. Opt. Soc. Am. 16, 375 and 433. Phys. Rev. 28, 425 (1926) 
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SPECTRUM OF THE RADIATION FROM A HIGH POTENTIAL 
X-RAY TUBE 
By C. C, LAURITSEN 
CALIFORNIA INSTITUTE OF TECHNOLOGY 


PASADENA 


(Received November 7, 1930) 


ABSTRACT 


A spectrograph of the Seemann type has been constructed for the purpose of 
investigating the radiation from the high potential x-ray tube at the California 
Institute. A typical spectrogram obtained with 600 kilovolts on the tube is presented. 
The photometer record shows a continuous spectrum with its maximum intensity at 
about 200 kilovolts and a short wave-length limit in the neighborhood of 600 kilo- 
volts. The range covered is roughly from 100 to 20 x-units. It is proposed to use the 
apparatus for determining absorption coefficients by photographing the spectrum of 
radiation which has passed through an absorbing screen. No anomalies of any kind 
have been observed so far. 


APPARATUS 


HE high potential x-ray tube at the California Institute of Technology 

has recently been rebuilt and equipped with a hot cathode and a tungsten 
target. A description of the tube and its housing has been presented by 
Lauritsen and Cassen.! The tube in its present form operates satisfactorily 
at 600 kilovolts and it has been deemed advisable to investigate the available 
radiation before attempting to go to higher potentials. 

In order to investigate the region of the x-ray spectrum from approxi- 
mately 150 kilovolts and up, a crystal spectrograph was constructed following 
the principle described by Seemann,? Siegbahn’ and others. High precision 
cannot be expected with any reasonable time of exposure by this method 
since the whole of this region of the spectrum lies within an angle of less 
than one degree and thick slits are required because of the great penetrating 
power of the hard radiation. On the other hand, the method is convenient 
and sufficiently precise for the approximate determination of the short 
wave-length limit as well as of the general distribution of intensity in the 
spectrum. Also, if any prominent lines, bands or absorption edges or other 
unexpected irregularities exist in this region, they should be found most 
conveniently by this method. 

The spectrograph consists essentially of a vertical slit 0.9 mm wide in a 
lead block 4 cm thick, in front of which a rock-salt crystal is placed. The 
crystal, slit and photographic plate are all rigidly mounted on a long arm 
which may be rotated through an angle of one degree on each side of the 
center, the rotation taking place about an axis through the vertical center 
line of the slit. Since the focal spot has the appearance of a thin horizontal 


1 C. C. Lauritsen and B. Cassen, Phys. Rev. 36, 988 (1930). 
2 Seemann, Phys. Zeits. 18, 242 (1917). 
* M. Siegbahn, Phil. Mag. 2, 639 (1919). 
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disk approximately 5 mm in diameter, it is clear that radiation will reach the 





central planes of the crystal under all angles from zero to 50 minutes, pro- 
vided that the distance from the focal spot to the crystal does not exceed 
35 cm and that the spectrograph arm is set at an angle of 25 minutes on either 
j side of the center line through the focal spot and the slit. Under these con- 
ditions the undeviated part of the light will produce an image of the focal 
spot as seen through the slit, the whole of the image being located on one side 
‘ of and adjacent to the center line of the photographic plate. The reflected 
portion of the light will appear as a spectrum on the other side of the center 
line. The short wave-length limit is thus given by the distance from the cen- 
ter line to the edge of the continuous spectrum, and lines of equal wave- 
length are parallel to the center line. Since different wave-lengths originate 
from different lateral regions of the focal spot and since the intensity is not 
absolutely uniform throughout the whole of the area, it is not possible to 
obtain the intensity distribution in the spectrum directly with any degree 
of precision, Correction can, however, be made for this lack of uniformity 





Fig. 1. Double spectrum 600 kv peak. 


if the direct beam is reduced by means of a filter to an intensity comparable 
to the intensity of the spectrum. This is so because the direct beam is a true 
image of the focal spot and the spectrum is a specular image of the direct 
beam, except for the wave-length selectivity, and there appears to be no good 
reason why the radiation from different paits of the target should differ 
appreciably in hardness if the intensity is the same. 


The photographic plate is placed 107.5 cm from the center of the slit. 
In spite of this comparatively great length, the resolution is not high on 
account of the wide slit, but it is sufficient for the present purpose. 

In order to determine the exact zero on the photographic plate it is most 
convenient to photograph both right and left hand spectrum on the same 
plate. This is done by adjusting the spectrograph as described except that 
the direct image is blocked out completely by a heavy lead screen while one 
spectrum is being photographed. The spectrograph arm is then rotated into 
position on the opposite side of the center line and the direct image which 
now appears on the part of the plate which has been exposed is blocked off. 
The second exposure thus gives a spectrum on the part of the plate which 
was covered up during the first exposure. The result is two spectra which are 

symmetrical with respect to the center line or true zero, Fig. 1 is a positive 





reproduction of a double spectrum taken in this manner. 
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A dark line may be seen running diagonally through the spectrum in 
each of the spectra of Fig. 1. The origin of this line becomes apparent if we 
consider the geometry of the arrangement. The spectrum shown is due to 
reflection from vertical (100) planes of the rock-salt crystal, but if the crystal 
is adjusted so that the horizontal (100) planes are parallel to a horizontal 
plane through the focal spot and some part of the photographic plate, then 
it is clear that from each point of the target there will be a ray of a given 
wave-length which will make the same angle with the horizontal as with the 
vertical planes. The intensity of this ray will therefore be divided between 
the spectrum shown and a spectrum which falls above, below, or within the 


direct image. 


~ 
on view 


Fig. 2. 


This second spectrum is of the type first obtained by Rutherford and 
Andrade and differs from the first in that the resolution is slightly less. The 
paths of a ray are shown by the lines a, ) and c in Fig. 2, and the heavy lines d 
in the end view are the intersections with the photographic plate of all the 
ravs which fulfill the foregoing condition. It is readily seen from the figure 
that the angle between the lines d and the vertical center line is given by 


b 
tané = . 
ad + b 
In the present case, this gives 
107.5 — 
tan@é = = ().755: 6 = 37.05 
142.5 


which agrees with the angles made by the dark lines in Fig. 1 with the center 
line. 
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The intensity distribution was obtained by means of a recording micro- 
photometer. In Fig. 3 the galvanometer deflections obtained from the 
photometer records are plotted as ordinates against distance as abscissas. 
The tube was operated at 600 kilovolts peak and, as indicated in the graph, 
the short wave-length limit corresponds closely to this value. The maximum 
intensity occurs somewhat below one half of this potential, as might be 
expected since the tube is operated with alternating current and a “thick” 


Defilec! 


target is used. It should be noted, however, that the softer radiation is 
decreased somewhat in passing through both the 6 mm steel wall of the tube 
and the crystal. The lack of symmetry which is apparent in the graph is 
due partly to the aforementioned non-uniformity of the focal spot and 
partly to a slight difference in exposure. 


RESULTS 


From the photograph as well as from the photometer record we may 
conclude that, within the limits of the resolution used, there are no unex- 
pected irregularities in the spectrum from tungsten in the region covered. 
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At the present time work is in progress to determine as accurately as possible 
how the spectrum is modified by different absorbers. The right half of Fig. 4 
shows the spectrum after passing through lead screens. One screen, 0.28 in. 
thick, covers the whole image, and in addition the central portion is covered 
by a second screen 0.14 in. thick. The direct beam was decreased toa suitable 
intensity by a steel block 2 in. thick. 





Fig. +. Left, direct image. Right, spectrum. 


The photometer records of spectra taken in this way show accurately how 
the spectrum is modified by any given absorber and it should therefore be 
possible to obtain from a single pair of records the absorption coefficient as 
a function of the wave-length. 

If nuclear absorption levels or other unexpected irregularities exist in this 
region of the spectrum, they should be readily detected by this method. 
The spectra obtained so far and also the modifications due to absorbing 
screens of aluminum, iron, copper and lead are roughly what would be 
expected. There is no indication of any sudden changes in intensity with 


wave-length. 
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MULTIPLE SCATTERING IN THE COMPTON EFFECT 


By Jesse W. M. Dumonp 
CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA 


(Received November 5, 1930) 


ABSTRACT 


In the experimental study of the spectral distribution of x-radiation scattered 
by light elements it has always, up to the present, been assumed that multiple scat- 
tering could be neglected. Recent improvements in experimental technique however 
make it possible to suppose that multiple (especially double) scattering may now be 
detectable from large scattering bodies. Multiple scattering may (1) affect the 
breadth of the modified line, (2) change the structure of the modified line, (3) distort 
the background in such a way as to render measurements of shift unreliable. It is 
therefore valuable to analyze the effect of multiple scattering in case some of the recent 
mutually discordant experimental results may be explained and harmonized in this 
way. 

Assuming initially monochromatic radiation, scattering of any degree of multi- 
plicity contributes a spectral band whose wave-length limits are here determined and 
discussed. The Breit, Compton, Jauncey formula taking polarization into account is 
adopted for the purposes of calculation of modified intensity. The results are thus fairly 
accurate for hard radiation scattered from very light elements. For softer radiation 

modified scattering at small scattering angles is greatly reduced below the values 
given by the Breit, Compton, Jauncey formula and unmodified scattering appears. 
The effect of this on the results of this paper is discussed in a qualitative way. 

The case of double scattering from a spherical scatterer is computed in complete 
detail and a formula for the ratio of double to single scattering is derived. Curves of 
the spectral distributions due to double scattering are shown. The dependence of the 

total doubly scattered intensity on the primary scattering angle is plotted. 

The natural width of the modified line is neglected throughout these calculations. 
Absorption in the scatterer is also neglected. The ratio of double to single scattering 
for a spherical scatterer observed under any given angle is proportional to the radius 
of the scatterer (neglecting absorption) and is given by 

Doubly scattered energy 9 


- = —orR(@) where 
Singly scattered energy 32 @) 





o is the linear scattering coefficient for the material of the scattering sphere, r the ra- 
dius of the sphere, @ the angle under which single scattering occurs. R(@) never differs 
greatly from 2.5. 

Triple scattering is negligible in comparison to single scattering. 

Twice modified doubly scattered radiation may contribute a faint asymmetric 
line or edge at the shifted position 


Ad =2(h/mc)(1+ cos 36) 


For hard radiation when the Thompson formula (with Breit correction) applies 
with fair accuracy to modified scattering alone twice modified doubly scattered radi- 
ation contributes a spectral band of breadth 

4(h/mc) cos 30 


Once modified doubly scattered radiation may cause a slight broadening of the 
Compton line except in regions near @=0 and @=180°, 
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INTRODUCTION 


N THE experimental study of scattered x-radiation it is impossible 

completely to eliminate multiple scattering. X-radiation whose direction 
is defined more or less closely within some solid angle 6Q; is incident upon 
a scattering body. The apparatus (a spectrograph, ionization spectrometer 
or what not) is arranged so as to receive scattered radiation whose direction 
again is more or less closely defined within some solid angle 6Q.. The angle 
6 between the initial and final directions has thus assignable limits of in- 
homogeneity and is called the angle of primary scattering. A ray cannot be 
prevented, however, from suffering any number of scattering processes be- 
tween its entry along the first direction and its exit along the last. It is the 
purpose of this paper to discuss theoretically the effect of such processes on 
the spectral distribution of the scattered radiation. The path of a multiply 
scattered ray is a broken line (in space) of (1/+1) straight segments, where 
M is the multiplicity of the scattering. The first and last segments have the 
fixed directions of the incident and scattered beams just mentioned. The 
other segments will be called intermediate rays. Their directions may be any 
whatever. In double scattering there is but one intermediate ray. The 
directions of the entry, exit and intermediate rays are to be plotted upon a 
sphere of unit radius by means of unit vectors originating at the center of 
the sphere and terminating on its surface. These vectors have the respective 
directions of propagation of the radiation along entry, exit and intermediate 
rays for their directions. The terminii of these vectors thus locate 1J+1 
points on the surface of the sphere and represent the (./+1) directions 
taken by the M-tuply scattered ray. The first and last of these points are 
fixed upon the sphere and separated by the angle of primary scattering @. 
The remaining (JJ—1) points may have any position on the entire surface 
of the sphere. The problem of J/-tuple scattering thus involves the (2.1/ —2) 
degrees of freedom of these points. 

Each process of scattering splits an initially monochromatic ray into 
modified and unmodified components as regards wave-length. The modified 
ray suffers an increase of wave-length given by AN=h/mc(1—cos 6;) where 
6; is the angle of scattering. There exist therefore (1/+1) cases in m-tuple 
scattering for the ray may be scattered M times but only modified 0, 1, 2, 3, 
etc. up to M times. Thus in double scattering the ray may be scattered twice 
without modifications, or modified at the first scattering but not at the 
second, or modified at the second scattering but not at the first, or modified 
twice successively. The two cases of singly modified doubly scattered radia- 
tion will be identical. 

Unfortunately the law governing the partition of the scattered radiation 
between modified and unmodified rays is complicated and not very well 
known. For sufficiently hard radiation scattered by sufficiently light atoms 
however the unmodified radiation is negligible at nearly all scattering angles. 
In this paper only this limiting case will be discussed analytically in detail 
and therefore the total shift of the original wave-length will be the sum of 
all the separate shifts occurring at the successive scattering processes. The 
breadth of the shifted line will be neglected. 





























MULTIPLE SCATTERING IN COMPTON EFFECT 


WAVE-LENGTH LIMITs OF MULTIPLY SCATTERED RADIATION 


It is immediately possible to make the following obvious generalizations 
if M is both the number of scatterings and the number of shifts of a given 
ray: 

1. When J is odd the maximum shift occurs for a primary scattering 
angle of 6=180°. 

2. When M is even the maximum shift occurs for @=0°. 

3. Minimum shift for any multiplicity is always zero and always occurs 
at a primary scattering angle # =0. 

4. Maximum shift is always 2M(h/mc) where M is the multiplicity of 
scattering. 

DOUBLE SCATTERING 
Referring to Fig. 1, on the sphere, A is the direction of the entry ray, 


B the direction of the intermediate ray and C the direction of the scattered 
ray. The arc AC measures the fixed primary scattering angle @ and is bi- 





Fig. 1. Illustrating spherical coordinate system (a, @) appropriate to 
the problem of double scattering. 


sected at P. P is taken as the pole of a system of spherical polar coordinates 
(a, @) which describe the location of the intermediate ray B whose two de- 
grees of freedom are the independant variables of the problem. 

That such a system of coordinates (a,) is appropriate to the problem 
will be seen from the following considerations: 


The total shift 
h 
Ad = —{2 — (cos 6: + cos 42) } (1) 
mc 


If the letters A, B and C represent unit vectors, then 
Ad = h/mc{2 — (A-B + B-C)} 
h/mc{2 — B-(A+C)} 
2h/mc(1 — cos 36 cos a) (2) 
(since| A + C| = 2 cos 36) 


Hence regions of constant total shift AX correspond to directions of B lying 
in zones on the sphere described about P as a pole. 
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WavE-LENGTH LIMITs OF DouBLY SCATTERED 
DOUBLY SHIFTED RADIATION 


Since cos a may range from 1 to —1, doubly scattered radiation may 
have a range of shifts defined by the inequalities 
2h/mc(1 — cos 30) < AX < 2h/me(1 + cos 364) (3) 
The limits of this range are evidently equidistant from the point A\ = 2h/mec 
and the total spectral width of the range is 4h/mc cos 30. Thus for a primary 
scattering angle of zero degrees double scattering contributes shifted radia- 
tion over a range of shifts from zero to 4h/mc while for a primary scattering 
angle of 180° double scattering contributes shifted radiation falling exactly 
at 2h/mc and extending over no range whatever. It is convenient to de- 
scribe the spectral distributions caused by double scattering in terms of a 
variable x, which is defined as x = (AX)/(2h/mc). This variable which is pro- 
portional to the shift in wave-lengths is a pure number. 


THe LAw GOVERNING INTENSITY OF SCATTERING 


In order to compute the distribution of intensity due to double scattering 
over the range determined for any particular primary scattering angle, 
some law must be adopted to express the scattered intensity as a function of 
the angle of scattering and of the polarization of the radiation. Kallman and 
Mark! have shown in a beautiful experiment that both shifted and unshifted 
radiations are completely polarized by scattering at 90°. In double scattering 
therefore the beam becomes partially polarized in the first scattering process 
and this fact is of importance in determining the intensity of scattering in 
the second scattering process. 

We therefore assume in this computation that each scattering process 
gives modified radiation having the intensity and polarization dictated by the 
classical Thompson theory of scattering and then apply the Breit correction 
to this by multiplying by a factor (1+ a vers 6;)~* (1+a vers 62)~* (c.f. Comp- 
ton “X-Rays and Electrons” pp. 304-305). Fortunately (1+a vers 6,)-% 
(1+a@ vers 6.)-* can be replaced by (Ao/A:)® (A:/A2)® and hence by (Ao/Az)3 
where Xo is the primary wave-length, A, the wave-length after the first scat- 
tering, and Az the wave-length after the second scattering. It is thus a very 
fortunate fact that the Breit correction factor depends only on the ratio of 
the initial and final wave-lengths in double scattering and not at all on the 
intermediate wave-length. Thus after the doubly scattered spectral dis- 
tribution has been computed by the classical Thompson formula and the 
distribution is obtained as a function of AX or x the Breit correction can be 
applied to the curves by simply multiplying each ordinate by the ratio 
(Ao/A2)* applicable to its abscissa. 

The above outlined method should give fairly accurate results for hard 
radiation and light scattering atoms where the presence of unmodified radia- 
tion and diminished modified radiation for small scattering angles can be 
neglected. 


1 Kallman and Mark, “Uber einige Eigenschaften der Comptonstrahlung Zeits. f. Physik 
36, 120-142 (1926). 
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Unfortunately the case of say MoKa radiation scattered by graphite 
is not such a one. The author has not succeeded in obtaining an analytic 
solution for double scattering rigorously applicable to this case taking ac- 
count of the unmodified scattering and diminished modified scattering. 
Moreover the really great analytical difficulties incident to such a solution 
seem scarcely worth overcoming as a good qualitative idea of the spectral 
intensity distribution can be obtained by a rough examination of the effect 
of reduced modified scattering on the case which is here solved. 

The shape of the scattering body also plays a part in determining the 
spectral distribution of doubly scattered radiation. It is of course practically 
impossible to take account of this factor and accordingly the solution here 
given can again only be applied to real cases qualitatively. In order to con- 
struct an ideal case in which the shape of the scatterer does not complicate 
the analysis we consider first only one electron at the center of a spherical 
scattering body. This electron receives primary radiation along-the direction 
A and scatters in various directions B to all the other electrons in the scat- 
terer which in turn scatter the radiation along the exit direction C. It is the 
spectral distribution of this radiation as affected by two wave-length modifi- 
cations which we propose to investigate. Subsequently the reasoning is 
readily extended to initial scattering by all the electrons in the spherical 
scattering body. 


COMPUTATION OF SPECTRAL DISTRIBUTIONS DUE TO DOUBLE SCATTERING 


Referring to Fig. 2, let E,, and E,,; be the primary component electric 
intensities E,; being normal to the plane of the first angle of scattering 6, 





Fig. 2. Illustrating the components of the electric vector of the incident, 
intermediate and exit rays. 
and £,; parallel to that plane. At a distance 7, from the first scattering elec- 
tron along the intermediate direction B the scattered radiation will have 
electric components E,2 and E,. again normal and parallel respectively to 
the plane of @; given by the Thompson classical theory as: 


E,.1e? tala 
E,: = — —— (4) 
rymc* 
: E,,e? 
Ewe = — —— cos 6. (S) 


rymc? 


A 
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Let B represent the supplement of the angle B of the spherical triangle 

A BC. We now resolve the electric vector whose components are E,. and 

ka. along two new directions normal and parallel respectively to the plane 
of the second angle of scattering #2 and obtain F,,3 and Eu; 

ku3 = Exe cos B + Eye sin B (6) 

E43 = Ea.cos B — E,esin B, (7) 


II 


The scattered radiation having these components E,;3; and F,3 is now scat- 
tered again by an electron at distance r, from the first, the new scattering 
angle being #. At distance 7, from the second scattering electron along the 
direction C the doubly scattered radiation will have electric components 
E,4 and Fas 


. E.,,3e? ; 

I _— : 2 (8) 
romc~ 
Ea3e" 

E., = - cos 6e. (9) 
rome 


We now express /,,,; and £,, in terms of FE, and E,; and evaluate the square 
of the final electric intensity - 


e 


E; = E..4* ae E.¢ = E, ” cos" B Ea\" cos” 6; sin? B 





r*ro*mic® 
-L 17 OY Dee cos B sin B COs 6; + cos” A4( E.\" cos” A; cos? B (10) 
+ E,,? sin? B — 2E,,E,, cos 0; cos B sin B)). 


For unpolarized primary radiation we average over all possible orienta- 
tions of £, about the primary beam giving all orientations equal weights. 
This means that 1En?} = | Ea?} =1Ef and jE x =(), the curly bracket 
here being used to indicate the average. The average doubly scattered 
intensity is then given by 
j » ol ] + 9 e8 9 
,£7; = —E;* ——— [cos? B 

2 ryre*m'c§ 


+ (cos? 6; + cos? @s) sin? B + cos? 6; cos? 62 cos? B] (11) 
Now from a well-known formula applying to spherical triangles we can 
express B in terms of 6;, 62, and @ 
cos 86 — cos 8; cos A. 


cos B = —— —- . (12) 
sin 6; sin @» 


By appropriate substitutions employing this formula we obtain the doubly 
scattered intensity in terms of 4), 02 and @ as 
{po 1 +9 e* 9 9 
E,?} = —E,? ————-|cos? A; aa CcOs* A» 


2 r\*ro*m'ic® 


+ cos? @ — 2 cos 6 cos 6; cos A. + cos? 6; cos? 62). (13) 
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We now transform this to an expression in a@ and @ the spherical coordi- 
nates of our problem by the use of the relations from spherical trigonometry 





cos 6; = cosacos}@+sinasin}6@cos¢ (14) 
cos 02 = cos a cos } @ — sina sin 36 cos ¢. (15) 
This gives 
{rol 1 7 4 é 9 * 1 ° 9 * 9 1 
ae, @ —o—— [2 cos? a cos? } 6 + 2 sin? a sin? } 6 cos? + cos? 6 
2 r\7re"mic 


— 2 cos 6 cos* a cos? } 6 + 2 cos @ sin? a sin? } 6 cos? 





16 
+ (cos? a cos? }@ — sin? a sin? 3 @ cos? ¢)?] (16) 
or: 
1 e ; " 
I; = — I; oe —|v(a, 6, ¢)]. (17) 
2 rire*mics 


I, is the doubly scattered intensity. It is the energy per cm.? per second 
at a distance 72 from the second scattering electron which in turn is at a 
distance r, from the first scattering electron. 

The exit ray is confined to some solid angle AQ; by the construction of 
the apparatus for studying the scattered radiation. If we multiply the last 
expression by the area r.°,AQ; we obtain the energy scattered per second into 
the solid angle AQ; 

1 eAQ 
I r2AQ; = —I ,——-|y(a, 8, ¢) J. (18) 
2 = r,?m4c8 

Consider now a spherical shell of radius r; and thickness dr, described 
about the first scattering electron as a center. Let the volume density of 
electrons in the scatterer be p. In the solid angle between a and a+da, 
¢@ and ¢+d¢ there will be pr,’ sin a dr;dad@ electrons. The energy per second 
entering the solid angle AQ; due to all the electrons just mentioned will be 

1 e®8AQspdr; 


T4r2?AQspr,? sin adrydadé = —I ————[y(a, ¢, 0) sin adad@]. (19) 
2 m'¢§ 





Integrating this from r=0 to r=r, the distance from the first scattering 
electron to the furthermost electron in the body in the direction (a, 9), 
we have 

1 CAQ3Pr ag 


—I,———[ y(a, ¢, 8) sin adado}. (20) 
2 m'¢8 


The total shift after two scatterings is independant of ¢ as has already 
been pointed out. We therefore obtain the total energy going into a zone 
between a and a+da, all of which contributes to one and the same shifted 
position in the doubly scattered spectrum by integrating around the zone 
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from ¢=0 to ¢=27. We thus obtain the energy per second passing out 
through the solid angle AQ; after first scattering by one electron at the 
center of the spherical scattering body of radius, 7, and then scattering by all 
other electrons in the body between the cones of half angle a and a+da. 
It is 

S8AQ, AQ. 
LPs da, 0)| = fy - oe | 2202 — 2 cos 0) cos? a cos? 3 8 


2 mic 2 mic¢8 


+ 2(2 + 2 cos 6) sin? a sin? $6 + 22 cos? 6 + 2x (cos? @ cos? } 6)? (21) 


> 
3 7 
— 2m cos? a cos? 3 @sin? asin? 36+ ri (sin? a sin? 3 0)? | sin ada. 


We must now transform this expression to one in terms of x and dx the 
independant variable used to express the shift in the description of spectral 
distribution. x=(AXA)/(2h/mc). The shift is given by 


AX = 2h/mc(1 — cos 5 6 cos a) 


hence 
x = (1 — cos } 6 cos a) (22) 
dx = cos} @sin ada (23) 
sin ada = dx/cos } 0 (24) 
. (1—.x)? _ Se ae i 
cos*a@ = ———-; sin* a sin* 3 @ = sin* >0— (1 — x)? tan? 3¢@. 
CcOs* ; ] 


Making these substitutions we obtain the total energy per second in the 
spectral region between x and x+dx; it is 


1 = we®AQspr. ’ 1 ec®&AQsprx 
Idx = BF rears [ yv(x, é)ldx = —I,— ma (2 cos? 6 





2 micS 2 mic8 
9 , 
3 3 
+ sin?@ + ry siné} 6) + (2sin? 36 —- > sin? +6 tan? 3 6)(1 — x)? (25) 


3 dx 
+ (2 + 2tan? 360+ rs tan’ 3 s) (1— | dees 


cos3 6 


This expression refers to one electron only at the center of the spherical 
scattering body as the first scattering agent. This same electron will by 
single scattering into the solid exit angle AQ; scatter energy 


ef*AQ 





~(1 + cos? 6). (26) 


2 m¢4 


1 
I. = —IJ, 


The ratio of double to single scattering is then 


u efor [y(x, 6] 





om aie = T _— neces dx (27) 
I. m*c* 1 + cos* 6 
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but the “scattering coefficient” is; 





8 e'p 
g¢nas- siren 
3 m*c4 
I] 3 (x, 6) 
ON (28) 
ip 8 1+ cos? a 
The ratio of total doubly scattered energy to singly scattered energy is 
1+cos 1/20 
| Idx . J [y(x, 6) |dx , 
—COSs 1/ 26 x “ 
— = —orR. (29) 














= or 
Ts 8 1 + cos? @ 8 


Evaluating the integral in this expression we obtain the ratio of double 
to single scattering, it is 





3 1 2 2 
or———_| 24 + —B cos? 36 + —C cost 3 | (30) 
8 1+ cos?é 3 5 


x 


2 cos?6+sin* 6+? sin* 36) 

2 sin? }@— 3 sin? 1 tan? 36) 
2+2 tan? }0+2 tan‘ 3 @) 

Both values of the bracket in the last expression and values of R are plotted 

in Fig. 3. 
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Fig. 3. Curves of the integral in Eq. (29) and also of the ratio R in that equation. 


For scattering from the central electron in the spherical scatterer to all 
other electrons, the value of rg is the radius r of the spherical scatterer. For 
initial scattering from all electrons in the scatterer however an average? 
value of {r.,} =%r must be taken. To obtain this we average the distance 
along any fixed arbitrary direction from each elementary volume in the 
spherical scatterer to the boundary of the sphere throughout the entire 
volume. Taking as volume elements of the sphere the space between two 
coaxial cylinders of radii r’ and r’+dr’ with axis through the center of the 


2 The bracket { } is used to denote “average value.” 
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sphere, the average distance in question is given by 1/(4/3)mr°f{(r?—r’2) 
4arr'dr’=2r. Now o=0.4 per cm for Mo radiation scattered from graphite. 

















o 
G=0 
Sak of wave -lergtls For A=708b XL 
708 732 7S 760 80+ 
, - 9 : 7 af . ? ’ to 7 
| | 
| 
oa 
3 + 
MI 
2 se | 
a RS Pr =. 
“> 7° = t~ ~4--4---F % 
, NK After apelynp Breik corres Tia? Yor 41 = 708 AY. — | _ 
NES 
>) ha) | 
x= 0 2 tT e oe 10 12 1# 16 1 Ze 


Fig. 4. Spectral distribution due to doubly modified double scattering 
for primary scattering angle of zero degrees. 
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Fig. 5. Spectral distribution due to doubly modified double scattering 
for primary scattering angle of 60°. 
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Fig. 6. Spectral distribution due to doubly modified double scattering 
for primary scattering angle of 90°, 
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It appears from Fig. 3 that R never differs very greatly from 2.5. The ratio 
then of total double scattering to single scattering for a spherical carbon 


A=20° 
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Fig. 7. Spectral distribution due to doubly modified double scattering 
for primary scattering angle of 120°. 
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Fig. 8. Spectral distribution due to doubly modified double scattering 
for primary scattering angle of 156°. 
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scatterer 0.5 cm in radius scattering Mo radiation is about 0.140, This would 
increase in direct proportion to the radius of the spherical scatterer were it 
not for absorption. Except at 6=180° however the doubly scattered energy 
is distributed over a considerable wave-length band. 

The spectral distribution of double scattering is obtained from expression 
(25) for scattering angles @ between initial and final beams of 8°, 60°, 90°, 
120°, 156° and 180°. These are exhibited in Figures 4 to 9. The curves both 
without and with the Breit correction computed for MoK radiation are 
shown. An interesting feature of these curves occurs when the angle of single 
scattering is 180°. At this angle the spectral distribution, due to doubly 
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Fig. 9. Spectral distribution (with reduced ordinate scale) of doubly modified 
double scattering for primary scattering angle of 178°. 


shifted scattering, narrows down to a sharp line (neglecting the natural 
breadth of the Compton shifted line itself). The presence of the term 
(cos 30)-! in the expression (25) makes that expression indefinitely large for 
6=180°. The energy associated with double scattering corresponding to the 
area under this infinitesimally narrow infinitely tall curve is of course finite. 


QUALITATIVE DISCUSSION OF THE EFFECT OF UNMODIFIED 
AND DECREASED MODIFIED SCATTERING 


Scattering of x-radiation of intermediate hardness (e.g. MoK radiation 
scattered by graphite) is of two well-known types; modified and unmodified. 
The total scattered energy follows the Thompson classical law with the 
Breit correction to a fair degree of approximation except at small scattering 
angles where the phenomenon of excess scattering occurs. Excess scattering 
is without doubt a phenomenon of interference of scattered x-radiation and 
applies therefore only to unmodified scattering (since modified scattering is 
incoherent). The total scattered energy is divided between the modified and 
unmodified types in proportions that depend on the scattering angle. Un- 
fortunately the exact analytic form of this dependence is unknown. It is 
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doubtless very complex since it depends in an intimate way on the mechanics 
of the scattering atom. Jauncey has treated the problem theoretically for an 
approximate Bohr atom model. In Fig. 10 the work of Ross and Woo 
relative to MoK radiation scattered from graphite is collected in one curve 
showing the ratio of modified scattered radiation to unmodified scattered 
radiation as a function of scattering angle. Fig. 11 shows Hewlett’s observa- 
tions on total scattering. Applying the ratio of the previous curve to this we 
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Fig. 10. Combined observations of P. A. Ross and Y. H. Woo of the ratio of modified to 
unmodified scattered intensity as a function of scattering angle. 
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Fig. 11. Hewlett’s observations of total scattering and the partition of this into modified 
and unmodified types in accord with Fig. 10. 


obtain the two other curves which give the modified and unmodified scat- 
tered energies as a function of scattering angle. Note that below 90° modified 
scattering diminishes rapidly and disappears completely at about 20°. In 
this same region unmodified scattering becomes very strong. 

To take account analytically of the effect of these phenomena on double 
scattering would doubtless be an extremely difficult task and one which 
would certainly not be warranted at present by the utility of its results. 
The main purpose of the present investigation is to obtain a qualitative idea 
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of the effect of double scattering regarded as an unavoidable impurity in the 
experimental study of single scattering. In this paper therefore we discuss 
only in a brief qualitative way the effect of excess unmodified scattering and 
diminished modified scattering on the results we have already obtained. 


OncE MOopIFIED DOUBLE SCATTERING 


Referring to Fig. 12 to which the same nomenclature of entry, exit and 
intermediate beams (respectively A, C, B) applies as before we obtain a 
rough approximation to the effect of once unmodified and once modified 
doubly scattered radiation by noting that the bulk of unmodified scattering 
will occur in the region of excess single scattering for 6, represented by the 
heavy shading around A. The zones of constant shift however are now no 
longer described around the bisector point P as before but around C the 
exit direction since the shift now occurs wholly in the second scattering 
process and depends only on 62. From these considerations it is evident that 





Fig. 12. Coordinate system for double scattering, unimiodined at 
angle @, but modified at angle @.. 


the only cases of importance are those where the intermediate ray point, B, 
lies in the shaded regions of excess unmodified scattering. The zones of 
constant shift have been so spaced in the drawing that they represent equal 
increments of shift. One can see immediately that the greatest intensity 
in the spectral distribution of once modified doubly scattered radiation occurs 
closely adjacent to and about symmetrically on either side of the position 
of the modified line for single scattering, and that there will also be a con- 
siderable contribution between these two positions. Since excess unmodified 
scattering has been observed to be many times as strong as ordinary un- 
modified scattering, there is just a chance that some of the observed excess 
breadth of the modified Compton line may be attributed to this effect. This 
argument however will not apply to explain the breadth of the Compton 
line for primary scattering angles of nearly 180° because here the entire 
shaded region is included in a very small range of shifts. Singly modified 
double scattering will be absent near 6=0° because in this case when B is 
in the shaded area @ will have values for which modified scattering is prac- 
tically absent. If then the effect in question is the sole cause of the observed 
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excess breadth of the Compton line we should expect this excess breadth 
to be least for small and large angles of primary scattering, 8, and greatest 
for 8@=90°. 


Twice MopIFIED DOUBLE SCATTERING 


The effect of decreased modified scattering at small angles on the already 
derived spectral distributions of twice modified doubly scattered radiation 
occurs principally in the left hand portions of the curves. We can approxi- 
mate the facts roughly by imagining circles each about 45° in radius described 
about A and also C from which the point B is supposed to be excluded (see 
Fig. 13). For scattering angles, 0, up to 90° these excluded regions will lie 
entirely in the low shift hemisphere. We conclude then that the right hand 
halves of the curves, Figs. 4 to 6 inclusive, will remain practically unchanged 
and that the left hand halves will diminish monotonically as we pass to the 
left, vanishing almost completely at their left extremities. 





Fig. 13. Coordinate system for doubly modified double scattering taking into account the 
absence of modified scattering at small angles. 


For scattering angles @ greater than 90° as the circles of exclusion move 
out of the low shift hemisphere into the high shift hemisphere, it is evident 
that energy will reappear in the left extremity of the spectral distributions 
and the suppressed portion of our curves will displace toward the right. 
The left hand regions or regions of low shift will however always be more 
depressed than the right hand regions of great shift. Fig. 14 shows the 
spectral distributions to be expected from the preceeding discussion. 


SUMMARY OF CONCLUSIONS 
1. The ratio of double to single scattering for a spherical scatterer is 
proportional to the radius of the scatterer (neglecting absorption) and is 
given by the formula 


Doubly scatteredenergy 9 
- — = —orR(@) 
Singly scattered energy 32 
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in which the R(@) is greatest for about 6=90° and has minima at 6=0 and 
6=180°, the latter being the smaller (Fig. 3). R never differs greatly from 2.5. 
2. For MoK radiation scattered by a sphere of graphite 1 cm in diameter 
double scattering is about 14°) of single scattering (neglecting absorption). 
3. From the last it seems probable that in most experiments the effect 
of triple and higher multiplicities of scattering can be completely neglected. 
4. Once modified doubly scattered radiation may account for a slight 
broadening of the Compton line except in regions near 6=0 and 180°. 
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Fig. 14. Spectral distributions due to doubly modified double scattering corrected for 
diminished modified scattering at small angles (shaded curves). The two curves R(@) and J 
are the same quantities as those shown in Fig. 3 after application of the Breit correction and 
correction for diminished modified scattering at small angles. 


5. Twice modified doubly scattered radiation when corrected qualita- 
tively for diminished modified scattering at small angles may contribute a 
faint asymmetric line or edge in the position 


h 
= 2— (1+ cos § 6) 


me 
6. For hard radiation where the Thompson formula with the Breit cor- 
rection applies with fair accuracy to modified scattering alone, twice modified 
doubly scattered radiation contributes a spectral band of breadth 


h : 
4— cos 56 
mc 
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List OF SYMBOLS 


AQ, 6Q—Solid angle measured in steradians. 

6—The angle of single scattering. 

6,—The first angle of double scattering. 

6.—The second angle of double scattering. 

M—The multiplicity of scattering. 

Ao—The wave-length of the primary radiation. 

Ai, Ax—The wave-lengths after first and second scatterings. 

h—Planck’s constant. 

m—The mass of the electron. 

c—The velocity of light. 

A, B, C—Unit vectors representing the directions of the entry, intermediate, 
and exit rays in double scattering. (Fig. 2). 

a, ¢6—The colatitude and longitude angles of the spherical polar coordinates 
of reference in double scattering. (Fig. 2.) 

x = (Ad)/(2h/mc)—wave-length variable used for convenience in describing 
spectral distributions due to double scattering. 

EnEaEn2EwEn3EauEnysEas—components of the electric intensity in the radi- 
ation. The subscript, ”, indicates the component normal to the plane 
of one scattering angle while the subscript, a, indicates the component 
parallel to that plane. The subscript 1 refers to the initial beam and the 
first scattering angle, subscript 2 refers to the intermediate beam and the 
first scattering angle, 3 to the intermediate beam and the second scatter- 
ing angle, 4 to the final or exit beam and the second scattering angle. See 
Fig. 2. 

B—The dihedral angle between the planes of 6; and 62. (Fig. 2.) 

r,—The distance from the point of first scattering to the point of second 
scattering. 

rx—The distance from the point of second scattering to the point of observa- 
tion. 

r—The radius of the spherical scattering body. 

o—The scattering coefficient for the material of the scatterer and the radia- 
tion used. 
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ABSTRACT 

The double crystal x-ray spectrometer is distinctly different in principle from 
the single crystal spectrometer and calls therefore for a completely new design rather 
than an adaptation of single spectrometer design. One such new design is described 
in this paper. Emphasis is laid on the fact that the wave-length selected by the 
process of two successive crystal reflections in any prescribed orders depends solely 
on the dihedral angle between the crystals. The advantages of rotating both crystals 
in equal and opposite directions with respect to a plane through their axes of rota- 
tion are pointed out and a spectrometer designed to accomplish this is described. 
In this spectrometer the turning of a single shaft drives the two crystals, the spec- 
trometer as a whole and the x-ray tube, at the proper angular rates about the proper 
centers to insure that the x-ray beam will at all times remain centered on the crystal 
faces and the window of the stationary ion chamber. It is pointed out that spurious 
fluctuations may be introduced in spectral curves by the x-ray beam migrating 
across small steps in the cleavage surface reflecting the radiation. Description of 
design covers design of spectrometer and turning mechanisms, the tube housing, lead 
shields, the detecting system consisting of an ion chamber with internal grid of 
special design to minimize natural leak and a Hoffmann vacuum electrometer con- 
nected to it through a short evacuated shield. Description of operating technique 
covers optical method of accurately orienting the crystal faces, methods of aligning 
the tube, precautions to eliminate background. 


HE Double crystal x-ray spectrometer first used in this country by A. H. 

Compton was developed by Bergen Davis! and his collaborators and has 
been very thoroughly investigated on the theoretical side by Schwarzschild? 
and on the experimental side by S. K. Allison* and others.‘ Its field of applica- 
bility seems to fall into two classes. First it may be used to study the proper- 
ties of the crystals themselves (reflecting power, lattice perfections, etc.*) 
Secondly it may be used as an x-ray spectrometer for exploring spectra (satel- 
lites, fine structure line breadths, Compton effect, etc.)® In the development 
stage it was natural not to differentiate sharply in the design of double 
spectrometers for these different classes of work. It was also perfectly natural 
to adapt the general features of single crystal spectrometer design to the 


1 Bergen Davis, Phys. Rev. 17, 608 (1921), Phys. Rev. 27, 18 (1926), Phys. Rev. 32, 331 
(1928). 

? Schwarzschild, Phys. Rev. 32, 162 (1928). 

3 Allison, Phys. Rev. 34, 176 (1929), 35, 1476 (1930). 

4 Ehrenberg and Mark, Zeits. f. Physik 42, 807 (1927). 

’ Davis & Purks, Phys. Rev. 34, 181 (1929). 

* DuMond & Hoyt, Phys. Rev. 36, 799 (1930), Bearden, Phys. Rev. 36, 791 (1930), 
Richtmyer & Taylor, Phys. Rev. 36, 1044, (1930). 
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double spectrometer. The authors feel however that now that so much is 
known about the possibilities and applications of this instrument the time has 
come to design it on its own merits in such a way as best to accomplish some 
particular function. 

The double spectrometer here described was developed for the second of 
the above mentioned uses—the study of x-ray spectra. It is designed to be 
used entirely in the antiparallel positions (”, ”) (Allison’s notation)’ though 
slight modifications would permit “straddling two different orders.” 

The authors feel that much of the apparent complexity arising in connec- 
tion with the double spectrometer comes from the introduction of such artifi- 
cial frames of reference as the axis of rotation of the pivots supporting the 
crystals, or the axis of the slits which latter in the double spectrometer play 
an entirely different role from the slits of a single spectrometer. 


Tue IMPORTANT THING Is THE DIHEDRAL ANGLE 
BETWEEN THE CRYSTALS 


It is of cardinal importance to note just what is measured in using the 
double spectrometer. A beam of x-rays reflected from a cleavage face of a 


A B 


a 


Fig. 1. Symmetrical ray selected by the bicrystalline reflection. 


calcite crystal falls upon the face of a second calcite crystal and is again re- 
flected into an ionization chamber. One or both of these crystals is turned 
slightly about an axis in the plane of the reflecting face and the ionization 
currents in the ion chamber are plotted as a function of the rocking angles of 
the crystal. For any given angle setting the two successive calcite reflections 
select out of the entire x-ray spectrum an exceedingly sharply defined narrow 
region of wave-lengths. 

Precisely what angle determines the wave-length selected by the double 
x-ray reflection? Is it the angle through which one crystal rotates with re- 
spect to the slit system? Obviously not, since the horizontal slit width has 
nothing to do with the resolution of the double spectrometer. Neither do the 
pivots about which the crystals rotate constitute in any sense the system of 
reference with respect to which the rotation must be measured which deter- 
mines the selected wave-length. The authors wish to emphasize most strongly 
that the important thing is the dihedral angle formed by the two crystal faces 
themselves. Let us consider for simplicity the case of first order reflection from 
two similar calcites in the anti-parallel position (1, 1). Fig. 1 represents the 
dihedral angle formed by the two crystals. Obviously the only rays which 
can be reflected twice selectively in the first order by both crystals are the 


7 Allison & Williams, Phys. Rev. 35, 149 (1930). 
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symmetrical rays which make the same glancing angle at crystal A as at 
crystal B. A symmetrical ray, ab, is shown in Fig. 1 passing from crystal A 
to crystal B. 

Define a plane normal to the dihedral angle formed by the crystal faces 
which for brevity and in accord with other writers we will call the “horizontal 
plane,” (plane of the paper in Fig. 1). Now the ray, a, 5, need not be parallel 
to the “horizontal plane” in order that it shall make equal glancing angles at 
each crystal and therefore be reflected. Such a ray oblique to the horizontal 
plane will in accord with Davis’ terminology be said to have an angle of 
vertical divergence namely the angle a of obliquity with the horizontal plane. 
The projection of such an oblique ray will in Fig. 1 still be represented by the 
line, a, b, but the ray will in space actually suffer reflection at both crystals 
under a smaller glancing angle than rays which are not oblique. See Fig. 2. 
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Fig. 2. Horizontal and oblique symmetrical rays. 


We can assert then that the glancing angle (Bragg angle) selected by the proc- 
ess of the double selective reflection is half of the supplement of the dihedral 
angle between the crystal faces for horizontal rays and less than this for 
oblique rays. The band of wave-lengths selected then by the double reflection 
has a sharp limit on the long wave-length side determined only by the 
dihedral angle between the crystals. The band selected can be made as nar- 
row as one pleases by limiting the vertical divergence. 

Nothing in the above argument is made to depend on or refer to the pivots 
or axes of rotation of the crystals. It is of course obviously desirable that 
these axes of rotation should be parallel to the edge of the dihedral angle 
formed by the two crystal faces in order readily to permit the measurement 
of that dihedral angle and also to maintain the “horizontal ray” in a per- 
manently fixed position with respect to the frame of the apparatus as the 
crystals rotate. The pivots, however, have nothing at all to do with the se- 
lective action of the double reflection and the dihedral angle between the 
crystals has everything to do with it. If this very fundamental fact is remem- 
bered all uncertainty as to the effect of various errors of alignment, etc., will 
vanish. 


BAND Pass FILTER WIDTH 


The theoretical spectral intensity curve of the wave-length band selected 
from a uniform or “flat” spectral energy distribution by the double crystal 
reflection in the (”,”) position with perfect crystals is shown in Fig. 3. It is 
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bounded on the long wave side by a vertical discontinuity and falls off on the 
short wave side more gradually, meeting the axis of abscissae so as to have 
a base breadth AA related to the maximum vertical divergence permitted by 
the system of lead stops limiting the x-ray beam. Referring to Fig. 4 let 4, 




















Fig. 3. Spectral intensity distribution selected by geometry of bicrystalline reflection with 
crystals fixed. This curve is not to be confused with the shape of a spectral line obtained by 
plotting reflected x-ray intensity against the dihedral angle between the crystals (or half of its 
supplement). Such line shapes are due to the overlapping of the “natural” line with the above 
curve as the latter is displaced spectrally with rotation of the crystals. The equation of the 
above curve is y=k/[2(Ao—Aa)Ao}"/2; do(1 — fa®) <Aa@<Ao calculated on the assumption that 
X-ray intensity is uniformly distributed over the permitted range of vertical divergences. The 
curve has a finite area in spite of the infinite ordinate at \o. In a spectral exploration as such a 
curve progressively overlaps a spectral line the intensity will therefore always remain finite. 























Fig. 4. Dependence of glancing angle on vertical divergence. 


be the Bragg angle for a horizontally reflected ray, 6, the Bragg angle for a 
ray of vertical divergence a. Then from the geometry of the case 


sin 8, = sin 6) cos @ 


and since 
n\/2d = sin 6 


we have 


Na 
AX = Xo — Aa = Ao(l — cos a) 


Xp COS a 
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and since in practice a is restricted to small values we have as a good approx- 
imation 
AX = 3rca?. 


DISTORTION OF LINEs BY VERTICAL DIVERGENCES 

The characteristic curve of Fig. 3 which may be regarded as the “fre- 
quency response curve” of the spectrograph for any given setting of the 
crystals displaces itself bodily toward longer wave-lengths as the dihedral 
angle between the crystal faces closes up. In this exploring process if a sharp 
line is encountered at some particular wave-length the ionization chamber 
will respond abruptly on the short wave side of the line and will diminish 
more gradually as the characteristic exploring curve moves to the Jong wave 
side of the lines. For this reason the distortion of spectral lines by the 
geometrical effect of vertical divergence alone is the reverse of Fig. 3, the 
discontinuous edge being on the short wave side. A good example of such 
distortions appears in Fig. 11 reproduced from a previous article’ published by 
the authors. It is an easy matter however to so limit the vertical divergence 
that its effect in broadening spectral lines is negligble in comparison to the 
effects of crystal imperfections and natural x-ray line breadths combined. 


OBJECTIONS TO ROCKING ONLY ONE CRYSTAL 


The fact that the rays selected by the mechanism of the double reflection 
are symmetrically orientated with respect to the dihedral angle between the 
crystal faces immediately suggests that in exploring the spectrum both 
crystals should be simultaneously rotated in equal and opposite directions 
with respect to a plane through their axes of rotation. 

As a matter of fact however most double spectrometers are built so that 
only one crystal is rotated in exploring the spectrum. This necessitates rotat- 
ing the other crystal from time to time and also changing the relative position 
of the two crystals to prevent the beam from falling partly or wholly off the 
reflecting faces. As one crystal alone is rotated the selected beam must of 
necessity progress across the faces of the two crystals in order that it shall 
continue to make equal glancing angles with both crystal faces. The ob- 
jections to this procedure are: 

1. The crystals may present different reflecting powers in different parts 
of their surfaces so that spurious intensity fluctuations would be introduced 
in a spectral exploration. In particular all calcite cleavage faces have multi- 
tudinous small “steps” where the fracture changes abruptly from one plane 
to another one parallel to it. At glancing angles such steps obviously cast 
shadows in the reflected ray which will diminish its intensity in different 
degrees according to whether the beam falls on a part of the surface with high 
steps or a part in which steps are nearly absent. If the beam is obliged to 
migrate across the crystal face these steps will surely introduce undesirable 
spurious fluctuations. Even though the beam remains centered at one point 
of the crystal surface the area of crystal surface employed in reflection will 


8 DuMond and Hoyt, Phys. Rev. 36, 799 (1930). 
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change with the glancing angle and some spurious fluctuations may thus be 
introduced by cleavage steps. These must however be regarded as an ir- 
reducible minimum. 

2. If the source of radiation is not a point but has “horizontal” extension 
so as to permit reflection of one wave-length over a considerable area of the 
crystal a part of the beam which would otherwise be reflected may fall off 
one edge of the crystal introducing again spurious intensity fluctuations in a 
spectral exploration. 

3. To avoid the above danger either very long crystals or a much reduced 
reflecting area must be used, both of which are rather serious difficulties. 

4. The initial task of aligning the double spectrometer and locating a 
given spectral line is much harder in the instruments whose crystals do not 
rotate simultaneously. 

As a corollary to what has just been said it is evidently desirable also to 
maintain the focal spot of the x-ray tube or other radiation source as well as 
the ionization chamber window always so situated with respect to the 
spectrometer that on one hand the selected beam from the tube will remain 
centered on the same reflecting region of the crystal faces for all glancing 
angles and on the other hand the selected radiation will always strike the 
ion chamber window at the same point. It is easy to see that large spurious 
fluctuations could be introduced by poor alignment of the beam entering a 
narrow ion chamber window and even with a broad window the transparancy 
of the window might not be uniform over its area. 

It seemed desirable to divorce the design of the double spectrometer 
completely from the single spectrometer and produce an instrument suited 
to the requirements of a large group of problems likely to be attacked with- 
out aspiring to achieve one so universal in its applicability that it would per- 
form none of its functions well. 


DESIGN OF THE SPECTROMETER 


Fig. 5 is a general view of the spectrometer. Two circular tables A turn- 
ing on conical pivots with accurately parallel axes support the crystals. The 
interpivotal distance is 5 inches. The pivot holes are bored in a single bed 
block of steel B, situated under the front end of the instrument. Two flat 
steel levers extending from the pivots to a similar steel bed block, C, under the 
rear end of the instrument can slide laterally, opening and closing like 
scissors, on top of this block so as to give equal and opposite rotations to the 
crystal tables. The levers are given their motion by a transverse shaft D, 
with right and left hand threads passing through nuts pivoted to the ends of 
the levers. The shaft is prevented from translating axially by a thrust bear- 
ing E half way between the levers. This bearing which has no axial play at 
all has sufficient lateral play to permit the screw shaft to translate in the di- 
rection of the crystal pivots as the levers spread apart. The action of the 
entire system is such that the crystal pivots and the driving nuts at the 
opposite ends of the two levers always stand at the four corners of a perfectly 
symmetrical isosceles trapezoid whose base is the screw shaft, whose sides are 
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the levers, and whose top is the line of centers of the crystal pivots. The top 
and sides are of constant length but the base and hence the altitude varies. 
It is for this reason that the thrust bearing on the screw provides for a slight 
lateral translation of the screw shaft. This same screw shaft carries a brass 
drum F on one projecting end on which the scale of glancing angles is en- 
graved. The smallest divisions on this drum represent ten seconds of are for 
the glancing x-ray reflection angles, and are separated by about three mm. 
A very important part of the design is a helical spring G running between two 
short posts on the under side of the two levers. The tension of this spring 
holds the levers together thus completely removing all backlash in the screw, 
nuts and pivots. The crystal tables are provided with angle scales engraved 
on their peripheries to facilitate setting. The levers can be clamped to and 
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Fig. 5. Double spectrometer providing equal rotation of both crystals with respect to the line 
of pivot centers. 


unclamped from the crystal tables in the same way as the tangent screws 
on a spectrometer or engineer’s transit. The crystal tables are provided with 
removable false tops H which can rock about a horizontal axis. This is 
accomplished by two steel ball bearings lying between the crystal table disk 
and the false top in slight conical depressions provided for them and two 
opposing screws which clamp the false top to the crystal table with any de- 
sired slight inclination to the level. This adjustment is provided to permit of 
making the crystal faces accurately parallel to the axes of the pivots during 
the initial optical adjustment. A carriage K bearing the vertical crystal 
support slides in ways provided on each false top. This adjustment permits 
the crystal face to be brought in coincidence with the axis of rotation during 
the initial adjustment. The crystals are each fastened to a brass pla‘e at 
three points with a very small quantity of beeswax. The brass plate can be 
fastened to the vertical crystal support LZ with two screws passing through 
holes large enough to permit of some adjustment in the position of the 
crystals. The vertical crystal support has a large hole bored through it at the 
point corresponding to the center of the crystal face. The reason for this 
hole is to permit sighting with a telescope and Gauss eyepiece at either side of 
a small interferometer plate which takes the place of the crystal during some 
of the initial adjustments. The entire spectrometer stands on leveling screws. 
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This instrument was very inexpensive to construct. The frame was built of 
pieces of cold rolled steel of rectangular section. No hand scraped surfaces 
were found necessary nor any precise circular divided scales. The one accur- 
ate requirement was that the pivot holes in the bed block B should have their 
axes as nearly parallel as possible. The spectrometer was built by one man in 
about five and one half days but in spite of its cheapness and simple con- 
struction a large amount of surprisingly precise work can be done with it. 


THE SUBSIDIARY ROTATING MECHANISM 


As before stated the quantities measured with the double spectrometer 
are: first, the angle between the crystals® and second, the intensity of the 
twice reflected x-rays. The spectrometer just described measures the first 
quantity with high precision and reproducibility. It is however an important 
refinement to provide motion of the ionization chamber and x-ray source 
relative to the spectrometer in order to avoid the errors caused by shifting 
of the beam across the crystal faces and kindred objections already men- 
tioned. It should be understood however that the requirement of precision 
in the motion of source and ion chamber is not as exacting as it is in the angu- 
lar motion of one crystal relative to the other since the precise wave-length 
selected by the double reflection does not depend on the position of source or 
ion chamber so long as these are so situated that some radiation can suffer the 
double reflection. The wave-length selected by the double reflection depends 
wholly and solely on the angle between the crystals and it is this quantity 
which must be and is measured precisely and reproducibly on the drum of the 
spectrometer. 

The motion of x-ray source and ion chamber relative to the spectrometer 
in the instrument here described is accomplished in such a way that the ion 
chamber remains stationary relative to the room. This means that relative 
to the room the spectrometer rotates about a center in line with the axis of 
rotation of the second crystal and at twice the rate of the second crystal. At 
the same time relative to the spectrometer the x-ray source rotates about a 
center in line with the axis of rotation of the first crystal at twice its rate of 
rotation. The scheme of holding the ion chamber stationary was adopted in 
order to permit the use of a short unarticulated vacuum sleeve surrounding a 
short lead wire of low and perfectly constant capacity connecting the ion 
collector to the electrometer. The ion chamber could probably be permitted 
to move relative to the room by using a vacuum tube amplifier such as has 
recently been developed for The General Electric Company in a shielded 
housing attached to the ion chamber to replace the electrometer. There is 
some question as to whether the enormous sensitivity, stability and repro- 
ducibility attainable with the Hoffmann type vacuum electrometer can be 
matched with the amplifier tube, and work is now being started to investigate 
the possibilities of this method in connection with a vacuum double crystal 
spectrometer now in process of design at this Institute. 


® The angles for which the drum is calculated are half the supplement of the dihedral 
angle between the crystals to facilitate their direct interpretation as glancing angles. 
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The above mentioned relative motions are obtained by means of two ro- 
tating false table tops of ply wood best shown in Figs. 6, 7. Ply wood was 
used instead of metal because the presence of the x-ray tube fed with high 
voltage made an insulating material convenient. The ply wood has proved to 
be quite free from warping and entirely satisfactory. The upper table top, F, 





Fig. 6. General view of spectrometer, rotating table tops, rotating mechanism, and x-ray 
tube housing. 
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Fig. 7. Schematic diagram of table tops. 


Figs. 6 and 7, carries the source of x-rays S. It rotates on a pivot centered at 
A consisting of a disk of ply wood 5 inches in diameter fixed to the second 
false table top G. This second table top in turn rotates on a stationary pivot 
centered at B and attached to the stationary table top. The centers A and B 
are separated by the same distance as the interpivotal distance on the spectro- 
meter. The spectrometer rests entirely on the second table top G, the position 
of its three feet being at K, L and M. (The reason for the large diameter 
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pivot A is here obvious.) The pivots of the crystals are directly above and 
in line with the table top pivots A and B. The table tops are given the re- 
quired angular rate of rotation by a second right and left hand screw with 
nuts and thrust bearing almost identical with the one on the spectrometer 
proper and having the same thread pitch. The two screws are connected 
through a simple type of universal joint transmission to gears having a ratio 
of two to one. The post supporting the gears is fastened to table top G and 
hence is stationary relative to the spectrometer. The entire mechanism is 
rotated by turning a shaft set at right angles to the screw shafts for accessi- 
bility the former driving the latter through bevel gears. This shaft terminates 
in a hollow brass tube T slotted axially on two diametrically opposite sides. 
A brass ball on the end of a long insulating handle can be inserted in this tube, 
and a rotary motion given to the tube by means of a diametral steel pin 
through the brass ball which engages in the slots of the tube. (The universal 
joints mentioned above are all of this simple ball and pin and slotted cylinder 
type and are very satisfactory.) 

The kinematics of the table top drive is identical to that of the spectro- 
meter. It is most easily grasped by taking as a reference, not the stationary 
table, but the false top G on which the spectrometer stands. The thrust bear- 
ing is attached to this top, and on turning the screw one nut causes the x-ray 
source on table F to rotate around the center A, relative to the spectrometer 
standing on G, at twice the rate of the crystal at A, while the other nut causes 
the table proper, and with it the room and the ion chamber, to rotate relative 
to the spectrometer standing on G around the center B, at twice the rate of 
the crystal at B. The same provision for lateral displacement of the shaft in 
the thrust bearing is provided for here as on the spectrometer. The whole 
device is of course open to the objection that it moves the source and ion 
chamber not accurately at twice the angular rate of the crystals, but rather at 
such a rate that the sine of one angle is twice the sine of the other angle. This 
however introduces far too small an error to cause trouble over the small 
rotations ordinarily required. 

The lever arms of the spectrometer have each a range of motion of a little 
over one degree on each side of their position of mutual parallelism. The 
tables have each a range of two degrees on each side of their central position. 
A great deal of work can be done with even so small a range as this but slight 
modifications would suffice to increase the range. It is our practise to pick 
some definite angle setting near the midpoint of the range of glancing angles 
to be studied and then orient the crystals and other parts so that this angle 
corresponds to the zero position of the drum F (Fig. 5) with the levers in their 
central position. Details of this operation are described below. 


MINOR SPECTROMETER DETAILS 


The tube housing shown in Fig. 6 consists of a micarta tube large enough 
to contain the x-ray tube covered with a wrapping of two layers of the leaded 
rubber 1/16” thick such as is used by Réntgenologists for protective aprons. 
This material which can be obtained from surgical supply houses has the 
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advantage of being opaque to x-rays, flexible, and a good insulator. Its uses 
are manifest to anyone who has worked in x-rays. Since it is an insulator one 
can construct small and therefore light boxes of it to contain x-ray tubes. 

Fig. 6 was taken with a copper target tube in the lead rubber housing 
provided with a projecting goose-neck tube with a charcoal trap dipping into 
a liquid air container. The latter is plainly visable standing on table F. An 
investigation over a period of two months!° was made with this copper target 
tube which had a mica window sealed on one side with DeKhotinsky cement. 
The liquid air and charcoal were found quite sufficient to maintain an ex- 
cellent x-ray vacuum throughout the entire period. 

The rotation required of the x-ray tube does not prevent the use of mer- 
cury pumps if necessary. These are installed under the table and are hung 
from table top F through holes bored in G and in the table top proper 
sufficiently large to permit the complete latitude of rotation. The pumps can 
be connected directly to the tube with an all-glass connection since no relative 
motion occurs between the two. 

In order to insure smooth and easy operation of the table tops the bulk of 
the weight of the x-ray tube, housing, pumps and spectrometer is relieved by 
a simple system of lifting ropes going up to levers and counter weights hang- 
ing from the ceiling. 

To give greater accessibility the entire table bearing all the parts so far 
described can be unclamped and moved away from the concrete wall shelf 
on which the detecting system consisting of electrometer and ion chamber 
stands. 


DESIGN OF THE INTENSITY MEASURING SYSTEM 


The essential requirements are that the x-ray intensity measuring system 
shall have high sensitivity coupled with great stability and reproducibility. 
It is also convenient to have as large an x-ray window in the ion chamber as 
possible to make the task of alignment less arduous. Fig. 8 shows a cross 
section of the entire ion chamber-electrometer system. The chamber is a large 
cylindrical brass can A nickel plated inside and out. A disk shaped cover B 
fastens on to the end of the can with six screws, but is electrically insulated 
from the can with a gasket and insulating washers under the screws. The 
cover is made gas tight by pressing a fillet of red universal laboratory wax 
into the groove C provided for it. The 1” X3” x-ray window in this cover is 
covered with a piece of paraffin impregnated Balsa wood—an extremely 
light yet stiff wood used in airplane model construction—with the grain run- 
ning across the short dimension of the opening. The wood is attached to the 
cover with beeswax. The inner surface of the wood is coated with aluminum 
beryllium alloy 0.0002 inches thick making contact with the metal cover. The 
purpose of this coating is completely to define the potential all over the win- 
dow so as to insure that the ions formed in this vicinity will not fail to be 
collected and also to prevent electrometer disturbances due to charges on 
the window inducing fluctuating charges on the collector. 


10 Dumond and Hoyt, Phys. Rev. 36, 799-809 (1930). 
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A grid D of No. 36 nickel wire wound on and spot welded to a light nickel 
wire frame extends back from the window and surrounds the collector. This 
grid is supported on the brass cover B and by means of an external wire con- 
nected to the cover is brought to the potential of about 135 volts above 
ground, furnished by radio B batteries. The rest of the can which stands in 
contact with the electrometer housing is grounded. 

The collector is led out of the can through a quartz insulator and as the 
can is itself grounded no guard ring is necessary. The collector of nickel wire 
can be fastened with a set screw into a brass cap on the top of the quartz 
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Fig. 8. Complete cross section of x-ray detecting system. 





insulator. This collector is put in place through the x-ray tube window just 
before sealing on the Balsa wood cover. A small piece of fine mesh nickel 
screen E is welded over the inside opening of the gas connection F. This 
saves having to shield the glass tubes leading into the ion chamber by the 
usual method of covering them with tin foil. 

The purpose of the internal grid is to minimize the fluctuations in the 
natural leak of the system caused by alpha-particles from the walls of the ion 
chamber. The space of about 5 cm from grid to chamber walls filled with 
about an atmosphere of methyl bromide accomplishes this very well. 

The gas outlet from the ion chamber leads to a mercury manometer on 
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which any change in the methyl bromide pressure can be read. There is a 
glass cock for filling and a drying tube in constant connection with the 
chamber. It is important to keep’the quartz insulator dry. 

A fine No. 40 wire passing through the quartz insulator connects the 
collector to the highly insulated suspension of the Hoffman type electrometer. 
The connection is made through a short evacuated sleeve G which is made 
vacuum tight against the top flange of the electrometer with a fillet 7 of red 
universal wax. 

Save for the flat flange just mentioned and for the use of universal wax 
instead of rubber gaskets to make the cylindrical case of the electrometer 
vacuum tight there are no changes in the design of the Hoffmann type 
vacuum electrometer. This remarkably satisfactory instrument has been 
completely described by its designer" and will not be given space in this 
article. 

SENSITIVITIES AND REPRODUCIBILITIES OBTAINED 

The spontaneous fluctuations in readings with this x-ray measuring system 
correspond to about 5X10-'* amperes. These fluctuations are the limiting 
feature at present on the smallest currents readable. They are only due in 
part to alpha-particles from the walls of the chamber. Their source is being 
studied and it is hoped that they can still be very greatly reduced. With these 
fluctuations at the present value currents of 6X10- amperes can be read 
with a reproducibility of one percent in a single reading of the order of a 
minute’s duration. Currents of 10- amperes can be read with fair accuracy 
by taking several somewhat longer readings and averaging. We have suc- 
ceeded in detecting the satellite Ka3, of copper at only 1.75 KV above its 
critical exitation voltage with this outfit. This satellite had spectral ordinates 
protruding above the background only one four hundred and fiftieth as high as 
the maximum ordinate of Ka, its parent line six X.U. away. 

The electrometer has a very stable zero that will remain constant within a 
few millimeters for months. The instrument as now used has a rather stiff 
suspension. Mr. Julius Pearson of this institute has succeeded recently in 
producing Wollaston wire suspensions for the Hoffmann electrometer with 
much longer period than the one now in use. The present rather low volt 
sensitivity could then be increased to about 19000 mm per volt. The capacity 
of the whole insulated system of electrometer and ion chamber is about 12 
electrostatic units. We therefore do not believe the possibilities of the above 
described system are exhausted with the results here reported. 


TECHNIQUE OF OPERATION 


Method of orientating the crystals initially. 

The procedure about to be described accomplishes in practice the ful- 
fillment of the following required conditions. (These conditions are not all 
equally exacting.) 

1. The crystal reflecting planes must be parallel to the axes of rotation of 
the crystal table pivots. 


11 G. Hoffmann, Ann. d. Physik 52, [7] 665-708 (1917). 
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2. The reflecting surfaces of each crystal must coincide with the axis 
of rotation of its pivot. 

3. The two crystal reflecting planes must make equal angles with the 
plane through the axes of the pivots. 

4. The dihedral angle between the two crystal planes must be set for 
some known value corresponding to the “zero” drum position of the screw 
and lever crystal turning mechanism to facilitate “finding” a spectral line. 

The first condition is quite important since the accuracy with which the 
dihedral angle between the crystal faces can be determined and hence the 
accuracy to which the wave-length scale is known depends on this. Also if 
condition No. 1 is not fulfilled the so called “horizontal” plane i.e. the plane 
normal to the dihedral angle formed by the two crystals may be quite oblique 
to the pivot axes and the general line of the x-ray apertures causing an un- 
suspected vertical divergence of the rays from the “horizontal” plane. Re- 
quirement No. 1 evidently becomes more rigid due to both of these last 
named effects as the crystal planes approach parallelism. 

The second condition cannot be fulfilled very rigidly since the plane of 
the face of the crystal is ill defined on account of the unavoidable cleavage 
steps. Its purpose is to minimize the gliding of the selected beam across the 
crystal face as the crystal rotates. 

The third condition merely insures that the radiation from the first to the 
second crystal shall pass parallel to the line of centers of the two pivots thus 
facilitating the initial alignment of the apparatus. It need only be fulfilled to 
within a few minutes of arc. 

The fourth condition is merely a convenience. We have however invari- 
ably succeeded by the technique about to be described in locating the crystals 
so closely that it was not necessary to hunt for a known spectral line over a 
region much greater than its own observed breadth. 

A small optical spectrometer with collimator removed and with a Gauss 
eyepiece in the telescope is used. Calcite cleavage faces as Allison has noted 
are poor reflectors of visible light at normal incidence. It is however, possible 
to obtain a sufficiently intense image of the cross hairs for visibility if one 
employs a strong but small incandescent light, well shielded from directly 
striking the eye to illuminate the Gauss eyepiece. A small black disk of 
paper, say 0.5 mm in diameter, should be carefully attached with a little 
shellac to the intersection of the cross hairs. The period from a set of gum- 
med letters such as are used to stick on drawings makes a good disk. This 
makes a spot large enough to see in the faint reflected image. Considerable 
patience is required in hunting for the faint image but the telescope and Gauss 
eyepiece furnish such an accurate method of alignment that it is worth 
while spending an hour or so, especially as it need not be done often. 

A glancing angle @ is first chosen as a setting for the crystals when the 
angle drum is set at zero and the levers and drum are then adjusted to that 
zero position. The optical spectrometer is then set on a bench in front of the 
x-ray spectrometer in such a position that in sighting the telescope alternately 
first at crystal one and then at crystal two the optical telescope must turn 
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through approximately the angle 26). The pivot axes of both spectrometers 
are then rendered accurately vertical with the leveling screws by taking 
reversed readings with a good level. A pair of light plumb lines of white silk 
thread are then adjusted so that they define a plane the perpendicular bisector 
of the plane through the two pivot axes. This can be done very precisely with 
a trammel which is engaged at one end in the conical center holes on which the 
crystal pivots were turned and is used as a “feeler” to strike four arcs in 
space centered each from a pivot and just making contact on either side 
with the plumb lines. The geometry is the simple Euclidian construction 
for a perpendicular bisector. The telescope is now adjusted to sight accurately 
down the line defined by the two plumb lines. The index of the divided circle 
on the optical spectrometer is then set on zero and clamped and the two 
spectrometers are rechecked for level and sealed to the bench with shellac or 
wax. The next step is to render the line of collimation of the telescope 
accurately normal to the pivot axes of the x-ray spectrometer. (Incidently 
when this is done the line of collimation of the telescope will be horizontal.) 
A small plane parallel glass interferometer plate is stuck on with wax over 
the hole in the vertical brass crystal support and trial reflections are made 
from both sides of this plate turning the crystal table through 180° each time 
and correcting half the error at the plate and half in the level adjustment of 
the telescope (not the spectrometer of course) until the Gauss cross hairs 
coincide with their reflected image equally well with the plate facing either 
way. This condition must be fulfilled for both crystal pivots. The telescope 
collimation line is now accurately normal to both the x-ray and optical 
spectrometer pivot axes and its level must not be subsequently disturbed. 
The plane parallel glass plate is now removed and the calcites ready mounted 
on their brass plates are attached to the vertical crystal supports with screws. 
Care should be taken to locate them symmetrically over their axes of rotation 
at the same height above the turntables. At zero reading on the divided circle 
of the optical spectrometer the telescope is now normal to the plane through 
the axes of the x-ray spectrometer pivots. Turn the telescope through the 
angle @, to right or left of this point as read on the divided circle of the optical 
spectrometer and sight at the crystal face. The crystal table is then un- 
clamped and orientated and the false top tilted until the cross hair image 
normally reflected from the calcite coincides with the cross hairs. The crystal 
reflecting face is brought into coincidence with its pivot axis by screwing the 
crystal support carriage diametrically across until a tangential view of the 
face with the crystal rotated to either side the stationary line of sight shows 
no displacement. This should be done after rendering the crystals vertical 
but before finally clamping them to the levers at the correct horizontal 
orientation defined by the telescope. 

The entire procedure requires three or four hours and need only be per- 
formed once for each pair of crystals used unless some accident disturbs 
them. The divided circles on the crystal turn tables should be carefully 
read and recorded.” This gives a fiducial reference from which the proper 


2 A scale capable of being read to within a few minutes of arc is sufficiently accurate. 
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settings for other glancing angles can be made without having further re- 
course to the Gauss eyepiece reflection method. 


Alignment of other parts 


It is necessary also to align the source of x-rays in the proper position 
with respect to the spectrometer. Furthermore the entire system of source 
and spectrometer must be aligned with respect to the ion chamber so that 
the rays will enter its window. We have found that this can all be done 
very conveniently after the crystals have been aligned by cutting a templet 
out of Bristol board such as is shown in Fig. (9). The edge segments A, B, 
of this templet are fitted against the crystals and the edges C, D, then point 
respectively in the direction of the source and the ion chamber window when 
their locations are correct. If the templet is carefully laid out with a precision 
vernier protractor and accurately cut with a razor blade held against a steel 





Fig. 9. Cardboard templet used in aligning table tops and x-ray tube. 


straight edge the method is accurate to within a few minutes of arc. It is 
helpful to lay out the normals to the crystal faces on the cardboard. When 
the cardboard is properly orientated against the crystal faces the reflected 
image of these normals should line up straight with the normals themselves. 

Needless to say the cardboard templet is laid out for the glancing angle 
6) corresponding to the mid value of the spectral range chosen for the par- 
ticular work to be done and the above described alignment procedure must 
be carried out with the spectrometer and the ply wood table tops set at the 
midpoint of their travel, i.e. at the point where the pivots of the nuts on the 
right and left hand screw shafts are separated by the same distance as the 
two crystal pivots. 


ELIMINATION OF BACKGROUND 


The detection and accuracy of measurement of faint lines is limited by 
the intensity of the x-ray background. Also the shapes and positions of lines 
are distorted by background structure. It is therefore essential to consider 
means of reducing background. The causes of background are: 
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1. Selective reflection of continuous or “white” radiation from the source. 

2. Natural leak of ion chamber and electrometer system. 

3. Amorphous or non-selective crystal scattering and fluorescent radia- 
tion excited in the crystals. 

4. Leakage of x-rays without suffering double reflection in the prescribed 
orders. 


It is not always desirable to eliminate all of the above factors. For example 
if one is interested in studying either the continuous “white” radiation or 
diffuse lines partaking of the nature of a continuous distribution No. 1 
should obviously be retained. Methods of controlling the above sources of 
background are in the above order as follows. 


1. The intensity of white radiation adjacent to a spectral line as measured 
by the ionization chamber relative to the intensity when the crystals are set 
for the peak of a line can be very widely controlled by controlling the cross- 
fire or vertical divergence of the x-ray beam. As long as the band of wave- 
lengths selected by the bicrystal reflection is of larger order of magnitude 
than the natural line breadth, diminution of this band width reduces the 
background ordinates faster than the line ordinates. The vertical divergence 
angle is controlled by the height of defining lead openings and their distance 
apart. In the present instrument one such opening provided with adjustable 
jaws is arranged in the heavy lead wall dividing the spectrometer and x-ray 
source assembly from the electrometer and ion chamber. 


2. Precautions for minimizing “natural leak” by means of the internal 
grid in the ion chamber have already been described as have also the use of 
quartz insulators kept dry with phosphorous pentoxide and the use of 
vacuum shields and a vacuum electrometer. 


3. It is well known that scattering of x-rays is of two types namely 
modified and unmodified scattering. It is also well established that selective 
crystal reflection is solely of the unmodified type since this is the only type 
in which the incident and reflected beams are coherent. As regards modified 
scattering however the two calcite crystals might as well be two pieces of 
wood or graphite. All wave-lengths of radiation from the x-ray tube can be 
scattered into the ion chamber by the Compton modified scattering process 
independent of the angle setting of the crystals. The ratio of modified 
amorphous scattering to unmodified selective scattering increases with the 
hardness of the radiation. It is evident therefore that a great advantage in 
background reduction will be gained by filtering out all radiations from the 
source harder than those which are in the range being studied. This is 
especially true when the x-ray tube is being operated at voltages consider- 
able higher than the critical voltage necessary to excite the radiation being 
studied. Thus in the study of the satellite Ka3;, of copper (parent, Ka) 
we have found that a filter made of a piece of nickel foil rolled out to 0.0005” 
thickness was extremely effective in reducing the ratio of background or- 
dinates to line ordinates. 

The same procedure evidently reduces the background due to char- 
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acteristic fluorescent radiations excited in the crystals themselves by the 
harder regions of the continuous spectrum from the x-ray tube. 

4. Finally a judicious disposition of lead shields will be found to reduce 
the background very materially. A vertical lead knife edge should be placed 
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Fig. 10. Scheme of lead shielding. 


normal to the reflecting face of each crystal, close enough to the face to 
prevent any radiation from reaching the ion chamber without suffering the 
double reflection. These knife edges if placed too close will cut down the 
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Fig. 11. Comparative curves of a;,2 before and after taking the precautions outlined in text 
for reducing relative background. 


available intensity. A little experimentation is the best guide in this matter. 
A shield with a rectangular window of the proper dimensions placed sym- 
metrically between both crystals is helpful as is also an extensive lead shield 
near the first crystal between it and the tube. Most effective of all is the 
use of a lead snout of rectangular cross section which hangs in front of the 
adjustable lead window in the lead partition between the spectrometer and 
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the detecting system and which projects forward to the second crystal. 
See X Fig. 6. The disposition of this snout if properly chosen will greatly 
reduce the solid angle from which stray radiation can enter the ionization 
chamber. Fig. 10 shows the system of lead shields. 

The degree to which relative background can be reduced is shown in 
Fig. 11 by a comparison of the spectral curves of Ka; of copper before and 





Fig. 12. General view of apparatus with lead shield between spectrometer and ionization 
chamber removed. 


after taking the various precautions discussed above for reducing back- 
ground. The ordinates of the curve with the low background have been 
plotted to a scale chosen so as to make both curves coincide at the maximum 
value of Ka;. The relative reduction of background is very evident. 

Fig. 12 is a general view of the apparatus with the lead shield removed 
from between ionization chamber and spectrometer. 
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THE VALUE OF e/m BY DEFLECTION EXPERIMENTS 


By G. E. UHLENBECK AND L. A. YouNG 
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(Received November 10, 1930) 


ABSTRACT 


To clarify the discussion as to whether or not a quantum mechanical treatment 
of the motion of electrons in magnetic fields will lead to a formula for e/m, differing 
from that obtained by classical electro-dynamics the following problem was solved. 
A uniform magnetic field in the z-direction exists in the half-space x>0. A plane 
monochromatic de Broglie wave, travelling in the positive x-direction representing 
electrons of arbitrary energy, impinges normally on the plane x =0. Solutions of the 
wave equation were found fulfilling appropriate boundary conditions at the plane 
x=0. Currents are calculated quantum mechanically and compared with the cor- 
responding classical expressions. It was found that for electrons possessing energies 
of the order of magnitude used in deflection experiments, no observable deviations 
from classical results are predicted. Another quantum mechanical effect is diffraction 
at slits. Simple approximate calculations show that this effect can produce a fractional 
error in e/m of the order of the de Broglie wave-length divided by the slit width. 
These results are opposite to the conclusions reached by Page. We may remark that 
he solved a problem of “stationary states” which does not represent the actual ex- 
periments. 


N A recent issue of the Physical Review Professor Leigh Page! tried to 

explain the difference in the value of e/m obtained from deflection experi- 
ments and from spectroscopic evidence on the basis of quantum mechanics. 
This was later criticized by Eckart,? first, he doubted certain approximate 
calculations of Page which may, however, be shown to be correct ;* second, 
he referred to articles of Kennard‘ and Darwin’ which prove quite con- 
clusively that a wave packet representing an electron in a uniform magnetic 
field moves as in the classical theory. To this we may remark that in the 
deflection method the measurements are not made with a single particle but 
are statistical since we work with beams of electrons. However the work of 
Page does not seem conclusive because he considers the solutions of the wave 
equation corresponding to the completely quantized motion of isolated 
electrons in a magnetic field of infinite extension; whereas in the deflection 
experiments we have to deal with a problem of “streaming.” 

§2. In trying to decide the question as to whether or not quantum 
mechanics will lead to results appreciably different from those of the classical 
theory we have considered the following problem. A homogeneous magnetic 
field in the x-direction fills the half-space x >0. Free electrons moving in the 


1L. Page, Phys. Rev. 36, 444 (1930). 

? C. Eckart, Phys. Rev. 36, 1014 (1930). 

’L. Page, Phys. Rev. 36, 1418 (1930). 

4 E. H. Kennard, Zeits. f Physik 44, 347 (1927). 
§C. G. Darwin, Proc. Royal Soc. A117, 258 (1927). 
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x-direction impinge on the plane x=0 and penetrate the magnetic field. 
Classically they will describe half-circles with centers on the plane x=0 
and radii 


a 


a (1) 


where p represents the momentum of the particles. We find easily that in the 
magnetic field the components of the current density S are 

: 2Ix 

S,=0 S,= ccmmnerenins (2) 

(r? 3 x2)" 2 
Tr . . 7 

and the total current /,S,dx=2Jr where J is the current density of the in- 
coming electrons. Quantum mechanically we have for x >0 the wave equation 


Vy + / : ~) (=) + #5 +(=)! 0 3 
4 Ae —{ x— — y—]} — { — } (2 2 —) \y,= ’ 
ie Ar "ay , Ox Xr . , d f mn (3) 


where A, is the Laplace operator for two dimensions and \=h/p is the de. 
Broglie wave-length of the incident electrons. The incoming electrons are 
represented by a wave function ¥;=A exp (27i/Ax) and the outgoing elec- 
trons by ¥,=B exp (—27i/Ax) so that for x <0 the total wave function 
Y2=Yity.. We must seek solutions of (3) which fulfill the boundary con- 
ditions at x=0: 


Vi = v2 (a) 
(mvz) Wi = (mvz)a2 (b) (4) 
(mvy) Wi = (mry) oe (c) 
where 
h @ 7 
mv, = —- —+—A,z 
2mri OX ¢ 
and 
h a € 
My = — — + —A,. 


2mri Oy C 


In our case the components of the vector potential A, and A, are given by 
A,=-—}Hy and A,y= Hx. Making use of results given in a recent paper 
by Landau® we find that 


Vi = o(xerizulr (5) 


is a solution of (3) if d(x) satisfies the equation 


2r\? x 
(x) + (=) (1 — ~) o(x) = 0. (6) 


*L. Landau, Zeits. f. Physik 64, 629 (1930). 
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We may therefore write our boundary conditions (4) (conditions 4c is iden- 
tically fulfilled) 


¢(0) =A+B (a) 
2ri 
¢’(0) = +.” — B) (b) (7) 
o(~) = 0 (d) 


It will be shown later that the condition ¢(*)=0 will, for a given X 


determine $’(0)/@(0) so that we may write for convenience 
Ae (8) 
2no(0) ' 


where g(A) is a real quantity. 
The first and second of the relations (7) then give 


g(A + B) = i(A — B) 


or 


B = ——A = Aer’. (9) 


The phase of the outgoing electrons therefore differs from that of the in- 
coming electrons by 77, where 


1 — 
cos TT = aa (10) 
q? 
Thus 
(0) = A(1 + e***), (11) 


§3. The current density is in general given by 


j dy _dy a 
S.=—— (— - i) -—AW 


4rim\ dx Ox mc 


with a similar expression for S,. Substituting equation (5) and the vector 
potentials we find immediately 


: 2. 
(S,)1 = 0 (S,)o =e —(BB _ A?*) = 0. (12) 
m 
This latter expression is zero due to the fact that A and B are equal in 
absolute value from (9). The intensities of the incoming and outgoing waves 
are equal. Further 


(8 « 
(Sy): = —— —xd$ (S,)2 = 0. (13) 
m Xr 
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At x=0 (S,)1=(Sy)2=0. Continuity of the currents is therefore a conse- 
quence of the boundary conditions (4). Finally we can calculate the total 
current 
72 ch a 
(Sy)\dx = -—-— xogd x. 
J 9 Arm J 9 

By integrating partially and using the differential equation (6) for ¢ and the 
fact that (©) =0 it is found that 


oz lL 1 - . 
f xogde = —F°(1 + 42)9(0)4(0). 
2 


Substituting (11) and noting that J=ehA*/md we find that 


/ 


[>] 
| (Sy)idx = 2Ir 
0 


which is just the classical value. 

§4. To get a clearer idea of (S,):, a quantity which would be measured 
in a deflection experiment, we must study the solutions of (6) in more detail. 
Introducing & = 2x(7/Ar)!”, v-+3 =7r/d (6) becomes 


d* 1 1 
+(>+ -+#)o=0. 
dé? 2 4 


This equation is treated extensively in Chapter 16, page 341 of Whittaker 
and Watson “Modern Analysis.” In their netation we obtain for a solution 
vanishing at&é=@: 











r ‘vy — 1) (vy — 1)(v — 2)(v — 3) 
o(é) = aD,(t) = ae rp) i — ~— wow oe id aa ey Benet 


for large & (a is a constant to be determined from ¢(0). Ifvis an integer the 
series breaks off and we obtain the well-known harmonic oscillator eigen- 
functions. This solution “connects” with the following solution for small &: 


( 7t\.. 

I pa DQvi2+i/4 

? £2 
= s 


$(é) = a} 7, 1 “', "*Myasarsan( =) 
i 7 ] 


1 j 

3 tide Qei2+ila 
(-)2 , 

+ caaalenibné 2h paa(=) 


2 ) 








where 


Mx.»(n) = nil2tme-a/2 J 1 +- _ ead , 





Ltm—k (g+tm—k)G+tm—k) , \ 


(oT 12m+1)" * ~ 21(2m + 1)(2m + 2) 
































VALUE OF e/m BY DEFLECTION EXPERIMENTS 
From this solution one can show easily that 
v 
ri—+1 
2 vO 
gq = 2(2v + 1)-1/2 tan — (14) 


r y 1 2 
(= +=) 


vr 
= tan ry forvy>1 








thus 
1 — ¢’ 
cos #r = ———— = COS px 
1+ q° 
or 
UT = TY. (15) 


Thus the phase shift between the incoming and outgoing electrons is com- 
pletely determined in terms of v and therefore in terms of the energy of the 
incident electrons. It may be worth while to remark that v is very simply 
related to the number WN of de Broglie wave-lengths on the circumference of 
the classical half circle 


y= N —} or wr = (N — 3) 
or since J is usually of the order of 10-* cm, v is of the order of magnitude 
10§ to 10'°. For quantum numbers of this order of magnitude one would 
hardly expect to find the differences between classical and quantum mechani- 
cal predictions to be measurable. 

§5. To add weight to this belief we may consider the character of the 
function ¢(x) in more detail. This is best done by considering the Wentzel- 
Brillouin-Kramers’ approximation. For x>r (x) behaves as an increasing 
or decreasing exponential and to satisfy the condition ¢(@)=0 we must 
choose the latter. We obtain 

{ 2s [* 


(x) : *— yirarl 
i OG (x? — ?)"*dx>. 
° (x? — r?)tl4 , Xr , iit f 


r 


This “connects” according to Kramers with the oscillatory solution for 
x<r 


(2) 2a 1 fe 2112 =t 
SS ee —— OF 4 wv Sa = oe? 
2) = pin Var J, aa’ 

We see immediately that this approximate solution determines g(A) =\¢’(0) / 


27(0) and we find for it again 


7 G. Wentzel, Zeits. f. Physik 38, 518 (1926); L. Brillouin, C. R., Juli, 1926; H. A. Kramers, 
Zeits. f. Physik 39, 828 (1926); A. Zwaan, Utrecht Dissertation, 1929; L. A. Young and G. E. 
Uhlenbeck, Phys. Rev. 36,,1154 (1930). 
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VT 
= tan — 
q ») 


The constant a has to be determined from the value (0). We find from (11) 


A(1 + e**)r' 2 
¢= . 
TV 
cos 
2 


Calculation of the current gives then for x <r 


47x a ft Bia T \ 
(Sy)1 = — - cos* ) " i) i = x*)'/dy — 1 ( . (16) 
rv, 


(r2— x7)! 2 


The current, therefore, is “oscillatory” with an amplitude just twice the 
classical current. The distance between the maxima is, for small x, just 
\/2 but increases as x approaches r. The distance A of the last maximum from 


x=r is given by 
2dr 7 ” ” » Tv 
-~ (r? — x*)'/2dy = — 
Ar r—A 4 


or 
> a . 
A= = 2/3877 /3 | (17) 
For x>r the current is given by 
(S;); = = exp bis = fw — r*)! ‘dx ; (18) 
) (Qe — rye" Loxrd, 


This current falls off very rapidly with increasing x, in fact in going a dis- 
tance nA beyond x=, the ratio of the current to the current at x=r—A is 
given by 
1 
——— exp (— 2n3/?), (19) 
2(n)'/2 
This ratio is equal to 0.001 for m=2.1. 
§6. These last results are perhaps doubtful since in the neighborhood of 
x=r the approximation from which they were derived becomes invalid. 
We can, however, derive them more rigorously by observing that near x =7, 


x? x 
1—-— — 22(1 -=), 
r? r 


Instead of (6) we have then 


2r\? x 
o”’(x) + 2(=) (1 7 =) o(x) = 0. 
oN r 








——— 
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This equation has been thoroughly discussed by Kramers (reference 7) and 
the solution which “connects” properly with the exponential and oscillatory 
solutions of (6) is found immediately to be 


Ho ob(2Y" (0-5) : 


w(z) has been tabulated by van der Held (see Kramers, reference 7) and we 
find from his table 


A = }n-2/3p2/871/3 (17’) 


and in traveling a distance 2.5A the current drops to 0.001 of its maximum 
value. 

§7. Returning now to the value of e/m we must distinguish between 
two quantum mechanical effects. The first is the one discussed above. If 
we observed the maximum current we would find it displaced a distance 
A from the classical position. If we calculated e/m from 

e cv 1 


m 7 H r 


we would obtain too large a value. This is in the right direction to explain 
the discrepancy between the deflection and spectroscopic values of e/m 
but quantitatively the relative error due to this is of the order (A/r)?/* which 
is about 10-*. The observed relative discrepancy is 4107. 

A second quantum mechanical effect in an actual deflection experiment 
would be diffraction at the entrance slit. A rough calculation based on the 
uncertainty principle shows that this can give a relative error in e/m no 
larger than \/d where d is the slit-width. For actual cases this is of the order 
10-*. We conclude therefore that the quantum mechanical analysis of the 
deflection experiments cannot explain the observed discrepancy. One might 
still think of relativistic or spin effects as a possible explanation of the e/m 
paradox, but we have a feeling that the resulting corrections will also be 
extremely small. 


Note added in proof: We regret that, in preparing this manuscript we 
overlooked the article of Charlotte T. Perry and E. L. Chaffee (Phys. Rev. 
36, 904; 1930) which seems to show quite conclusively that the discrepancy 
between the two values of e/m is due to an error in the experimental deter- 
mination of the velocity of the electrons in the deflection experiments. 
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RELATIVISTIC WAVE MECHANICS OF ELECTRONS 
DEFLECTED BY A MAGNETIC FIELD 


By Mitton S. PLEsset1* 
SLOANE Puysics LABORATORY, YALE UNIVERSITY 


(Received November 3, 1930) 


ABSTRACT 


It is shown that the relativistic wave equation for electrons in a uniform magnetic 
field leads to the same wave function as that already deduced by Page from the non- 
relativistic equation. As in the latter case the motion at right angles to the field is 
quantized. 

An expression is found for the current density from the relativistic wave equation. 
The relativistic expression differs from the non-relativistic only by a constant factor 
which does not affect the calculation of the mean radii of curvature of the electron 
current. 

Hence, for the relativistic case, as for the non-relativistic, the mean radius of 
curvature is less than that expected on the classical theory. It follows that the classi- 
cal relativistic relation between e/u and the mean radius of curvature upon deflection 
gives a value of e/u which is too large. 


HE non-relativistic Schrédinger wave equation which governs the 
motion of electrons in a uniform magnetic field has been solved by Page.! 
The solution which he obtains shows that the mean radius of curvature of 
the electron current is less than that given by the classical theory. An 
appreciable error is thus introduced in applying the classical relation to the 
determination of €/u from deflection experiments. 
Since the energies of the electrons used in deflection experiments are so 
large, it was thought that a treatment of the problem by the relativistic wave 
mechanics would be of interest. 


THE RELATIV!sTIC WAVE EQUATION 
The relativistic wave equation? for electrons in a magnetic field is 


1 OW 4rie 4x?/ .. é 
O09 me ee eh -Y — i ee’ + — FW oO (1) 
c? Ot? he h? c? 





where A is the vector potential and yu is the rest mass of the electron. Fora 
uniform field H in the X direction 


HA H 
A=- J>3 + ky. 


Eq. (1) then becomes in cylindrical coordinates r, 4, x, 


* Charles A. Coffin Fellow. 
1 Leigh Page, Phys. Rev. 36, 444 (1930). 
2 L. Brillouin, J. de Phys. et Le Radium 8, 74 (1927). 
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1 OW 2iell dy 


1 @a ( ~) 1 oy 0*y 
pom Sieoe 5 


r or\ Or r? 060? Ox? c* oat? hc 00 





4x? / | 24 rT) y=0 (2) 
h? (. ° 4c? a 
We put 
VY = Riv) X (ale imee 2rilel nyt (3) 


owe 


where € is the total energy, and m must be an integer for W to be single- 
valued. Eq. (2) then leads to the ordinary differential equations 





a*xX 
| + aX =0, (4) 
dx? 
@R 1 dR j4n2E?  4w*y2c? Qe m 
a a ee ee 
dr? r dr | h*¢* h? he 
rer? =m? ? (5) 
— ——_ - —-a}R=0, 5 
hrc? 2 j 
| in which a@ is the constant of separation. Eq. (4) states that the electrons 
| moving parallel to the field act like free particles. If we introduce the quanity 
| ee 
W = —(E* — ye‘) (6) 
2uc? 


into Eq. (5), the differential equation for R becomes identical in form with 
that deduced by Page from the non-relativistic wave equation. We have 
therefore, as in Page’s development, 


1 el] 
w= (s+) =), s= m+ b= 01,200. (7) 
2 2muc 


The energy € is thus completely quantized.* We have also for the radius of 
curvature r of the deflected electrons* 


mew x 


@H? +1 








rT = 


(8) 


in which x=(27elI/hc)r?. Relativistic electrodynamics gives for the radius 

of curvature in terms of W as defined in Eq. (6) 
Que? 
e*H? 


9 


r 


(9) 


so that the same correction is introduced here by the wave mechanics treat- 
ment as in the non-relativistic case. 


3 If the restriction of m to integral values is removed, the energy is no longer quantized. 
However, k remains integral in order that the series for R shall terminate, and the validity of 
subsequent calculations in both the non-relativistic and the relativistic case is not affected. 

‘ Page, Eq. (33). 
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THE CURRENT 


It is necessary to deduce from Eq. (1) an expression for the current 
density. If the differential equation satisfied by y is multiplied by WY and the 
result subtracted from the product of Eq. (1) and y, we get 


- a 1 ay a°y 4rie - - 
iv-W — wrew — — (7 my ) ———A-(W + 9TH) = 0. (10) 
¢* 1€ 


ot? ar 
Also 
YV-VY — WV VY = 0-H — YY), 
and 
A: (WVY + ¥VY) = V- (AW). 
In addition we have 
_0y ay d( _¢@ y 
oo eS Lyi (tog “) ; 
Or? or? atl Ot y 


4 € 0 Ww 
— 4n1—- — ) 
” h oat 


using Eq. (3). Eq. (10) may now be written 


ehc? = = oe i) + # —_— - 
mie” (Iv yvy ie vy Pride = 0. (11) 


We identify Eq. (11) with the equation of continuity so that the current 
density is 





ehc? /_ = 4dr ie 
j= (Fev - Wwe - 


he 
This result is to be compared with the non-relativistic expression® 
eh — 4m _ 

j = ——| vy — Wy — ——Apy }. (13) 

dri he 

The factor ehc?/4ai€ appearing in Eq. (12) is, like that in Eq. (13), constant 
in the coordinates so that it does not affect the calculation of the mean radii 
of curvature averaged with respect to the current density. Hence the 
relativistic wave equation gives the same values for the mean radii of 
curvature as those calculated by Page.® 


CONCLUSION 


The relativistic wave equation leads to the same wave function for electrons 
in a uniform magnetic field as that determined from the Schrédinger equation. 


® Condon and Morse, Quantum Mechanics, p. 30. 
® Page, Eqs. (45) and (48). 
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The electron motion at right angles to the field is found to be completely 
quantized as before. The relativistic expression for the current density differs 
from the non-relativistic only in that E/c? replaces uw; this difference does not 
affect the calculation of the mean radii of curvature. We are led to a con- 
clusion similar to that reached by the non-relativistic treatment: the classical 
relativistic relation between €/u and the mean radius of curvature upon de- 
flection gives a value for €/u which is too large. 

In conclusion the writer wishes to acknowledge his indebtedness to Pro- 


fessor Page for suggesting this problem and for his helpful criticism during its 
consideration. 
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THE HYPERFINE STRUCTURE OF S AND P TERMS OF 
TWO ELECTRON ATOMS WITH SPECIAL 
REFERENCE TO Lit* 


By G. Breit AND F. W. DoERMANN 
DEPARTMENT OF Puysics, NEW YORK UNIVERSITY 


(Received October 28, 1930) 


ABSTRACT 


The proper form of the interaction energy between the nuclear magnetic moment 
and the electronic system of a many electron atom is discussed. The results of the 
Dirac equation for a single electron are taken as the guiding principle. A form for the 
interaction energy is set up as an expression involving Pauli’s spin matrices. No 
convergence difficulties occur in the form here given. 

The interaction Hamiltonian is applied to the *S and *P terms of two electron 
atoms for the case of Russell-Saunder’s coupling. An exact formula is derived for the 
resultant hyperfine structure of *S terms and corrections to the Goudsmit Bacher 
formulas for *P terms are given. It is shown that the Landé interval rule for *S 
hyperfine structure levels is exact, and that therefore the ratio of intervals can be 
used to determine nuclear spin moments. 

The formulas are applied to the Lit 5485A line. Proper functions for S levels 
are worked out by the variational method and applied to the calculation of the 
magnetic moment of Li;. With Schiiler’s wave-length data the nuclear g factor is 
2.13 on the assumption that the nuclear spin is 3/2. 

The accuracy of the calculation is discussed. It is likely to be good to at least 
2% in g. As a by-product of the calculation the lowest energy level of ortholithium 
has been computed as S —1.1354 in units of the ionization potential of Litt. The 
empirical value is —1.1358. 


INTRODUCTION 


HE effect of nuclear spin on optical spectra has been treated mathemat- 

ically by Casimir,! Fermi,? and Hargreaves.’ A somewhat more qualita- 
tive discussion has been given by Goudsmit and Bacher.‘ In the first three 
of these discussions only one valence electron has been taken into account 
while in the last the more general problem of several valence electrons has 
been considered. Qualitatively the conclusions of Goudsmit and Bacher are 
doubtless correct. They do not pretend, however, to be exact quantitatively 
since the coupling between the nuclear spin and the electron system is sup- 
posed to be taking place only through one of the electrons—the most closely 
bound s electron. It is desirable to have a more exact quantitative theory 
taking into account the coupling of all electrons. This is particularly neces- 
sary if, say, two of them are s electrons. The desirability of having such a 
theory lies in its application to the absolute value of nuclear magnetic mo- 


1 See Goudsmit and Young, Nature, March 22, 1930. 

2 E. Fermi, Zeits. f. Physik 60, 320 (1930). 

* J. Hargreaves, Proc. Roy. Soc. 124, 568 (1929); 127, 141 (1930); 127, 407 (1930). 
*S. Goudsmit and R. F. Bacher, Phys. Rev. 34, 1501 (1929). 
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ments. The observations of Schiiler' and of Granath® show that the *S level 
of Li* is split into three components under the action of the nuclear magnetic 
field and the separation between the two extreme components can be said to 
be 1.03 cm.—'. There seems to be no doubt that the three components in 
question are single. They are well and definitely resolved and so there is 
every reason to expect eventually a reliable determination of the nuclear mag- 
netic moment from the determination of their frequency difference. The 
whole electron system consists in this case of only two electrons. The solu- 
tion of the Schroedinger wave equation offers in this case a simpler problem 
than that of the heavy alkalies. An attempt at its approximate solution and 
its application to the calculation of the magnetic moment of the Li; nucleus 
is one of the objects of this paper. The other object is to formulate a theory 
for the interaction between the nucleus and the extranuclear electrons. A 
few words must be said as to why the formulation of such a theory is still 
subject to speculation. 

For one electron an unambiguous result can be obtained by means of 
Dirac’s equation. This has been done by Fermi and Casimir. The problem 
is not so clear when two or more electrons are discussed. There exists no satis- 
factory mathematical treatment of the relativistic two electron problem. It 
becomes necessary to use approximate equations. The simplest form for 
these is obtained through the introduction of Pauli’s two row, two column 
spin matrices. It is known that interactions between electrons can be satis- 
factorily represented by means of these as long as the electronic velocities 
are not excessive. It might be expected that the interaction of nuclei and elec- 
trons can be reduced to the same basis. There is one important distinction 
between these two cases. The force between a nucleus and an electron is 
attractive while that between two electrons is repulsive. When an electron 
is close to the nucleus its velocity cannot be treated as small. For s electrons 
an important contribution to the interaction energy is due to the influence 
of the nucleus when the electron is very close to it. We have no reason to ex- 
pect the interaction energy of the nucleus and the electron to be of the same 
form as that used for the approximate interaction energy of two electron spins. 
In fact it will be seen below that the two forms are different. 

In order to derive the proper form of the interaction energy we consider 
first of all the one electron problem from the point of view of Dirac’s equation, 
i.e. we write down the formulas found already in Fermi’s paper. The two 
smaller Dirac wave-function components are then eliminated and an equation 
in the two large components is obtained. This equation involves Pauli’s 
spin matrices and is equivalent to Dirac’s for the discussion of the lighter 
nuclei. Having thus derived the proper interaction energy for one electron 
the result is generalized to the case of two or more. This forms the first sec- 
tion of the paper. In the second section the separations of spectral terms caused 


5 H. Schiiler, Zeits. f. Physik 42, 487; (1927). 

*L. P. Granath, Phys. Rev. 36, 1018, (1930), (letter). See also the forthcoming paper by 
P. Giittinger in the Zeits. f. Physik in which the feasibility of i=3/2 is explained and an 
approximate g value for the nucleus is derived. 
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by nuclear interactions are worked out. In the third the results of the second 
section are applied to Li.* 


(1). THE INTERACTION ENERGY FOR ONE ELECTRON 


We consider one electron of charge—e under the influence of a single 
nucleus. The electrostatic potential due to the nucleus is Ao. The electrosta- 
tic field due to the nucleus is central so that Ao is a function only of 7, the 
distance from the nucleus. The vector potential due to the nucleus is A. 
In special applications we take 


as. a) 
With the abbreviations 
p=P+(e/c)A, P= (h/2ri)V, po = E/c + (e/c)Ao (2) 
Dirac’s equation for a stationary state is 
(Po + arp: + acpe + asps + aymc)y = 0. (3) 


We suppose the a, to be those given in Dirac’s original paper. Equation (3) 
is, of course, a set of four equations in four y,. It is convenient to eliminate 
¥1, We as has been done by Darwin.’ The reader will verify that without 
approximations 

\ 


) po — mc — (p8)— _— ‘y=0 (4) 

| Po + me f 
where now the o;,(k=1, 2, 3) are Pauli’s matrices and the column matrix 
contains only two rows: w3, ¥4. Equation (4) must be satisfied for any solution 
of (3). It is not altogether equivalent to (3) if it is used for the determination 
of the eigenwerte E. In (3) it is required that ¥y* Yitye*Wotys*Wst+Ws* Wo 
be integrable. Equation (4) is equivalent to (3) only if this last condition 
is added. If ordinary perturbation methods of quantum theory are applied 
to (4) it is more convenient to use the integrability of ¥s* Ws +yW4* Ws as the 
restricting condition. As will be seen presently this circumstance is only of 
secondary importance. For the present we use (4) and derive its conse- 
quences as though y3*¥3+yW4*Ws4 were integrable. Performing the operations 
indicated in the second term, using (2) and introducing the electric and mag- 
netic fields by means of 





€ = — grad.Ao, H = curlA (5) 
we obtain 
c ¢ e* he 
E = mc? — eAyp + — —__| +—(PA + AP) + —A?+ (aa) | 
E+ mc? + eAg ¢ c? 2nc 
h ce e€ 
+ — — EP) +—(€A) + i[E x P]¢ 6 
2mi (E + me? + =I ) o ) (6) 


+ iA [x ale| 


7C. G. Darwin, Proc. Roy. Soc. 118, 654 (1928). 
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this is still an exact consequence of (2). At this point we make an approxima- 
tion. The magnetic moment yu enters (6) only through A. It is sufficient to 
consider first order effects in uw and therefore in A. Quadratic terms in A 
are disregarded from here on. The perturbation energy may therefore be 
taken to be 


c he 
H’ = “|< (PA + AP) + — (39) | 
E + mc? + eAgLe dae 
(7) 
hc*e Je 
ii se heisinncenesttamiennsiiie [le x ale = “(EA)|. 
2r( E + mc? + eAy)? € 


In this expression the first bracketed term corresponds to the ordinary inter- 
action energy which is used in the ordinary discussions of two electrons. In 
the major portion of the configuration space eAo is negligible in comparison 
with E+mc which itself is very nearly 2mc.? In this region the first bracketed 
term of (7) is on using (1) 


” e eh 
HH” = —- = — (V6) ce) — (8) 
mec ry ag sen 


= |r x P| (8’) 


where 


is the angular momentum operator and Y operates only on r~' in the < >° 
Only the expression (8) can be derived from the analogy of the electron and 
the nucleus to little magnets. For s terms the second term of (8) becomes in- 
determinate. A definite theory can be formulated by means of (8) if w is 
defined as the limit of a spatial distribution of magnetic moment of very 
small extension in the limiting case of the extension approaching zero. It will 
be seen that the result of this limiting process is different from the applica- 
tion of the complete expression (7). For s electrons the result of (7) is (—2) 
times the result of (8). Hargreaves obtains zero for the perturbing effect of 
(8) in the case of s terms. He could have obtained any desired result since 
he has not defined the limiting process for the definition of the divergent in- 
tegrals. Nevertheless his conclusion is correct. The exact effect of the first 
bracketed term of (7) is zero for s terms, the presence of eA,» in the denomi- 
nator making the radial integrals convergent. We see that the picture of 
the nucleus and electron as little magnets can be applied in this instance only 
with extreme caution. 

If the field is central E is perpendicular to A and the imaginary part of 
the second bracketed term of (7) disappears. Its real part containing [EX A] 
is then combined advantageously with the second part of the first bracketed 
term containing (3Cé). The perturbation energy becomes 
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as (e/mc) (Mu) 1 
1+ (EF — me? + cAop)/2mec? Fr 
hec/2x (usd) 3(ré) (ru) 
t= =, 
E + mc? + eAg ( g* r? 
(hec/2x)(e| E}) (us) (ré)( ru) 
+ 1 ~ (9) 
(E+ me* + eAo)* \ r‘ 


The perturbation energy (9) can now be applied to the calculation of the 
splitting of an s term. The first term of (9) can be disregarded because 
My =0 for s terms and because the radial integration gives a finite result 
on account of eA, in the denominator. It is well known that the conser- 
vation of angular momentum makes it sufficient to calculate the average 
value of the coefficient of u.c- in the second line of (9). If this average be w 
and the nuclear spin & the s level is split into two components displaced by 
the amounts w, (—1—7')w from the normal position.’ The first component 
is the one of higher fine quantum number f=k+4. The first term of the sec- 
ond line of (9) contributes nothing to the average, again on account of the 
presence of eA» in the denominator which makes the integration over r 
give a finite result. Using (5) we have on partial integration 


* S8ru hec yy’ Sau he 
a a a 
. 3 Iw E+me* + eAy 3 \4armec 


0 


in agreement with Fermi. 

The final result is not quite correct, but is very nearly so as long as the 
region in which eA» < <mc? is negligible in comparison with the spatial exten- 
tion of YW’ and as long as E—mc? < <mc*. Both of these conditions are satis- 
fied by light nuclei. The difference between (10) and the correct result is of 
the order of the square of the fine structure constant.’ This example shows 
that the substitution of the integrability of Y3*W3+w W4*W, for the integrability 
of PitYityo*vet+ Ys*vs+ Yas is safe for the present purpose. We see also 
that the “empirical terms” used by Hargreaves are essentially the same as the 
terms derived above from Dirac’s equation. 

If instead of using Dirac’s equation we were to use the perturbation en- 
ergy (8) we should run into convergence difficulties. These may be avoided 
by endowing the magnetization of the nucleus with a finite spatial extension, 
say, by supposing that the magnetization is spread uniformly through the 
volume of a very small sphere. The first term then has no effect. The second 
gives rise to 


beh v(P) eee 1 = ‘ 
oo a f ;, fu Par. --) dV p-dV p 
4amc 3 pp 


8 See e.g. G. Breit, Phys. Rev. 35, 1447 Equations (2), (3) (1930). 
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where P’ refers to a point of the magnetization distribution, and f(P’)u 
dV, is the amount of magnetization in the element of volume dV ,,. The above 
double integral is easily evaluated if I= /f(P’)(1/rep-)dVp- is interpreted 
as the electrostatic potential at the point P due to a charge distribution of 
density f(P’). The total charge is 1=/f(P’)dVp-. Outside the nucleus 
ApI=0. Inside the nucleus Poisson’s equation gives 


Apl + 4nf(P) = 0, 


[ener = — dr. 


In the limit of a very concentrated distribution f(P’) we may set therefore 


and therefore 





Ap faery = — 4ri(P) (11) 
rpp 
and we obtain then 
w= _, a -~y?(0). (12) 
3. 4nrmc 


This result must be multiplied by (—2) in order to give Fermi’s result (10). 
The magnetic doublet model and Dirac’s equation give in this case widely 
different results. The present evidence is that Dirac’s equation is to be pre- 
ferred so that (10) is the correct result. 

The perturbation energy (9) must now be generalized to include the case 
of two or more electrons. As has been explained in the introduction we are 
unable to start with a relativistic treatment of two electrons and we must 
instead use (9) as an indication for the proper generalization. A number of 
possible generalizations could be made. Thus one could go back to (7) and 
use the field quantities E, 5¢ taking E, to be the electric intensity at x1, 41, 
z,, due to the nucleus and also due to electron 2. Or else one may take &, 
5; to be due only to the nucleus. For the first choice we would have to prove 
that the imaginary term in (7) contributes nothing to the result. Also some 
generalization of Ay would have to be made. Having no basis for forming 
such a generalization we make instead 


H’ = Hy’ + Hi,’ (13) 


where /7/,’ is the result of using x1, V1, 2:1, 61, instead of x,y,z, 6 in all of the terms 
of (9). The coupling between the nucleus and each electron is thus represented 
by a term in the Hamiltonian which depends only on the coordinates of that 
electron. 

It mayalso be shown that for *S terms the alternative generalization based 
on (7) leads to the same final result as (13) provided &;= —grad Ao; (¢=1, 2). 
The proof is not of sufficient interest to be quoted. 
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II. CALCULATION OF PERTURBATIONS IN TERMS OF SCHROEDINGER 
v IN CONFIGURATION SPACE 


We use the method of sums.’ Consider first the 2S terms. It is convenient 
to define the following spin functions 


S; == S,'S,*, So = 2-12(S,'S 3? “+ Sa?7Sg'), J.) = Ss'Sp?. (14) 


The upper indices refer to electrons 1 and 2. The functions S., Ss are as 
usual the column matrices (§), (7). The upper row corresponds to the spin 
being in the positive direction of the z axis. In the case of Russell-Saunders 
coupling there are three unperturbed functions for the *S state. Each of these 
is the product of a coordinate function and of one of the three functions (14). 
The coordinate function Y depends only on the shape and size of the elec- 
tronic triangle and not on its orientation.'!® The wave-function of the cou- 
pled system nucleus+electrons contains also the angular momentum speci- 
fication of the nucleus. Denoting the angular momentum of the nucleus by 
k there are 2k+1 nuclear functions 


fs 0) 0 
() 1 0 
, 0) : 0 . ‘ 7 
Ne=]-} Nenr=]°-b--: Nee =]: (15) 
‘ 0) 
0) 0 1 


The vector matrix uw operates on these. The matrices for the components of 
u we take to be the same as Fermi’s. Among these the only matrix having 
diagonal elements is w,. It is a diagonal matrix with elements u(1, k—1/k, 

-++—1). The complete system of 3(2k+1) unperturbed eigenfunctions 
can be represented by means of 


sey (0 Oo ) (16) 
u;™ = §;N™ 5 
meh k-1,---—8& 











= a - 
Angular momentum; k+1 k k—1 | —------- —k+1 —k —k-1 
u,* u,*-! “2 | ------- u,* 
Uo* Nok 1 —_— — Uy **1 uy * 
u_| ——--—-— | u_y*t? 4_y-*tt u_,—* 











In Slater’s method we are concerned only with diagonal elements of the per- 
turbation energy matrix. The functions N™ are normal orthogonal and so 
are the S;. We are therefore concerned only with coefficients of uw... This 
shows that the combined effect of the perturbations due to My and Mby is 
zero. In fact this perturbation is of the form (B,+B.)u. It, by itself, com- 
mutes with the total angular momentum. When it operates on 1," it gives 
(B,.+Be,)uu,;* as well as other terms which do not count in the calculation 


9 J.C. Slater, Phys. Rev. 34, 1293 (1929). 
10 E. Wigner, Zeits. f. Physik 43, 624 (1927). 
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of the diagonal element fu;**H’u,*. The operation on uy gives —(B,.+ 
Bz,)uu_.. Thus fu**H'u,* = —fu;**H'u,*=0. The term of highest fine 
structure quantum number is therefore undisturbed. Considering the other 
functions in turn we see that there is no splitting due to terms not involving 
6;, 6. This may also be verified by a short calculation of the average value 
of A(r:)0/d¢: +A (re)d/dd2 in the six dimensional configuration space x; 
- ++ g/o for the case of Y=yW(r, 72, cosO2). The terms in 6, 6: are easily 
treated if it is observed that, so far as the calculation of diagonal matrix 
elements is concerned, the symmetric operator A;é,+ A.é. gives rise only 
to the following terms 


(Ai6, + Asé2)(S1, So, S-1) > (Au: + Azz)(Si, 0, — Si). (17) 


Letting now the energy perturbations of the three levels be w;, we, Ws, in 
the order of decreasing fine quantum numbers, we have from the first three 
columns of the table 





k—1 
w= fvu. +Ax)p, Wit we = —_ ee W + We + Ws 
k—72 
= ( —_ 1 aa )m 
k 

1 k+1 
weo,il-—» -——j}. 18 
( k k - 


Here 








hec , (321° 1 2Ze? 
Ay, = . — 2mc?{ —- — —) + Zi" pe 
24(2mc? + Ze?/r;)? re 73 re) 


Performing one of the integrations over 1 for A;. and over 2 for Az, and then 
using the symmetry of ¥” in 1 and 2 we have 


w, = 1610/3 uno [ v0, 0,0; x, y,z)dx dydz 


(19) 
(uo = eh/4amc > 0). 


The only approximations made in deriving (19) are: (1) a definite case of 
Russell-Saunders coupling, (2) neglect of the square of the fine structure 
constant. The integral in (19) represents the probability of one electron 
being in a unit volume near the nucleus. The connection of this result with 
Goudsmit and Bacher is seen if one writes 


W = (1/2"/2) [bn(qudm(q2) — on(q2)m(qi) |. (20) 


If m refers to the inner electron state 


[vo q)dq © 3¢,7(0), and 


w, = (82/3)uH0n7(0). 
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In this approximation (19) (18) are identical with the corresponding formulas 
of Goudsmit and Bacher. A discussion of the difference between (19) and 
(21) will be given in the third section. 

If k=1/2 formula (18) remains correct if it is read from left to right and 
the last term is omitted. This is obvious from the order in which 7, we, ws 
have been derived. 

We may now consider *P terms. There are three linearly independent 
coordinate eigenfunctions, each of these combining with the three spin func- 
tions (14) and thus forming a set of nine functions. The coordinate functions 
may be called yi, Yo, Y_1 corresponding to components 1, 0, —1 of the orbital 
angular momentum in the z direction. The functions ¥;S; do not themselves 
correspond to definite parts of the *P term but certain linear combinations 
of them do. These linear combinations are 


3Py : 3-1/2(W1S1 — WoSo + ¥-1S-1) 
3P, : 2-12(WiS9 — WoS1), 27-1/2(WiS_1 — W181), 2-2/2(— WiiSo + WoS_1) (22) 
3P os WiSi1, 2-/2(WoSi + WiSo), 67'/2(WiS_a + 2WoSo + ¥-iSi), 

2-1/2 oS_1 + P_-1S0), Y-1S-1. 


The functions Y, Yo, Y_1 are supposed to be chosen so as to transform them- 
selves under rotations as 2-!/2(x+iy), —2z,271/*(—x+7y). The set of functions 
(22) is convenient because each one of them corresponds to a definite value 
of the angular momentum in the direction of the z axis. There is no splitting 
for the *P») and we consider first the *P; term. The nucleus is brought in by 
means of the NV, just as in the case of *S terms. The table which led to (18) 
is still correct, though now the functions «,” are different. Equation (18) still 
is valid although (19) is different. The validity of (18) is of course a conse- 
quence of the “cosine” law of interaction with the nucleus. It is very easy 
to prove this by writing #;."=7,N". Denoting the coefficient of uw. in the in- 
teraction energy by a, we have on applying the diagonal sum theorem from 
both ends of the table w,;=J,;= —J_1; we+wi=(Rk—1)/kIi+]0= —((k—1)/k) 
I_4—Io where J,;=pfv,*a.2;. This shows that J) =0 and (18) follows at once. 
It remains to determine w;. We may omit the correction in the denominator 
of the first term of (9). Thesecond two terms are easily combined and we 
have 


WW, = (eu/mc) foter 3M. a ro? M:)?1 at b for@e. + B.62)2, (23) 


where 











hec { ry c 2Ze° \ 
B, = — Poa vb 2mc? 32 2 + 21fi¢. (24) 
2m, 2mc? + (Ze?, r:)}2 ] j 


Here Z is a unit vector along the z axis. B» is obtained from B; on substitut- 
ing 22, 72 everywhere for 2, 7;. Performing the operations ¢,, 6: in (23) we 
have on letting 
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B=B,+B; (25) 


eu eu 


footer 3M i. + r2°*M22)Wo 
(26) 
+ (u/2) [ {Yo Bao — Ws*2-1(Bs + iB,)o — Yot2-!™(Bs — iB,)Ws} 


force 3My. + ro 3 Me.)yi + 


2me 2mc 


This form for w,; is correct quite regardless of any approximations which may 
be made as to the representation of the coordinate wave-functions Y, Po, Y-1. 
The most general forms of these functions are known and formula (26) is 
capable of giving exact results as long as the hypothesis of Russell-Saunders 
coupling applies."' In practical applications it appears to be unnecessary to 
use the general expressions for Yi, Yo, Y_1 and it will usually be sufficient to 
use the approximate representation of these coordinate functions as antisym- 
metric combinations of products of functions involving the coordinates of 
only one particle. We may then write 


Wi = (3'/2/4r)(F sin Oe: — F sin 62e'*:) 


vo = — (6'/2/4r)(F cos 0; — F cos 62) (27) 
Yui = — (3'/2/4r)(F sin 0,e7**: — F sin 0se~ #2) 
with 
P= frdflrd; 4 [ Prtntdrdrs = 1 o 
F = filre)fo(ri). 


The functions in (27) are chosen so as to give (22). The radial function /; is 
characteristic of the p state while f. describes the s state. If a single electron 
were in that s state it would be described by a normalized function y,. 


Clearly 
4ry.(0)? = f2(0)? / ( J srorsar.). (28) 
/ 


It is now found on substitution of (27) (27’) 


(e/2mc) fvrocan. + ro 4Me:)~, = (eh ‘4armc)(ry~8) (29) 


where 


. 


(r3) = 4 [rr 37 °ro"dridre (29’) 


is the average of r~* for a p electron described by a radial function fi. Also 
in 


3 [vorBan 


11 E, Wigner, reference 10; G. Breit, Phys. Rev. 35, 569 (1930). 
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we need 


A numerical estimate shows that the term in 7;-4 is much too small ordinarily 
to be worried about and this as well as similar terms will be neglected here. 
The most important contribution to (26) comes from terms of the type 


(3/16n2) fF cost O,Bae = (2n/3)noba(0) 
using (28). The remaining contributions are from 
ni sf (v2 2(B, + iByWo + Wo*2-"*(Bz — iByyi} = (3/5)uo(ri*). 
It is thus found on substituting all the integrals into (26) that 


wW1 = puyol (44/3)y.(0)? + 2(r-4),}. (30) 


This completes the calculation for the *P, level. 

For the *P, term we have five electronic functions given by the last row 
of (22). These may be designated as v; and the complete specification of the 
system is given by u;,"=v,N". Only now /=(2, 1, 0, —1, —2) and the 1, 
Vo, V-1 are, of course, now different from those just used for the *P; term. In 
this case the table following equation (16) must be extended so as to have 
five rows. The method of sums does not suffice now to give all of the separa- 
tion ratios. Denoting as before the coefficient of uw. by az and letting J/= 
Jvi*a.v, the method of sums shows directly that J2=—J2, = —TJ-1, 

o=0 but it does not establish a relation between J; and J2. A calculation is 
easily made however for both J; and J2. We have: 


e 
I, = fotos + ro Mo.)~i + wf viBab. (31) 


me 


Substitution of the integrals used in deriving (30) gives 


I = ppol(8x/3)¥.(0)? + (8/5)(r-*) 5]. (32) 
Also 


Ty = (ue/2mc) forccan. + ro 3M. 
(33) 


+ (u/2) J verBve + yit2-2(B, + iBy Wo + o*2-!!"(Bz — iB, vi} 


so that J; can be obtained directly from (30) by subtracting (6/5) (r3)5 
inside the parenthesis. We have 














HYPERFINE STRUCTURE 1743 
Ty = ppol(4r/3)p.(0)? + (4/5)(r-*)p] = T2/2. (35) 
Using the sum theorem we obtain for the first order energy perturbations 
3P. = pol (4r/3)¥.(0)? + (4/5)(r-), | 


k—2 3 k+ 3 (k + 1) 
(2,828, 3, 48, 4m) ag 
h k k k 


which may be compared with 





3P, 


a 1 ok+t 
puol(4n/3W(0)8 + 20Fo](1,-—» -———) G6” 


3S l 


—_ 1 k+1 
upuo(8a/3)(1 + ova(oy*(1, = ~—). (36”) 


Here ¢€ is a fractional correction due to the difference between (19) and (21). 


III. APPLICATION TO THE SPECTRUM OF LI* 


The pattern of the *S—*P line of Lit which has been observed by Schiiler 
and Granath will now be considered using formulas (36) (36’) (36’’). The 
quantities (7%), and ¢ are of the nature of correction terms and neglecting 
them these formulas become identical with those of Goudsmit and Young. The 
y.(0)? for a hydrogenic s electron is related to the fine structure constant by 


2ry5(0)*u0? = Ra®Z?/n3 (37) 


which is convenient for calculation. For Z=3 and n=1 we have in agree- 
ment with Goudsmit and Young (8/3) Ra?Z*/1840n' =0.228 cm=. Letting 
(1840u/u0) =g(k) we have therefore the approximate formulas: 


ee 9(2 k—2 3 k+3 nes, 

= ae 7 yy Z, - 9 — . — P — 

2 (g(k, bk ; A A 
1 k+ 1 

IP, = 0.228(¢(#)/4) (1, on eel an tape -) (37) 
k kk, 

™ 1 k+1 
S, = 0.228(¢(8)/2) (1, a eda ) 


which are essentially the same as those of Goudsmit and Bacher. An applica- 
tion of (37) and the calculation of relative intensities of the different compo- 
nents of the pattern can be made for different values of k. As has been 
pointed out by Granath the experimental pattern seems to agree best with 
k=3/2. Most of the pattern appears to have a complicated structure which 
is due to unresolved components. Only the group of lines due to transitions 
from *Po to *S; is simple and consists of three fairly well resolved lines. It 
seems difficult at present to draw any definite conclusions from the other 
parts of the pattern on account of lack of resolution. We consider it more im- 
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portant therefore to have an accurate calculation of the correction factor 
1+€ in (36’’) than the exact evaluation of (36) (36’). There is, of course, no 
difficulty in evaluating (y-*),. In the present instance the correction due to 
(7-8), amounts to about 1 percent change of scale for *P: and 2 percent for 
3P. 

In order to calculate the factor 1+ € in formula (36’’) we must use the exact 
formula (19) which differs from (37) by the correction factor 


ite =2 { v0, g)dq/is(0)?. (38) 


The integral in the numerator is here the same as the integral in (19) the 
letter g having been written collectively for (v, y, s). The denominator of 
(38) is the ¥,(0)* of (37’) for n=1, Z=3. 

We have tried a number of approximate forms for ¥(g:, gz) using the 
variational method for the calculation of the term value of the *S state and 
then, having adjusted the constants to satisfy the minimizing requirement, 
calculated by means of the minimized functions the corresponding value of 
1+e. In all cases € is a small quantity. Some of the functions give eigen- 
werte in good agreement with the spectroscopic term value while others give 
considerably too small absolute values of the eigenwert. It is of interest to 
see the results of the bad as well as the good functions. A very rough ap- 
proximation to the solution can be obtained by tentatively expressing y as 
an antisymmetric combination of products of hydrogenic functions in the 
field of a certain nuclear charge. This charge is deliberately taken to be the 
same for both functions. The two functions entering are orthogonal and the 
calculation is very easy, only the common scale of the two functions being 
varied. In terms of the ionizing potential of Lit* the energy of the 2°S state 
is known to be A= —1.1358. In the same units the above function gives 
A= —1.1279 and 1+¢€=0.96. This is a deliberately poor function, possessing 
only one adjustable constant and giving a relatively poor agreement with the 
experimental \. Hylleraas’ has considered the *S terms of two electron 
spectra particularly for He. In his paper functions are tried for the ortho 
state of He which may be called the one and the two term approximations. 
These give for He X= — 1.0855 and X= —1.0871. The values of 1+ € for He 
would be 1.009 and 1.04 respectively. For Lit we have used only the one 
term function of Hylleraas giving \= —1.1337, 1+¢€=1.02(6). From the 
analogy to He one would suppose that the two term function of Hylleraas 
would give an appreciably larger «. The calculations with the two term 
function of Hylleraas are complicated, however, and we have found it more 
convenient to use another three constant function. It has also been possible 
to correct the first function of Hylleraas by a numerical method. In the 
latter case there is no definite way of using the eigenwert agreement as a 
check on the calculations and the procedure itself is somewhat arbitrary. 
Nevertheless it appears to be definite enough to show that 1+¢€>1.02 and 


12 Hylleraas, Zeits. f. Physik 54, 347 (1929). 
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that 1+¢ 1.06. The numerical method is based on the following considera- 
tion. In the units of length (h?/87°mZe*) and energy (RhcZ*) used by Hyl- 
leraas the differential equation to be satisfied by y is 


0" 2 @ 0° z @ 1 1 1 a) , 0 
ss ce cao oe ame me res Fons ok 
Or,’ r, Or, Ors ro OFe ry = re"/ sin 6 06 06 


(39) 


The angle @ is the angle between the radii 7, 72 and ry». is the distance between 
electrons 1 and 2. A certain approximation to the exact solution can be ob- 
tained by supposing that y is a function only of 7; and rz the angle 6 being not 
very important for this particular state. The one term function of Hylleraas 
does not involve @ or ry» and is therefore an approximation of the above type. 
A systematic neglect of the dependence on @ leads of course to a definite 
eigenwert and to a definite eigenfunction if the variational method of solu- 
tion is used. Substituting a general function ¥(m, re) into the variational 
equation it is easy to obtain a differential equation for this Y. This equation 


0? 2 2a 0? 2 fa) Xr 1 1 1 
eae ig Ae come: lp seems ie eee samen ae ears ay ane lp came aoe oe BE Oe 
or," al or, Ore” reo Ore 4 r) re 
p = 7; if Lal > ro 
- (39’) 

p=fPrell ry) < Pe. 
We assume for the present that the eigenvalue of \ and the function yw of 
(39’) is a sufficient approximation to \ and W of (39). For the 4S state y 
is antisymmetric in 7“, 7. and vanishes therefore at r,;=72. It may be repre- 
sented in two dimensions. Along the axis of re =0, =r. The last three terms 
behave as 1/ro. If W and its first and second derivatives are finite on this 
axis we must require in order to satisfy (39’) that the term 2/r20W/0r2+ W/te 


should be finite i.e. 
0 log y 1 
rev) 2d - 
Ore ro=0 2 


Any trial Y can be checked qualitatively by means of (40) and it is important 
to require that (40) be satisfied in the region of large numerical values of y. 
This is especially advisable in the present instance because the result we wish 
to obtain is proportional to fw(n, 0)*7°dr, and depends primarily on the 
relative values of y? at re=0 as compared with its value at other points of 
the , 72 plane. 

The one term function of Hylleraas was tested by means of (40) and it 
was found that although the logarithmic derivative was fairly close to its 
correct value the differences were greater than could be desired. At a certain 
point of the 7; axis (40) is satisfied exactly. This point lies at a higher value 
of 7; =719 than that which corresponds to the maximum of r,°y*. At the maxi- 
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mum yw does not decrease rapidly enough in absolute value with increasing 
ro. The dependence of wy on rz has therefore been corrected by a numerical 
procedure so as to satisfy (40) in most of the important region. For values 
of 71;<110 this was done by fitting the exponential curve Ce’? to the curve 
of ¥ against 7. The fit was made quite accurate by calculating logarithmic 
derivatives 0 log ¥/Orz of the Hylleraas function and choosing such an f2 


for every 7, that 0 log ¥/dr2=—}4. This procedure is of course somewhat 
arbitrary because e~’*? is not an exact solution of the wave equation and 
because 0 log ¥/dr.= —} only at r2=0. Nevertheless, as long as the cor- 


rection applies only to a small region in the neighborhood of r2=0 it is es- 
sentially correct. For values of 7; >/7yo it is impossible to correct the function 








rm ir ni h“/Anin > > 


Fig. 1. The dotted curve is computed from the corrected, the full drawn curve from the 
uncorrected Hylleraas function. The values of d log y/dre are indicated along the curve. 


by means of an exponential curve. The correction was made, therefore, by 
choosing for a given ™=7,'’>?io a corresponding value of 7=17’<rjo for 
which 0 log ¥/0r2 deviates from —1/2 by a numerically equal amount. The 
known corrections to for 7;=7,' and different values of r2 were now applied 
tor,=7,’’ for the same values of 72. This means that the deviation of 0 log 
y /dr. was taken as the criterion of the magnitude of the correction to y. The 
whole procedure has of course a meaning only as long as the corrections are 
small and can be applied therefore only in the region around the maximum 
of nW(n, 0). This maximum is at 7;=12.5 in the units of length used by 
Hylleraas. The correction becomes zero at 7;=18.5. The values of (n(n, 
0))? at these points are 397.5 and 270 on an arbitrary scale. The relation 
of the corrected and the uncorrected (7:~(1, 0))? curves is shown in the Fig. 1. 
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Using these corrections the changes in /r;°¥7(n, 0) dr and ffr2rc%y*(n, 
re) dridr2 were calculated. The corrections were taken into account only 
from r,=0 to 7, =24. The resultant change in 


i] nyr(n, 0)dr; /f freee, r2)dridre 


came out to be 4 percent. The value of 1+. for the one term Hylleraas func- 
tion is 1.02; the corrected value of 1+ € is 1.06. 

Our correction procedure may of course be criticized on account of its 
arbitrariness. It must be remarked however that the sign and the order of 
magnitude of the result depend primarily on the fact that an increase in the 
absolute value of is required at the maximum by the logarithmic derivative 
condition, Thus even though mathematically the procedure lacks rigor we 
believe that for the present problem it is essentially correct. 

Finally we have calculated the result by the variational method usine 
the following trial function. 


y = o(kry, kre); o _ (r; =— c)e~ a!2)r,—(b/2) 1, —_ (re _ Cem (e!2)r2— (0/2), | 


This function is easily adapted to the variational procedure. It contains 
three constants a, b, c. The constants c and k can be treated algebraically 
in the variational procedure. Extensive arithmetic need be done only in 
connection with different values of a. The reason for this simplification is 
that the minimized value of k is always close to k=1 if b=1, the difference 
from 1 being usually of the order of 1 percent or less. The first step is to 
choose a value of a and to take }=1. Supposing that k=1 the minimizing 
is done for c analytically. Taking this value of c the result is minimized for 
k analytically. It is then usually unnecessary to minimize for c again because 
k=1. The procedure is repeated for different values of a and the absolute 
minimum of} is fixed by interpolation.’ We reproduce the final formulas for 
the minimizing procedure because they may be found useful in other con- 
nections. For k=1 


/4 = (M — Ly) + L,)/N (41) 
where 
M = 8a-*h-* + 246-1a-5 — 12a—1b-18-7(38-! — 2a7) 
— c|4a-%b-? + 12a-4b-! — 8abB-*(38-! — a-)] (42) 
+ ¢[2a-'b-$ + 2a-%5-! — 4ab8-*]; B = (a + b)/2 
Ly = 24a~*b-* + 48a~5b-? — 488-7 — c[16a-%b-3 + 24a-4b-? — 408-8] 
+ ¢?[4a-*b-3 + 4a-*h-? — 88-5] (2) 


18 The trick of minimizing for k is the same as that used by Hylleraas, ref. 12, and amounts 
physically to satisfying the requirement of having the kinetic energy = — (1/2) times the poten- 
tial energy in a system with Coulomb forces. 
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ZL, = 24[a~4b-*(a + b)-* + a-8b-2(a + B)-2 + a *-"(a + )-4 
+ 2a—'b-"(a + b)->] — (33/2)8-7 — c}8[2a-3b-2(a + 8)! 
+ a-%b-2(a + b)-2 + a 2b-"(a + b)-3 + 3a "'-“"(a + 0)-4] 
— (25/2)8-®} + c2{4[a-%-2(a + b)-! + aba + b)-3] — (5/2)8-7}. 


(42°) 


Z=nuclear charge. 
N/4 = 24a-5)-3 — 188-8 — 12c(a-4b-* — B-*) + 2c2(a-3b-3 — B-S) (42’”’) 
With L =L,—JL, the minimum for k comes at 
k = L/2M, d/4 = L2/(MN). 
The correction factor 1+ € is obtained as 


1+e= k*J/(N/4) (43) 
where 


I = 24a-* — 12c(a~* — B-*) + 2c%(a-* + O-* — 28-*). (43’) 


Using these formulas we obtain for a=1/3; the minimized values c = 1.446, 
k=1.0057, —A/4=1.1344; 1+¢€=1.082. The minimized value of a is 0.380 
giving —A/4=1.1354; 1+¢€=1.063. On the other side of the best value of a 
say at a=0.395, 1+€=1.055 and —A/4=1.1353. The calculations have 
been carried out to several more significant places and for several more values 
of a than given here. It is seen that 1+ € varies slowly for the functions 
minimized for changes in c and k. We believe therefore that 1.063 is a fairly 
accurate result. It will be noted that this result is in agreement with that 
obtained by correcting Hylleraas’ one term function, and that the agreement 
between the empirical value of —A/4=1.1358 and the theoretical —\/4 
= 1.1354 indicates that the trial function used by us is good. 

It is of course usually said that an eigenfunction may be quite inaccurate 
even if a good eignwert is obtained. We must consider this in somewhat more 
detail. Eckart'™ has already given a criterion for the accuracy of eigenfunc- 
tions. His reasoning may be extended so as to apply to the specific calculation 
of a given quantity. 

With Eckart we let the solution obtained by the variational method with 
a given trial function be ¢ and the true eigenfunctions (in this case of the 
3S system) we denote in order of decreasing term values by Yi, yo, ---; 
the corresponding negatives of the energy we write as Wi, W2,--- The 
functions Yi, Ye, -** are supposed to be normalized and ¢ =ayWi+ae2x2+ - 
where a;>+a.?+ --- =1. Eckart shows that if the trial function @ leads to 


4 Carl Eckart, Phys. Rev. 36, 878 (1930). In this very useful and interesting paper will be 
found also many other eigenfunctions. The one used for *S states of two electron systems is 
very similar to ours. The eigenwert for Lit obtained by Eckart is almost as good as the one 
obtained by us. For our work we considered it important to minimize for c as well as the other 
constants because the result is fairly sensitive to c. See also D. S. Hughes and C. Eckart, Phys. 
Rev. 34, 694 (1930) for the calculation of the isotope effect in the Li* spectrum. 
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an energy value of absolute amount .W then 1—a,;2’<(W,— W)/(Wi-— W2). 
We wish to calculate the average of some quantity {(q) over the configuration 
spaceg. The function ¢ gives the approximate result fo=/f(q)¢*(q)dq while 
the correct result is f= /f(q)W2(q)dg. The error is 


fo- f = (a7 — 1)fir + a2*fee + as*fa3 + +++ + 2aidefig +--> 


where fix = /Wiwifdg. The terms in f;; are of the order of 1—a,2 which in the 
present instance is of the order of 4X10-*. The relative error due to these 
terms can be therefore neglected. The only important terms are those con- 
taining fi;(74#1). Thus 


2 


fo 7 f 


IIS 
IIS 


ai fii. (44) 


. 


: 


[Ms 


= a) 
2a, >) aif: 
t=2 


Any a? for i>1 is <1—a?<(W,—W)/W,—W.) and the general order of 
magnitude of (44) is that of ( (Wi— W)/(Wi-— We) )."*. The numbers fi; are 
however also of influence. In order to estimate their order of magnitude we 
must use the special form of f. In our problem 


f = 3(6(0, gz) + 6(0, qi) 


g: and ge being the collective coordinates of electrons 1 and 2. We have then 


fre = Jvc, q)¥2(0, g)dq. (45) 
Approximately 
Yi S 2-1/2[G16(qi)bas(q2) — bis(G2)b2s(q1) | 
2 S2-/2[G16(qi)bae(qe) — b1s(G2)b30(q1) | 


IIS 


where the ¢;, is a hydrogenic function in a central field of nuclear charge Z 
and ¢25, $3, are hydrogenic functions in a central field of nuclear charge Z —1. 
The hydrogenic functions ¢2,, ¢3, are orthogonal to each other. In this approx- 
imation fi2=(1/2)d2,(0)¢3,(0). It will now be remembered that ¢,,(0)* in a 
field of charge Z is proportional to Z*/n*. Thus ¢25(0)?/¢1,(0)? 2(2/3)8(1/2) 
=1/27 in the present instance. We may take fi2<(1/30)f. It is now seen 
that we are likely to over-estimate the error by setting ad3;=a4--- =0, and 
attributing all the error to a:. Doing this, fo—f=2a2fi2=2(a2/15)f and as 
itself is of the order of 1/15. Thus the accuracy is likely to be about (1/2) 
percent. 

This estimate of the probable error together with the agreement of this 
calculation with the numerically corrected Hylleraas one term function makes 
us think that the value 1.06 which we calculated by both methods is probably 
correct to 1 percent. We thus conclude 


1+e=1.06. (46) 


Using this value in (36’’) (37) and the separation between components (1) 
and (3) observed by Schiiler we can derive the value of u. According to Gran. 
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ath’ a slight correction must be applied to Schiiler’s value of vyae for compon- 
ent (3). Applying this correction and supposing that the nulear spin is k = 3/2 
we get 

1840u/k = g = 2.13. (46’) 


Thus Schiiler’s and Granath’s present data and the present calculation 
indicate that the magnetic moment of the Li7 nucleus is approximately equal 
to three times the theoretical magnetic moment of a proton. If the magnetic 
moment were exactly three times the theoretical protonic moment the g 
value in (46’) would be g=2. We have at present no explanation to offer 
for the result g=2.13. It appears that the difference of 2 and 2.13 is a greater 
one than can be accounted for by experimental errors or by errors in our cal- 
culation. It must be remembered however that the measured position of 
component (1) may be slightly affected by the proximity of component (2) 
and that the whole separation between (1) and (3) is of the order of 1 cm™. 

Note added in proof: If the coupling of the electronic orbits and spins is 
of the Russell-Saunders type and if the electronic interactions are not too 
large and yet considerably larger than the interactions with the nucleus it 
is possible to give general and simple formulas applying to the case of one 
electron being in an s state. The other electron is then supposed to be in an 
L state of azimuthal number /. We may write the energy perturbation as: 


w = (u/k)(b/2)[f(f + 1) —k(k +1) — JG +1)] (47) 


where the perturbation Hamiltonian is taken to be H’ =(Bu) = (u/k)(Bk) 
and k denotes the angular momentum matrix vector of the nucleus in units 
h/2r. If the electronic angular momentum J is made diagonal in the sense 
of containing 6(j’, 7’) as a factor for every matrix element and if for a given 
j the matrix J, is also diagonal then the diagonal elements of B. for that j 
mb; m=(j,j—1,:-:,-—j). The proof of this for B, is exactly similar to 
the proof for ¢, in the usual derivation of Lande’s g factor. The momentum 
J, remains conserved under the perturbation B,. We now imagine the 
coupling between the electron spins and between the two orbital momenta 
to be removed. The matrix elements may be referred to pairs of states such 
as (7S1/2)m, and (#Z741/2)m, A canonical transformation does not change the 
sum of diagonal matrix elements belonging to a given m=m,+m.. The 
values of b for *Z and 'Z are thus related to the values for °S, 7 by sum rule 
equations: 


(1 + 1)b(*Zi41) = (2)0CS) + U + 2)bCPLi41/2) 
L{b(%Liy) + B(%L1) + bCL)} 
= (3)b(?S) + 21b?Lisi2) + (@ — 3)bCLi-1:2) ete. 
The values of } for the one electron terms are known: 
b(2Lj) = Yo L(E + 1)/j(j + 1-9); BCS) = (1627/3)uo ¥s(0). 


The value of b('L,) is easily computed directly and the remaining Db are 
determined by the above relations. Thus: 
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b(®Lia1)/mo = (8x/3)Ws(0)/(l + 1) + 41(7-3)/(21 + 3) 
b(3Li)/mo = (82/3)Ws(0)/(2 + 1) + 2(r>) 

b(3Li1)/mo = — (8x/3)Ws(0)/l + 400+ 1)()/(22 — 1) (47’) 
b(*L1)/po = 2(r). 


Formulas (36), (36’) are special cases of (47), (47’). The fact that sum rela- 
tions such as used here exist has already been mentioned on p. 210 of Pauling 
and Goudsmit, Structure of Line Spectra. According to an informal com- 
munication of Professor Goudsmit these sum rules have been used by him 
in the derivation of the general formula in Phys. Rev. 35, 440 (1930). 
Formulas (47), (47’) apply only to the case of Russell-Saunders coupling 
but have the advantage of simplicity. 
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ABSTRACT 
Part I. Tue THERMOANALYsIS OF Br SINGLE CRysTALs 

Production of crystals..-Crystals of Bi of any desired orientation were grown 
by the method of Goetz, one half of each normally, the other half within a transversal 
magnetic field. The orientations, predetermined by a seed crystal, were not affected 
by this process. Two methods of growth-—-the continuous and the discontinuous— 
were used. 

The thermoanalysis of a crystal.—A method and its experimental realization— 
the thermoanalyzer—were developed to measure and to localize any changes of the 
thermoelectric properties along the lengths of crystals without applying mechanical 
stresses to them, by progressive local heating of the specimens. Thus very small dis- 
tortions and imperfections were detectable due to their thermoelectric asymmetry. 

Types of thermoanalytic diagrams.—The types of diagrams to be expected in 
the cases of a perfect single crystal, a double-, and a triple-crystal are discussed and 
a method is developed to analyze simple diagrams by means of the theory of heat- 
conduction. 


Part I]. EXPERIMENTAL RESULTS OBTAINED FOR NORMAL AND MAGNETIC 
CRYSTALS WITH A DISCUSSION OF THE SAME 


The thermoanalysis of normal crystals.—The application of the thermoanalysis 
to normal single crystals leads to the discovery of very small (0.1 mm?) regions 
within a crystal which are distorted due to very slight variations of the cooling con- 
ditions during the growth of the crystal. It is however possible to avoid these faults 
and to produce comparatively perfect specimens. Hence it was possible to test the 
results of different methods of crystal production and thus to refine the method used 
as well as to measure quantitatively the influence of a magnetic field applied to the 
crystal during the time of its formation. 

The thermoanalysis of ‘‘magnetic’”’ crystals.—It was found that the normal half 
of a crystal has a thermoelectric e.m.f. against the “magnetic” half. The sign and size 
of this e.m.f. depend on many circumstances, though mainly on the orientation the 
growing crystal has with regard to the direction of the field lines. 

The effect as a function of the orientation.—The effect is a maximum if the 
principal axis of the crystal grows normal to the lines of force and it is very small 
(probably zero) if the axis grows parallel. It depends furthermore on the orientation 
of the crystal with regard to the direction of the thermoelectric current, since the two 
orientations in which the principal axis is normal to the field show different effects: 
a small one if the axis is normal to the current, and a very large one if it is parallel. 

The thermal e.m.f. obtained for the latter case is 4.3 microvolts ‘degree which 
would correspond to a change in orientation of ca. 21° though no actual change 
could be observed. 

The effect as a function of the method of growth.—The thermal e.m.f. depends 
largely upon the method of growth, i.e., whether the crystal is grown by the con- 
tinuous or the discontinuous method. 
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The effect as a function of impurities.—The thermal e.m.f. for one and the same 
method of growth depends also on the amount of chemical impurities. Four kinds of 
purest Bi of different provenience were used in which the total amount of other sub- 
stances was less than 0.2%. The metals were spectroscopically examined and it 
appeared that the purest (electrolytic) metal showed the smallest effect. It was ca. 
40 times smaller than the effect of another metal. Traces of Ag and Pb seem to 
affect the “magnetic” sensibility of Bi most whereas Sb is rather ineffectual. 

The effect as function of the field strength.—The effect is influenced by the 
strength of the applied field through this problem is not yet settled quantitatively. 
The influence of a small field (10? Gauss) is comparatively large, whereas fields of more 
than 13000 Gauss seem to decrease the effect. The maximum seems to depend largely 
on the kind of the impurity. 

The effect as a function of temperature.—The Peltier-effect (extrapolated from 
the E—f(@) curve) shows a sharp discontinuity between 75° and 90°C, indicating a 
different relation with regard to the transformation of Bi in this region. No indica- 
tion could be found that annealing above this point (16 hours) and aging (1 month) 
destroys the difference between the normal and the “magnetic” half of a crystal. 

Discussion.—The results obtained are brought into relation with the diamagnetic 
anisotropy of Bi. They are discussed with regard to the investigations of other authors 
and it appears possible to describe the effect as due to a change in the secondary 
(mosaic) lattice of the Bi crystal. 


Part I. THE THERMOANALYsIS OF Br SINGLE CRYSTALS 


INTRODUCTION 


N A previous paper,' one of the authors (G) described a method of pro- 
ducing Bi-single crystals which permits the application of strong, trans- 
verse, magnetic fields at the zone of crystallization. The influence of the 
field with regard to the orientation of the growing crystal was studied. It 
was stated that the orientation frequently obtained was such that the tri- 
gonal axis of the growing crystal was parallel to the lines of force, in case the 
crystal formed its first center of crystallization within the field. This result 
was partly to be expected if one considers the anisotropic nature of the dia- 
magnetic susceptibility in Bi-crystals. Since this susceptibility is at a mini- 
mum along the trigonal axis and at a maximum in a direction normal to it, 
the crystal assumes an orientation which corresponds to a minimum of 
free energy. This result is in agreement with the early observations of Pliicker? 
and Leduc.* It was, however, unexpected that the orienting forces of the 
strongest magnetic fields (ca. 22000 Gauss) did not show any effect as soon 
as the orientation of the growing crystal was already predetermined before 
entering the field. Yet for several reasons, it seemed to be interesting to 
investigate whether or not a crystal, one half of which was grown without, 
the other half within a strong transverse magnetic field showed any differ- 
ences at all between portions, despite the fact that the orientation remained 

unchanged. 
Such differences were to be expected when, for instance, one considers 

1 A, Goetz, Phys. Rev. 35, 193 (1930). 


2 S. Pliicker, Pogg. Ann. 76, 583 (1849). 
3M. A. Leduc, C. R. 140, 1022 (1905). 
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such observations as those by Tieri* concerning the Hall-coefficient of Bi 
specimens crystallized within a strong magnetic field. A large difference of 
the Hall-e.m.f. was found with regard to different directions of the field 
lines crossing the solidifying specimen, but no interpretation was possible as 
these experiments were only of a qualitative kind and made on polycrystals. 

Production of crystals. Since the method of producing bismuth single 
crystals introducing a minimum of mechanical stresses has already been 
described! as well as the way the magnetic field was applied to the growing 
crystal only a few additional details need be mentioned. 

The orientation of the crystal was as usual predetermined by a seed. 
Three main orientations were grown, which represent the possible primary 
relations between trigonal axis, the (111) plane, and the lines of force (see 
below). Two methods of growth were employed which differ principally 
both in the forces applied to the growing crystal at the moment when the 
field was energized and in the results obtained. 

The first, entitled the discontinuous method, consisted of two distinct 
processes. A crystal was first grown completely without any magnetic field 
present (the residual field of the magnet was eliminated by removing the 
pole pieces). It was then removed from the growing trough, etched, and 
carefully examined with regard to irregularities. In case none could be de- 
tected, the crystal was put back into the trough, which in turn was put 
back into the furnace so far that only one half of the crystal melted again, 
the temperature of the furnace being regulated so as to bring the border 
between solid and liquid-crystal exactly in the middle of the pole-pieces of 
the magnet. This process had to be done very carefully in order to avoid 
irregularities along the molten part of the crystal. As soon as thermal equi- 
librium was reached, which condition could easily be recognized since the 
progress of the molten region into the space between the pole-pieces stopped 
as soon as the heat distribution along the crystal became stable, then the 
magnetic field was excited and the driving mechanism of the crystal appara- 
tus started, thus the second half of the crystal recrystallized within a mag- 
netic field. 

The second, entitled the continuous method, allowed the growth of the 
unmagnetic and magnetic portion of the crystal in one operation. This 
was done by growing the crystals in the usual manner but with the pole- 
pieces in position, a precaution that was essential as pole-piece moving with 
its attendent jarring was prohibitive during crystal growth. This method had 
the advantage over the other process in that the growing forces remained 
undisturbed when the magnetic field was applied. In the early work, how- 
ever, it had the disadvantage in that the residual field of the magnet (200- 
300 gauss) was present while the first portion of the crystal, the so-called 
“unmagnetic” half, was growing. Later this was eliminated by providing 
a magnetic by-pass of soft iron across the pole-pieces and above the trough, 
which was removed just before the magnet was energized. All these crystals 


4 L. Tieri and E. Persico, Linc. Rend. 30, 464 (1921). 
5 A. Goetz and M. F. Hasler, Proc. Nat. Acad. 15, 646 (1929). 
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were grown in an atmosphere of CO, to prevent oxidation so that in this re- 
spect they also differ from those grown by the discontinuous method. 

The first observation which was made on crystals produced by the dis- 
continuous method with fairly high fields was the difference in reflectivity 
of light between the two halves. This was most noticeable in the P, orienta- 
tion, i.e. where (111) plane is parallel to the lines of force. It seems improbable 
that this difference was due entirely to the second recrystallisation, as etch- 
ing off a layer one mm or more thick does not destroy the phenomenon. 
This effect must indicate a change of the structure of the crystal-faces within 
the magnitude of the wave-length of visible light and it thus seems possible 
that it concerns the mosaic lattice as described by Zwicky’ and Goetz’. This 
possibility will be treated at length in another paper. 

The thermoanalysis of a crystal. Concerning a sensitive and effective 
method for determining possible differences between these two halves of 
one and the same crystal, it seemed that the measurement of the thermo- 
electric effect fulfilled the requirements the most easily. Nevertheless, it 
was experimentally rather difficult to perform, since it was not only necessary 
to measure the thermal e.m.f., but also, to localize it as exactly as possible. 
Several methods had to be tried until it was possible to obtain reliable re- 
sults and since such methods, allowing one to measure as well as to localize 
a thermal e.m.f. within a metal rod, are excellent indicators of inhomoge- 
neities undetectable by other means, their applicability is not limited to 
this special problem. Thus it seems worth while to give a more detailed 
description of them. (A similar idea has been put into realization inde- 
pendently by Terada and his collaborators.’:!°."!_ Their experiments had been 
started for testing the reality of the Benedix effect and the method used was 
fundamentally the same as ours, though quite different in its experimental 
procedure. It is excellent for the detection of the “residual” thermal e.m.f.’s. 
but its sensitivity and reproducibility does not meet the requirements of 
the thermoanalysis of single crystals). 

A simple method first tried was as follows: The crystal was first protected 
by a coat of Duco paint. Leads from the ends were connected to a galvano- 
meter, the contacts with the crystal being kept at constant temperature. 
The crystal held horizontally, had a small portion heated locally by the touch 
of a mercury meniscus which topped a column of that metal heated by an 
electric furnace. If the meniscus was moved slowly along the crystal and the 
deflections of the galvanometer were observed with relation to the position 
of the heated point on the crystal, then a thermoelectric analysis of each 
increment of length could be made. As soon as the heated region came near 
the border between the “magnetic” and the normal part of the crystal, an 


6 A. B. Focke, Phys. Rev. 36, 319 (1930). 

7F. Zwicky, Proc. Nat. Acad. 15, 253, 816 (1929); 16, 211 (1930). 

8 A. Goetz, Proc. Nat. Acad. 16, 99 (1930). 

® T. Terada and T. Tsutsui, Proc. Imp. Acad. Tokio. 5, 132 (1929). 

10 T, Terada, T. Tsutsuiand M. Tamano, Sc. Pap. Inst. Res. Tokio, 7, 201 (1927). 
tT. Terada, T. Tsutsui and M. Tamano, Proc. Imp. Acad. Tokoi, 3, 507 (1927). 
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e.m.f. was indicated, but its size could not be measured accurately enough 
inasmuch as it depended on the contact conditions of the meniscus which 
could not be made reproducible. This method was therefore soon abandoned. 

The next method that was tried was to heat the crystal locally by focus- 
ing radiation of a powerful incandescant lamp upon it. This method had the 
advantage over the previous one in that it was not necessary to protect the 
surface of the crystal against amalgamation, with a covering that was neces- 
sarily a poor heat conductor. A large number of observations were made 
with this apparatus, but due to certain inherent difficulties of the method, 
it was abandoned. These difficulties were, for instance, providing air cooling 


































































































Fig. 1. Diagram of the thermoanalyzer. 


of the crystal as a whole while preventing local convection currents at the 
point of heating, measuring the temperature at the heated point, and pre- 
venting surface irregularities of the crystal from nullifying results. 

The final arrangement which avoided the above difficulties and made 
possible quantitative measurements is sketched in principle in Fig. 1. A 
cylindrical glass-container A, 18 cm long and 9 cm wide is held in the 
center of a large glass-container B, 24 cm wide, 22 cm long, by means of 
the bridge C. The glass tube D, is fixed in A and brings water at a definite 
temperature and under constant pressure into A. The vessel B is also filled 
partially with water, the level of which determines the position of a ring- 
shaped float E. Two tubes F, and F, fixed in B permit the raising or lowering 
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of the water level of B, which in turn changes the position of the float E. 
Fixed to this float E is a frame arrangement which carried several devices 
among which are the syphons //; and //2 and the scale L. The latter serves 
to measure the position of the float (incl. oil-layer, furnace, etc). with regard 
to the crystal by means of the reading microscope K attached to the frame 
C. The syphons J/7, and JJ, have one end in A just underneath the water level, 
the other end of /7, going outside of B, the corresponding end of JH, ending 
beneath the water level in B. The horizontal part of both syphons is shaped 
to separate air bubbles from the flowing water, while adjustments for bring- 
ing the tubes into correct positions with regard to the water levels are pro- 
vided. Beside the syphons, the float carries an electric heater consisting of 
a horseshoe-shaped thin mica sheet around which chromel wire is coiled. 
Other mica sheets above and below insulate this heater from two thin copper 
plates which bind the whole, making for equal heat-distribution and rigidity. 
This whole heater is only 1 mm thick. Its position with regard to the frame 
is so adjusted that it is slightly above the water level in A when the ends of 
IT, and I7Tz are just beneath that level. On top of the water in A, an oil layer 
of 5-10 mm thickness is put in order to surround the heater completely. 

The arrangement works as follows: A constant flow of water is sent into 
A filling it up to a certain point, any surplus being syphoned off. The water 
level in A is therefore determined by the position of the syphons JJ, and J7. 
or what amounts to the same thing, by the position of the float E in B. Thus 
a permanent circulation of the water underneath the heated oil layer results 
and consequently a constant temperature. If a certain amount of water is 
added through F;, the water in B rises and with it the float E holding the 
syphons and the heater. The oil layer in turn maintains its relative position 
with regard to the heater since its supporting water column rises as the 
syphons rise. It follows that the water has to be added into B at a rate suffi- 
ciently slow to permit hydrostatic equilibrium in A. Therefore, the cross 
sections of JZ, and J/, are large compared with the cross-sections of Fi, F2 
and D,. As an indicator of this equilibrium, a differential thermocouple is 
also attached to the frame of EF. It consists of two copper wires L; and Ly» 
the ends of which, ending one in the middle of the oil-layer, the other within 
the water flow, are connected by a very thin constantan wire. This arrange- 
ment is extremely sensitive because the smallest displacement of the heater 
with respect to the oil and water level changes the temperature of the L,- 
junction and thus produces a change in the thermal e.m.f. 

The reason the syphon JJ, is used in addition to J/, is first, to tie the water 
systems of B and A together so that no appreciable differences in levels can 
result, though when everything is properly adjusted this link is largely static, 
and second to minimize the irregularities of the flow in J7,; caused by small 
differences in the surface tension of the water at the outside opening of 
IT,, since IT, allows the large volume of B to be used as a “shock-absorber.” 

Thus an arrangement is obtained which permits the production of a 
constant gradient of temperature which can progress regularly along a crys- 
tal X hung perpendicularly in the vessel A as indicated in Fig. 1. To com- 
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plete the ensemble, it is only necessary to mount the crystal in a holder which 
keeps it in a fixed position with regard to A, to provide connections at both 
ends with a galvanometer, and to keep the upper contact at a constant tem- 
perature, a temperature as nearly as possible equal to that of the water 
in A. All this was conveniently done with the so-called “crystal holder.” 

The bakelite tube J in Fig. 2 carries at one end the adjustable fork K, 
while close to the other, the clamp Z. Through J sealed in, runs a copper 
wire, the upper end of which is connected directly with the galvanometer. 
The other end, bent, has a silver wire (0.1 mm) soldered to it. Furthermore, 
J carries an adjustable clamping device 1J which permits one, with the aid 
of a similar one on the bridge C (Fig. 1), to fix the holder into any position 
desired with respect to A. To cool the upper end of the crystal, a copper 
cylinder JN is used, provided with two small side-tubings V, and Ne, N; being 
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Fig. 2. Diagram of the crystal-holder of the thermoanalyzer. 























connected with a water reservoir of constant temperature and hydrostatic 
pressure. In this way, a water flow is sent into N; and leaves by N2. To fix 
the upper end of the crystal X in a position within the water flow without 
applying any fatal strain to it, a thin rubber membrane O (toy balloon), 
with a small central hole burned into it is stretched over one end of N. 
After a silver wire similar to the one previously mentioned is spark-welded 
to the crystal, the hole in O is momentarily widened, the crystal is intro- 
duced and the membrane released, thus forming a tight sleeve around X 
able to withstand the water pressure in NV. The silver wire protruding now 
from the other side of N is soldered onto the copper wire P which is fixed in 
the bakelite plate Q. After a contact is thus made between P and X the 
plate Q is sealed tight onto N, thereby closing the water chamber. WN and 
X are fixed into the clamps L and K respectively and the silver wire at the 
lower end is also spark-welded onto X. The whole arrangement is then put 
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vertically into A in such a way as to bring X exactly through the central 
opening of the heater. The wire P is connected through a resistance box to 
the galvanometer, thus completing the circuit. Fig. 3 shows a photograph of 
the apparatus for thermoelectric analysis and Fig. 4 a picture of a crystal 
in its holder. It is necessary to mention that the task of mounting crystals 
without distortion is a very difficult one and requires a special apparatus 
which avoids critical stresses by facilitating manipulation. 








Fig. 3. Photographic view of the thermoanalyzer. 


After the crystal, mounted in its holder, is fixed in the vessel A, the dif- 
ferent water flows are started and the furnace-plate is heated, in general 
to a temperature 18° above the temperature of the running water. After 
thermal equilibrium is reached measurements are taken, i.e. the deflections 
of the galvanometer are measured as a function of the position of the oil- 
layer relative to the crystal. Simultaneously the e.m.f. of the differential- 
thermocouple is measured to assure the maintenance of thermal equilibrium. 
Schematically, this is shown in Fig. 5 where X is the crystal, W, is the water 
cooling the lower end of the crystal and W, that cooling the upper end, G, 
is the galvanometer measuring the thermal e.m.f. of the crystal, and G, one 
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measuring the difference of temperature between the heated oil layer O 
and the water W). 

Types of thermoanalytic diagrams. First it will be considered what type 
of curve one would expect to obtain with this kind of an arrangement. The 
obvious way of plotting the diagrams is to represent the thermal e.m.f. 
indicated by G, (Fig. 5) as ordinates and the position X of the heated oil 
layer with respect to the crystal as abscissas. If one considers first, the diagram 
that would be obtained by the thermoanalysis of a perfectly homogene- 
ous crystal, one realizes immediately that it should not show any thermo- 








WJ 





Fig. 4. Photographic view of the crystal-holder. 


electric potential (Magnus’ Law). However, certain deviations from this 
purely theoretical concept would be expected as a result of imperfect cooling 
conditions and the heat-conductivity of the crystal. Hence, a curve of the 
type (Fig. 6a) should be expected where the dotted line would represent a 
constant thermal e.m.f. (7) due to a slight difference of temperature between 
the two different end cooling systems; while the curve drawn in full with 
its deviations from linearity would be obtained under normal conditions, 
where the cooling of the ends of the crystal is not sufficient to prevent a 
slight heating by conductivity through the metal. To make this “end-effect” 
as small as possible, the crystal used has to have a small cross-section. 
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The second case to be considered is shown in Fig. 6b, which represents 
an unhomogeneous crystal consisting of two halves (I and II) which have 
a thermal e.m.f. against each other. In this case, a curve as shown would be 
obtained, where a maximum would permit the exact localization of the junc- 
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Fig. 5. The electrical connections of the thermoanalyzer. 


tion of the two halves. Superimposed upon this would be the effects pre- 
viously mentioned (Fig. 6a) though the constant thermal e.m.f. due to 
unequal end cooling would have a different value since the thermoelectric 
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Fig. 6. a. b. c. Ideal thermoanalytic diagrams of a perfect single-crystal (a), 
a perfect double-crystal (b), and a perfect triple-crystal (c). 


power of the two halves against the contact metal would be different from 
that in the case of a single crystal. Also, the end effects due to heat conduc- 
tion would be slightly asymmetrical since the two different halves are in- 
volved. The sharpness of the maximum would depend on the gradient of 
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temperature along the crystal and the extension of the heated zone, thus 
the former should be as large while the latter as small as possible. 

The third case in which the crystal consists of three sections, the second 
of which has the same thermal e.m.f. against the first and third (I: III, 
Fig. 6c) would result in a curve which is more complicated. With the heated 
oil layer at the junction of part I and II, there would be a maximum similar 
to that of Fig. 6b. Furthermore, on reaching the transition from II to III, 
the same thing would happen but in the opposite direction. Of course, the 
shape of the whole curve would depend very much on the relative sizes of 
the three sections. For a given thermoelectric difference of II against 
lil, and a given difference of temperature in the analyzing apparatus, the 
maximum and minimum would not change their size as long as II is large 
enough to contain the whole drop of temperature given by the analyser. 
If II is smaller, the maxima would decrease and finally there is a size below 
which no effect whatsoever would be observed. This threshold determines 
the resolving power of the arrangement and depends on the construction of 
the analyser and the cross-section, thermal conductivity, etc., of the crystal 
as follows: The thermal conditions along the crystal within the thermo- 
analyser are easily treated by the general equation of the temporal and local 
distribution of temperature in a linear conductor :” 


— = Q (1) 


where: @ is the temperature for (x) and (t); r the radius of the crystal; 
U the circumference of a cross-section of the crystal; S the cross-section of 
the crystal; x the length-coordinate of the crystal; 4, the thermal emissivity 
of the crystal surface; c the specific heat; p the density; and o the specific 
heat-conductivity. 

Since the movement of the oil-layer along the crystal within the ther- 
moanalyser is very slow we assume stationary condition, i.e. 06/d¢=0. 
Thus Eq. (1) becomes: 


0°0 Uh h 
a an Get nd (2) 
0x7 S a ro 
and: 
@ = A exp [— (2h/ra)'/2x] + B exp [(2h/ro)*/?x]. (3) 


For the assumption of x= ©, i.e. in case the crystal cools down so fast 
that the gradient is zero before the end of the crystal (which fact is indicated 
by the absence of the end-effects mentioned above), B and therewith the 
second term of Eq. (3) equals zero; hence we obtain: 


6=A exp [— (2h/ra)'!?x] (4) 
where A equals @nax, the temperature of the oil layer, if one assumes the 


absence of any gradient across the crystal as was already done for the entire 
crystal implicitly in Eq. (1). 


12 Enzyklopiidie d. Math. Wiss. 5, I, 181. 
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The graphical interpretation of the records of the thermoanalyzer is 
more convenient if Eq. (4) is transformed into a logarithmic system: 
In @ = In A—(a@h)!!2x (5) 


where a=2/re. If In@ is chosen as ordinate and (h'x)-' as abscissa, the cool- 
ing curve becomes a straight line as shown in Fig. 7a, its inclination (tg @) 











\ 
QA. © 
.\ c 
\ £ . 
4 \ \ i‘, J im i 
oe ee Iv AIR = WATER ~ 
\ \ . : F \ 
a : . ——— 
Wha) (When) (VE eter*) 


Fig. 7. a. b. Ideal logarithmical diagrams of a perfect double-crystal. The difference be- 
tween the thermal “emissivity” for water (/water) and air (Mair) of the crystal surface is com- 
pensated to obtain a symmetrical diagram. (Read +1 for —1 as abscissa-exponent!) 


is for the same crystal directly proportional to In§,.x and h~!, since r as well 
as o is constant for a perfect single crystal over its whole length. However 
h differs with the surrounding medium, for instance is /gir</water. Hence 
the diagram is simplified if the abscissa indicates = (hx) instead of x, since 
in this case tg @ is the same for the cooling curve above and below the oil 
layer (Fig. 7b). 
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Fig. 8. a. b.c. The geometric construction of an ideal thermoanalytic diagram of a double 
crystal (a), a symmetrical triple-crystal (b), and an asymmetrical triple-crystal (c). 8b and 
8c show each two different cases, the diagram of a large intersection (D; full-drawn lines) and 
of a narrow one (d; dotted lines). 


If the rather valid assumption be made that EZ, the thermoelectric force 
is directly proportional to the (small) difference of temperature between oil 
and water (resp. air), then equation (5) may be written: 


In E + In Emax — a! (6) 


where Emax represents the value of E for 6=@max. 
The e.m.f. produced at the ends of a double crystal (Fig. 6b) in a certain 
position (G) of the oil layer (C) (Fig. 8a) is given by the intersection with 
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the In E- & line (a). Hence it is evident that the thermoanalytic curve E = 
f(x) will follow that line in the coordination used if one moves the intersec- 
tion gradually through the oil layer; i.e. the diagram (Fig. 8a) thus obtained 
is identical with the measured curve to be expected from an ideal double 
crystal if one considers (according to the medium surrounding the inter- 
section) the two different values of / which are easy to find empirically. 

It goes without saying that the assumption of a two-dimensional oil layer 
is not justified by the actual conditions, since the oil layer is of considerable 


thickness (d), its temperature is not uniform, and d@/dx has a definite value 
as well as being a function of x9j1. 
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Fig. 9. The temperature diagram of the thermoanalyzer. The full drawn 
curve shows T=f(x), the dotted lines show dT/dx=f(x). 


However the real curves are not fundamentally different from the con-' 








structed ones as soon as the oil layer is comparatively thin and the thermal 
conductivity of the crystal is large. If those conditions are fullfilled one can 
account for the oil layer by moving the starting point of the second half of 
the diagram the distance d (dotted in Fig. 8a). This also is only a first ap- 
proximation on account of the neglection of 0°0/dx?# const. 

Fig. 9 shows the measured temperature-curve within the thermoanalyzer 
(full line). The latter can be found readily by putting a very thin thermocou- 
ple into the thermoanalyzer instead of the crystal. The dotted line in Fig. 9. 
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indicates the temperature gradient obtained by graphical differentiation of 
the first curve. It shows that very large gradients could be obtained at the 
boundary layer of oil and water and furthermore that the theoretical as- 
sumptions are not very far from the facts. 

In the case of a triple-crystal (Fig. 6c) the construction of the curve is 
more complicated though fundamentally the same. Here it is necessary to 
subtract the intersections of the In E-é lines crossing one transition from 
those crossing the other. The diagram obtained thus has ordinates 
(In E,—1n Es), while what is actually desired is In (E,— 2). Since the latter, 
however, give curves it was considered advisable from the view-point of a 
convenient graphical synthesis of the analytic curves to plot the former, 
which give the correct position of the maxima though not the correct 
amplitudes. The ordinates of the two elementary curves have to be plotted 
into the opposite direction since the sign of the e.m.f. is negative in one case. 
Fig. 8b. shows the diagrams of two different triple-crystals, one with a large 
II-section (full lines) (distance D) the other with a small section (dotted 
lines) (distance d). The thermal e.m.f. between I and II, and II and III is 
the same in both cases. 

It is apparent that the size of the maxima depends on the length of the 
heterogeneous section II as soon as d <x; (xo means the distance from the 
heated intersection to where @ is so small that the thermal e.m.f. cannot 
be measured anymore). 

The de. ase .s max is then simply given by: 


'~ - = (> — d) tgd for d < (h'?- x9) = Sp. (7) 


Moreover it is evident that the width of the maxima increases with the de- 
crease of their sharpness. 

Fig. 8c. shows another case of a triple crystal, for different lengths of 
section II. Hence where all three sections are thermoelectrically different, 
one obtains asymmetrical curves as well as an asymmetrical decrease of the 
maxima with the length of section II. 

The dependence of the size of the maxima indicated by the thermo- 
analyzer on the length of the straight section II determines the resolving 
power of the analyser in each case since it is related to the indicated e.m.f. 
in the limit by the relation: 


h\ te 
de = — Ens (=) exp [— (2h/ra)'!?x]dx (8) 
ro 
which, if the oil layer is at section IT reduces to: 
2hk\*!2 
de = — Emax (—) dx. (8a) 
ro 


The resolving limit of the analyzer is then given by a value dx which 
cannot produce a measurable de. 

The above considerations show that the analyzer, since it does not 
necessitate the slightest stress to the specimen, is a very useful instrument 
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for the detection of all kinds of invisible imperfections and inhomogenities 
in a metal rod, especially in a delicate single crystal. The size of the smallest 
imperfection still resolved distinctively was ca. 0.1 mm. It is evident that 
thin rodes (2 mm diameter in the present case) give better results than those 
of large diameter since the method integrates necessarily over the whole 
cross-section. 


Part II. ExXxprerRIMENTAL RESULTS OBTAINED FOR NORMAL AND MAGNETIC 
CRYSTALS WITH A DISCUSSION OF THE SAME 


The Thermo-analysis of normal crystals. It is obvious that the described 
method of thermoanalysis can serve the purpose of detecting and locating 
imperfections within a crystal only when the size of these imperfections is 
above the resolving power of the analyzer. Furthermore, it is only possible 
to detect imperfections which cause a thermoelectric force, though it seems 
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Fig. 10. The thermoanalytic curve of a Bi crystal grown in a glass tube. 


very probable that any imperfection due to a distortion or change in the 
orientation within the crystal causes one, the detection of which depends 
only on the sensitivity of the galvanometer used. Thus this method is an 
excellent one to detect and locate heterogeneous inclosures and local plastic 
deformations. The latter became particularly evident in the case of bismuth 
since they result in the production of twin lamellae, a phenomenon treated 
at length in other papers.!:® 

Because of the effectiveness of this method, the different means used for 
the production of single-crystals were examined as to their results, and the 
development of a method of crystal growth—a method to be used in the 
final experiments—was guided step by step by this kind of an analysis. 
First crystals grown in glass tubes by the ordinary method (Tammann, 
Bridgman) were analysed. The deviations of the thermoanalysis curves from 
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linearity were to be expected especially as the taking off of the glass cover 
caused many twin lamellae. Fig. 10 shows a diagram of such a crystal. 
The lower curve in Fig. 11 shows a crystal with one maximum (.\/) due 
toa magnetic effect (see later). After the analysis had been taken, the crystal 
was plastically deformed just enough to produce two small sets of twin 
lamellae visible under the magnifving glass. The diagram taken afterwards 








Fig. 11. The thermoanalytic curve of a crystal before (black points) and after artificial de 


formations at the points 7 —T (white points). 1/ is the starting point of the magnetic field. 


is shown as the upper curve where now two new maxima (7', 7) occur, the 
position of which coincided exactly with the location of the twin-sets. 

ig. 12 shows the diagram of a crystal of the same orientation (Ps) as 
Fig. 11 as homogenous as it could be produced by means of the graphite 
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Fig. 12. Thermoanalytie diagram of a “perfect” Bi cyrstal grown after the 


method of Goetz. Orientation: Ps; Kind “B.” 


trough, protecting atmosphere, etc., as previously described.’ Since the 
scale of this diagram is the same as that of the previous ones it is quite evident 
that the method used for the production of the crystals gave quite uniform 
results. 

Among the large number of diagrams which have been taken, it occurred 
very often that an apparently perfect crystal (as far as its examination after 
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etching was concerned) showed one or more distinct maxima. If the crystal 
then was cleaved at the point where a maximum had occurred, a microscopic 
investigation led always to the detection of the inclusion of a crystal of dif- 
ferent orientation, which orientation was usually that of a twin. Fig. 13 is a 
microphotograph of an inclosure of that kind which shows distinctly a 
twinned region within a normal (111) plane characterized by an alternating 
penetration of (111) planes through the original (111) plane bordered by 
(110) planes. The cause of these imperfections was, in general, a disturbance 





Fig. 13.  Microphotograph of an etched (111) plane (ca. 120), illuminated by polarized 


light, showing the inclosure of a region of periodical twinning. 


(vibration, shock, ete.) happening at the time when this section of the crystal 
was forming. 

The Thermoanalysis of “magnetic” crystals. The described method of 
thermoanalysis was applied to crystals, one half of which was grown without, 
the other half within a magnetic field. Jt was found that the border between the 
magnetic and the normal half of the same crystal was the origin of a thermo- 
electric force. Furthermore, the thermoanalysis gave evidence that the origin 
of this force is exactly at the point where the magnetic field was applied. 

Nevertheless, it proved to be very difficult to measure this force quantita- 
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tively and not even the sign of the force could at first be reproduced although 
each curve showed a maximum if it was not fogged by other maxima due to 
distortions of the kind mentioned in the previous chapter. 

To obtain a higher degree of reproducibility, certain definite conditions 
were imposed, both on the growth and on the analysis of the crystals. In 
the case of the former, four factors were considered in detail; orientation, 
method of growth, impurities in the metal, and field strength. For the latter, 
two were considered: temperature of the analyser, and annealing. 

The actual resuits obtained indicate that at least for the first group of 
factors the thermoelectric effect is a very complicated function. Despite this 
complexity, however, duplication of results, for any given set of conditions, 
with less than ten percent deviation was always possible which indicates 
that all the principal variants were being controlled. 

In presenting the results, the effect of each of the factors on the thermo- 
electric e.m.f. will be considered separately, the stationary values of the 





Fig. 14. Scheme of the three principal orientations used. The hatched plane indicates (111), 
the arrow shows the direction of the field lines. 


other factors being given in each case. Up to date, no attempt has been made 
systematically to evaluate the effect for a large number of values of each of 
the factors involved, thus the results given can be considered as merely 
illustrating the type of variation to be expected. 

For all the curves to be presented unless specifically stated, the tempera- 
ture difference in the analyzer was 18°C measured above tap water tempera- 
ture of 19° to 21°C. 


TaBLE I. The three orientations of the crystal. The third column entitled “vector” shows the 
direction of the heat flow and the electric current during the thermoanalysis. 





Trigonal axis to (111) plane to 
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The ordinate scale for the curves unless designated otherwise is 7.8 X 10-6 
volts per division. It indicates positive volts for increasing values. The 
units of the abscissa (length of the crystal) are centimeters. 

The solid line, mid-way along the abscissa scale on the curves, indicates 
the position of the furnace at the intersection between magnetic and normal 
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portions of the crystal.* The dotted lines indicate the extent of the oil layer 
to each side of the furnace. 

The Effect as a function of the orientation. As stated already in a pre- 
vious paper,! there are three different orientations the crystal can have with 
regard to the rod, or to the direction of the lines of force, which are essentially 
different from each other.** Fig. 14 shows the trough G dotted, and the crys- 
tal white, the latter cut parallel to its main cleavage plane. These three ori- 
entations are given in Table I. 
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Fig. 15. The thermoanalytic curves of three principal orientations of crystals 
(kind “B”) grown discontinuously. Field: 21,000 Gauss. 


Fig. 15 shows analysis curves typical for crystals of these three orienta- 
tions. The method of growth employed was the discontinuous (see Part I), 
the metal of purity “‘B” (see later), and the field strength 21,000 Gauss. 


* The thickness of the center line is misleading in so far as the furnace plate was 1 mm 
thick and since it was impossible to measure exactly its location within the oil layer. 
** The possibilities of azimuthal variation were neglected in these experiments. There 
are: 
two for P,; with regard to field and vector; 
two for P, with regard to vector; 
two for P; with regard to field; 


which differ essentially from each other. 
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The diagram shows that the thermal e.m.f. depends very much on the 
orientation of the principal axis of the crystal to the field lines. The effect 
is very large for P; and very small for P2; it is even very probable that the 
latter effect does not exist at all, and that the deviations of the P.-diagram 
are due to inaccuracies in orientation. 

The thermal e.m.f. is positive for P; and negative for P;.*** Furthermore 
it is apparent from Fig. 15 that the thermal e.m.f. is different from what 
should be expected, (see Part I) since the maximum of the P3-curve is not at 
the intersection with the center line, i.e. within the region where @ is largest 
(see Fig. 9), but at the border between oil and water, where d6/dx is a maxi- 
mum, a fact which seems to be of fundamental importance. 

The size of the thermal e.m.f. of the P3 curve seems remarkable (7.8 x 
10> volt or 4.33 x 10-* volt. degree~'), if one calculates the change in orien- 
tation of the crystal necessary to produce this e.m.f. due to the thermoelectric 
anisotropy of the Bi crystal. Taking the data from Bridgman’s last paper" 
one finds that change would have to be 21°**** though it is certain that there 
is no larger change of orientation at the intersection than 0.5°, which fact 
makes it quite impossible that the thermal e.m.f. is just caused by a change 
of orientation at the intersection. 

The Effect as a Function of the method of growth. The methods of 
growth employed have been described in detail in Part I of this paper. They 
were the so-called continuous and discontinuous methods. Fig. 16 shows 
curves obtained from crystals of P; orientation produced in these different 
ways. “1” represents the former, “3” the latter; in both these cases, no appre- 
ciable residual field was present over the so-called unmagnetic portion of the 
crystal. “2” on the other hand, represents a crystal of continuous growth 
with the residual field present during the growth of the first portion of the 
crystal. 

It was entirely unexpected to find that the thermal e.m.f. depended on 
the method of growth to such an extent as shown in Fig. 16. The fundamen- 
tal difference between “3” and “1” is that there exists only one maximum in 
the former (the small second maximum at “db” has a different cause and will 
be discussed later) and a maximum and a minimum in the latter. The sign 
of these maxima is opposite in “1” and “3” but their positions coincide again 
with the extreme values of d@/dx in the thermoanalyzer, for “1” the thermo- 
analytic curve is even an exact repetition of the d0/dx curve in Fig. 9. The 
case of “2” shows the large effect of a small field strength which seems to 
suppress entirely the first maximum. 


*** The thermal e.m.f. is called positive, if the end of the “magnetic” part of the crystal 
had a potential positive with regard to the end of the normal half, if the temperature of the in- 
tersection was higher than the temperature of the ends. 

**** This diagram (Fig. 16, page 384) shows the thermal e.m.f. as a function of cos? ¢ 
(angle of orientation). The units of the ordinate are 3.10 X10~* volt for a difference of tem- 
perature of 68°. Since the curve is no straight line it is necessary to measure at an orienta- 
tion corresponding to P3, i.e., starting from cos? ¢=1.0. Since the e.m.f. unit in that diagram 
equals 3.10 X10-*/68°=4.56 X10-* volt. degree“ our e.m.f. corresponds to 4.33/4.56=.95 
units. For the P;-orientation (cos ¢=1) .95 ordinate-units correspond to 1.3 units of the 
abscissa, i.e. cos ¢= (0.87)¥2, ¢—21°. 
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The effect as a function of the impurities. It became evident that the 
magnitude of the effect depended on the purity of metal used to a very large 
extent. The crystals were therefore prepared from four kinds of Bi, obtained 
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Fig. 16. The thermoanalytic curves of three different types of growth: 

1. represents the continuous method without residual field. Field: 13,800 Gauss; kind 
“B;” orientation: P3. 

2. represents the continuous method with residual field. 

3. represents the discontinuous method. 


from different sources which were designated as “A”, “B”, “C”, and “D” as 
shown in the following table: 


TABLE III. Sources of the bismuth used. 














Designation Characterization Source 
“— Bismuth C. P. Brown Corp., Philadelphia, Pa. 
— Bismuth purissimum Hartmann & Braun, Frankfurt a, M. 
m. “ys Electrolytic bismuth ” as 
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| >" Bismuth “Kahlbaum” Kahlbaum A. G., Adlershof. 
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Fig. 17 shows curves of the P; orientation made from these four different 
kinds of Bi by the discontinuous method, whereas Fig. 18 represents the cor- 
responding crystals of continuous growth. 

Since from this work as well from papers of Bridgman™ and Kapitza" it 
was well-known how large the influence of small impurities on the electric 
effects of Bi can be, an accurate chemical analysis of the different kinds of 
bismuth was made. The results obtained showed that a very small trace of 
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Fig. 17. The thermoanalytic curves for P; crystals of four different degrees of 
chemical purity, grown discontinuously; field: 21,000 Gauss. 


silver was present in “A” and “B” though the amount was too small to 
measure accurately. Since this indicated that the impurities involved were 
only present in extremely small amounts, a spectral analysis seemed neces- 
sary. Dr. R. M. Badger was kind enough to perform this type of analysis 


18 P, W. Bridgman, Proc. Am. Acad. 63, 351 (1929). 
144 P, Kapitza, Proc. Roy. Soc. A119, 358 (1929). 
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on our different metals. The amount of any impurity was measured by com- 
paring the intensities of its different lines as observed for the specimen, with 
the intensity of the lines obtained from these metals by themselves in a spark. 
The bismuth electrodes were prepared by putting the metal into Pryex tubes 
after it had been ascertained that the glass did not contain any of the metals 
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Fig. 18. The thermoanalytic curves for P; crystals of three different degrees 
chemical purity, grown continuously; field: 13,800 Gauss. 


involved in the analysis. The results are given in the following table, where 
the numbers give the approximate relative intensities of the indicated lines 
of the different impurities. 

In general, it can be said that the amount of impurity is probably on the 
whole, smaller than 0.2 percent in the worst case which is “B”. It is smaller 
for “A” which contains less Pb. “D” differs from “B” in its smaller content 
of Pb, though in its Ag content, it compares badly whereas “C” is by far the 
purest metal. 
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TABLE IV. Spectroscopic analysis of Bi samples. 
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Sample Ag Cu Sn Pb Zn Sb 
Le) 

_ * ™~ rN Lao) io 2) * nN wy NN * i) NN wy ™~ N ie 9) oo 
co 0 + rico S| HR O© HH ~ CO wel a NN es om AN 
N * NS NIN i> 2) | Oo Oo N io 2) o i> @) ™ * NN —N ~ wy 
Lae) * * m1 N 7 ~*~ Lae) * N “nN i b> al * * Lae) N N 

A 4 1 1 ? 5 1 ? ? 

B, 5 2 1 ?}10 4 3 

Be 5 2 1 ?1/10 4 2 

C, 1 5 3 ? 1 ? 

C, 1 3 3 1 ? ? ? 

D, 7 3 1-0 1 6 

Dz 7 #4 1 5 ? 























Thus it is apparent from the diagrams that the purest metal shows the 
smallest effect and also that the shape of the curve is but little affected. 
However the secondary maximum (“b” in Fig. 16) depends apparently on 
certain impurities, since it is very large for “A” (Fig. 17) and extremly small 
for “C”. The method of continuous growth shows very little difference 
between “A” and “B” and a large one for “D”. It is however doubtful 
whether the latter measurement is a correct one, since it is based only on a 
few crystals. 
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Fig. 19. The size of the maximum thermal e.m.f. as a function of the field 
strength for P; crystals, kind “B,” grown discontinuously. 


The effect as function of the field strength. Although it is very im- 
portant to investigate what influence the strength of the applied field has 
upon the thermoelectric effect, it appeared very difficult to obtain reliable 
results, since the effect does not depend merely on the orientation of the 
crystal but also on the method of growth, the degree of impurity, etc. Thus 
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it has not yet been possible to obtain results of sufficient generality in spite 
of a great number of experiments. It seems however quite safe to state that 
there exists a large influence of the field strength upon the effect for crystals 
of the P3-orientation. (Other orientations have not yet been investigated.) 
The influence of weak fields is very large whereas strong fields (20,000 Gauss) 
decrease the effect considerably in the case of the discontinuous method of 
growth. 

Fig. 19 shows the variation with field strength obtained for crystals grown 
with the P; orientation, composed from a series of curves of the type of 
Fig. 15. 

The field which produces the maximum thermoelectric effect seems to 
depend on the degree of impurity. 
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Fig. 20. The size of the maximum thermal e.m.f. as a function of the temperature (The 
hypothetical transformation of @ into 8-bismuth.) P3-crystal, grown discontinuously, kind 
“B,” field: 21,000 Gauss. The dotted curve gives dE/dT=f(T). 


The effect as a function of the temperature. Fig. 20 shows the thermo- 
electric effect between the magnetic and the normal half of a crystal as a 
function of the temperature of the analyzer and its differential curve (dotted). 
The curve was taken by first finding the exact position of the maximum e.m.f. 
with regard to the crystal (crystal of the type shown in Fig. 15) and then 
keeping the oil layer in this position while its temperature was gradually 
increased. 

This experiment performed on many crystals of different types showed 
that the temperature function is independent of the impurity as far as the 
general shape of the curve is concerned. It seems quite interesting that the 
differential curve dE/d@ shows a large discontinuity between 75° and 90°C. 
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Since it is well known that this type of curve is obtained in case of an allo- 
tropic transformation (as for instance typically represented by the thermo- 
electric effects of pure iron") the shape of the curve obtained indicates that 
the discontinuity in that region of temperature affects the magnetic half of 
the crystal in a different way than the normal.* 

In order to decide whether or not the magnetic state is just a more or 
less unstable allotropic state caused by strains, such as magnetostriction for 
instance, specimens showing large effects were reinvestigated after being 
annealed. The temperature of annealing was chosen sufficiently high, above 
90°C, to see whether the transition into the normal state would occur at all. 
Since this annealing temperature lies above the point of the mentioned 
“allotropic” transformation, it is probable that the instability of the magnetic 
state should increase, or the transformation speed into the normal state 
should be larger. 

Hence a crystal (P; of “B”) was analyzed the day it had been grown, it 
was then annealed for 16 hours at 100°C and analyzed again a month later. 
No change of any kind could be detected which indicates that the magnetic 
state is either perfectly stable or that it has an extremly slow speed of trans- 
formation. It is certainly not affected by passing the critical temperature of 
normal Bi. (75°-90°C) 

This experiment shows also that the effect of “aging” Bi crystals does not 
occur in specimens of good crystallographic perfection, either for normal or 
for magnetic crystals. 


* The question of the existence of an allotropic state of Bi is still unsettled; Cohen" and 
Cohen and Moesveld'’ concluded from their pyknometric and dilatometric measurements 
the existence of an a and a 8 modification, of which the a-state is stable below 75°C having a 
larger density than the 8-state. However those authors found other transformations at 81° 
and 90° where the volume of the metal changed suddenly according to its thermal history. 
Wiirschmidt!* found also the transformation point at 75° indicated by a large contraction, 
whereas the temperature coefficient of the conductivity of Bi as investigated by Bridgman" 
and Holborn* does not indicate such transformation at all.—Furthermore the thermoelectric 
effects as a function of temperature of Bi single crystals were investigated by Boydston;* 
his E=f(@) curves show a slight indication of an irregularity in this region of temperature 
which becomes much more evident in the case of the Peltier effect curve: I= T dE/d0=f(@). 
A recent investigation of the Peltier and Thomson effect by Fagan and Collins” would have 
cleared this question, if the authors had extended their measurements into this region of tem- 
perature. The exact coincidence of the critical temperatures in our measurements with the 
transformation-point of Cohen, Moesveld, and Wiirschmidt seem to render the existence more 
probable. It may be however another example of “pseudo-allotropy” due to traces of im- 
purities as shown recently in the case of extremely pure Zn by Guertler and Anastasiadis.* 
% A, Goetz, Phys. Zeits. 25, 562 (1924) ; 26, 260 (1925). 

6 E, Cohen, Akad. Amsterdam Verls. 23, 1224 (1914/15). 

17 E, Cohen and A. L. Th. Moesveld, Zeits. f. Phys. Chem. 85, 420 (1913). 

18 J. Wiirschmidt, Verh. d. Dt. Phys. Ges. 16, 799 (1914). 

19 P. W. Bridgman, Proc. Am. Acad. 52, 636 (1916/17). 

20 L. Holborn, Ann. d. Physik (4) 59, 152 (1919). 

21 R. W. Boydston, Phys. Rev. 30, 911 (1927). 

2 H. D. Fagan and T. R. D. Collins, Phys. Rev. 35, 421 (1930). 

#3 W. Guertler and L. Anastasiadis, Zeits. f. Metallkd. 21, 338 (1929). 








1778 A. GOETZ AND M. F. HASLER 


DISCUSSION 


Generally speaking the above results are far from sufficient to start an 
explanation of the effect on a theoretical basis. However the results seem to 
us complete enough to describe the phenomenology of the new effect approxi- 
mately. 

In the first place, it becomes quite apparent that the forces acting at 
the moment of crystallization affect the physical properties of the final crys- 
tal, though neither orientation nor any other macroscopic quality of the 
crystal is changed. 

The effect as a function of the orientation of the crystal with regard to 
the field seems the most obvious of all the different variations of the experi- 
mental conditions. Here we obtain the largest thermal e.m.f. if the principal 
axis is normal to the field lines and parallel to the heat flow, i.e. the direction 
of the thermoelectric current (P3); whereas the effect is extremly small if the 
axis is parallel to the lines of force and normal to the flow (P2). 

If one considers the anisotropy of the diamagnetic susceptibility of Bi 
crystals as measured recently by Focke® a relation between the two effects is 
quite evident, since in case of P; the crystal is in a position where the direc- 
tion of the largest diamagnetism is parallel to the lines of force, and in the 
case of P, the most “paramagnetic” direction is parallel to the field. Hence 
there is a difference between the content of free energy in these cases, for the 
stability of the latter orientation is larger than of the former. If thus a crys- 
tal is grown first under normal conditions and then in a field the stability does 
not change in case of P, but is does considerably in case of P;. One knows 
from other observations that a crystal prefers the P: orientation when form- 
ing its first center of crystallization in a transverse field and also that the 
field does not even change the opposite orientation if the latter is enforced 
on it by inoculation with a seed crystal. The fact that the point where the 
field (i.e. the instability) started is indicated by a thermal e.m.f. shows that 
the crystal undergoes a change in spite of the constancy of the orientation. 
Thus the question about the nature of this change becomes quite interest- 
ing. 

Regarding the shape of the thermoanalytic curves it seems as if the curves 
of discontinuous growth (Fig. 15) were due to a double-crystal, and the curves 
of continuous growth (Fig. 18) represent a triple crystal, as has been discussed 
in Part I of this paper. The interpretation however of the latter diagram 
as being produced by a triple-crystal seems to be inadequate in view of the 
fact that the length of sandwiched section (d in Fig. 8b and 8c) would have 
to coincide with the thickness of the oil layer to produce the maxima and 
minima at the borders of the layer, which coincidence would seem highly 
improbable; i.e. the maxima and minima of the diagrams are too close to- 
gether to be sufficiently resolved if caused by a section smaller than the thick- 
ness of the oil layer. Furthermore those points are not situated at the 
hottest part of the oil layer of the thermoanalyzer (Fig. 9), which means that 
the e.m.f. does not reach its maximum if the intersection of the crystal passes 
Omax- This concerns also the simpler curves of discontinuously grown crystals. 
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It is very remarkable that the maximum of the thermoanalytic curves of P3 
crystals of that kind does not coincide with the position of the maximum of 
temperature, it does however with the maximum of the d@/dx curve (Fig. 9), 
i.e. at the border of water and oil. The only simple interpretation is thus that 
the thermal e.m.f. must depend on the temperature gradient or the density of the 
heat flow crossing the intersection between the normal and the magnetic half 
of the crystal. 

It is apparent that such effect cannot be caused simply by a double or a 
triple crystal, since the normal thermoelectric effect between two different 
metals or orientations of a crystal depends on the difference of the tempera- 
ture and is quite independent of the gradient. 

Although the assumption of such a type of a thermoelectric effect seems 
to be quite uncertain in view of the fact that the passing of the opposite gradi- 
ent (air-oil, see Fig. 9) is much less expressed by the thermoanalytic curves of 
discontinuously grown crystals, such interpretation gains some probability 
in the case of continuously grown crystals (Fig. 18). Here we have curves 
which are the exact duplicates of the d@/dx-curve in Fig. 9. The sign of the 
effect is opposite in these two cases as follows: If the heat flow passes the inter- 
section from the magnetic to the unmagnetic half of the crystal, the magnetic 
half is positive for continuous growth and negative for discontinuous growth, 
it is opposite if the heat flows the other way. In case of the discontinuous kind 
of growth the flow in the latter direction seems to be quite uneffective. 

It is apparent that the influence of the impurities makes the effect very 
complex since the amount as well as the kind of impurity affects the thermal 
e.m.f. considerably as is to be expected from a consideration of the depend- 
ence of the anisotropy of the electric qualities of the crystal on the impurities 
present. Hence the mentioned instability of a growing crystal for certain 
orientations to the field should increase with an increasing anisotropy due to 
impurities and one should expect a kind of proportionality of the thermal 
e.m.f. with the impurity. This however is misleading, since it seems highly 
probable from the behavior of the metal “C” that an ideally pure Bi would 
not show any effect at all though the electric anisotropy still exists. Thus 
one is forced to assume that the whole effect is due to disturbances of the odd 
atoms in the Bi lattice, which disturbance is the more distinct the better 
the crystal is grown. The disturbance for a given amount and kind of odd 
atoms depends apparently on their arrangement and is easily affected by 
external influences during the act of crystallization. 

Very little can be said as yet about the actual nature of the disturbance 
though the fact that an extremely small amount of impurities has such a large 
influence seems to indicate a surface effect, i.e. an effect due to a formation of 
mono- or bimolecular layers of odd atoms sandwiched between perfect sec- 
tions of the crystal. Such arrangement similar to the “II-planes” of the mosaic 
structure of Zwicky’** would require an extremely small amount of impurity 
to affect the electronic relations inside the lattice. In case of such an inter- 
pretation one has to assume that the arrangement of the odd atoms is changed 
under the influence of the magnetic field. It is quite plausible that such a 
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change differs with different crystallographic planes as well as that its effect 
works differently upon an electric current, which might explain the difference 
between the thermoanalytic curves of P; and P;. (As has been mentioned 
already, these two orientations do not differ with regard to the field though 
they do with regard to the direction of the thermoelectric current and the 
heat flow.) 

The arrangement of the odd atoms within the crystal must depend 
furthermore on the growth of the crystal, since any disturbance during the 
formation will affect the structure of the impurity layers to a certain extent 
and therewith the final influence of the magnetic field. This would explain 
the difference between the thermoanalytic curves of continuous and discon- 
tinuous growth. The highly metastable state of the configuration of the odd 
atoms within the II-planes during the formation of the mosaic structure is 
probably the cause of the extreme experimental difficulties in obtaining 
reproducable results. 

Before closing, some comparisons with Kapitza’s work'* may be made. 
Concerning his observations of the change of the resistance of Bi crystals 
in a strong magnetic field, it is to be expected from the above hypothesis 
that the Bi which produced an extremly small effect (K.’s crystal “B” in 
his Fig. 13, p. 407, l.c.) would produce in our case one of medium magnitude, 
since his metal corresponds to our “D” whereas his crystal “A” which pro- 
duced large effects was of metal of the same source as our “C” giving—for 
our work—the smallest values. Since the influence of an external field upon 
the change of the resistance depends on the size of the disturbances within 
the lattice, our experience with “D” and “C” compares very well with Kapit- 
za's result, his effect being fogged by impurities and ours depending apparently 
wholly on their presence. 

With regard to the orientation effects, a few comparisons can also be 
made. Being limited to only one direction of both the field to the rod, and the 
thermoelectric current to the rod, only one kind of Kapitza’s experiments 
can be considered for comparison, i.e. where the field is normal to the cur- 
rent. Neglecting the azimuthal effects which seem in our case small, at 
best, comparison shows that the largest values obtained of the ratio R/Ro is 
for our so-called P; orientation, while the value obtained for P; was only a 
fractional part of that for P3, both facts being analogous to our results. In 
the case of P, however, a larger value of R/Ry than the P; case was obtained 
which might be interpreted as suggested above, because our effect shows 
no appreciable influence of the odd atoms upon this orientation. 


CONCLUSION 


Although these experiments were begun several years ago during which 
time a very large number of measurements have been taken, it is not yet 
possible to give more than a rather vague hypothetical explanation as out- 
lined above. Since the experimental difficulties proved to be very large the 


*% F, Zwicky, Helv. Phys. Acta 3, 269 (1930). 
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actual progress has been slow and the present paper cannot be more than a 
first report about the new effect. 

Systematic measurements concerning the electrical conductivity, the 
specific gravity, and the magnetic susceptibility were performed and are 
still going on, as well as a detailed x-ray analysis of the “magnetic” crystals. 
Papers concerning these subjects will be published in near future. 

In conclusion, the authors wish to express their great obligation to Dr. 
R. A. Millikan for the interest he has taken in the work, to Dr. R. M. 
Badger to whom they owe the spectroscopical analysis of the bismuths, to 
Mr. J. Pearson for his helpful technical advice and to Mr. A. Focke for his 
general assistance. 
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THE SECOND VIRIAL COEFFICIENT FOR GASES: A 
CRITICAL COMPARISON BETWEEN THEORETICAL 
AND EXPERIMENTAL RESULTS 
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(Received November 5, 1930) 


ABSTRACT 

On the basis of a recent theory developed by London the second virial coefficient 
of the following gases is computed: He, Ne, A, He, Ne, Or, COz, NH; HO. For the 
latter two, which have a dipole moment, the contribution of the dipole interaction is 
also calculated. The results are compared graphically with experimental values. 
Except in the case of He and Hz, where the disagreement can be accounted for by the 
presence of zero point energy, the accordance is as good as is compatible with the 
inaccuracies involved in the calculation. 


N CLASSICAL mechanics it was necessary to attribute the attractive 

forces between molecules, demanded by the empirical equation of state, 
to the action of permanent or induced electric poles. Quadrupoles of con- 
siderable moment had to be assigned to molecules such as Hz and the rare 
gases, which were known to carry no permanent dipoles. On the basis of the 
old atomic models the required assumption of quadrupoles was entirely 
reasonable, but it is inconsistent with the more uniform distribution of 
molecular charges calculated by wave mechanics. Moreover, molecular 
quadrupoles are not yet accessible to measurement, so that a theory based on 
their presence is unsatisfactory for its lack of verification. This difficulty was 
very successfully removed when Ejisenschitz and London! recognized the 
significance of attractive forces appearing in the second order perturbation in 
the interaction problem of two H atoms. They showed that Van der Waals’ 
constant a derived from these forces had the correct order of magnitude. 
In a recent paper London? made a more comprehensive study of these 
forces (arising from what he terms the “dispersion effect”) and discussed in 
particular their relation to two other types of attractive forces existing be- 
tween molecules having electric moments (resulting from Keesom’s effect of 
“alignment” and Debye’s induction effect). He also calculates the con- 
tribution to the molecular interaction energy of the dispersion effect and 
finds that this is for many gases the dominant constituent of the interaction, 
resulting in a van der Waals’ a which is in good agreement with the value for 
a computed from critical data. This is indeed a most interesting and signifi- 
cant observation which goes far in confirming the validity of the theory. Now 
it is well known that van der Waals’ equation with a independent of the 

1 R. Eisenschitz and F. London, Zeits. f. Physik 60, 491 (1930). 


2 F. London, Zeits. f. Physik 63, 245 (1930). 
3 Cf.2 for literature. 
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temperature does not represent the empirical equation of state, nor does the a 
value computed from critical data coincide of necessity with the factor of 
1/T in the equation of state. It is the purpose of this communication to make 
a closer comparison between the predictions of London’s theory and experi- 
mental data, and also to include in this comparison the results of a recent 
paper* dealing with forces between dipole molecules. We shall express the 
second virial coefficient, B, as a function of the temperature and examine its 
correspondence to the empirically known values of B for various tempera- 
tures. The results can best be shown by graphs. 

It will appear that the calculation of B is of necessity very crude because 
idealized assumptions and lack of accurate knowledge of certain quantities 
involved; hence one is justified to look for qualitative agreement only. The 
method used is essentially the same as that employed by London to compute 
a and is subject to the same errors. 


Non-PoLar GASES 


It may be recalled that B is the coefficient of 1/V in an expansion of 
bV/RT as a power series in 1/V. It measures the deviation from Boyle’s 
law for small concentrations (large V) and is given analytically by: 


B= en [ (e~e/AT — 1)r2dr, (1) 
0 





where € is the interaction energy of two molecules expressed as a function of r, 
the distance between the molecules. If N is Avogadro’s number, B is in ce 
per mole. The result of London’s® investigation concerning the dispersion 
effect is: 

3 a®V, 3 a®V; 

~- <-es— (2) 

4 


ro ~ 4g 7 








where a is the polarizability, V, the excitation potential, and V; the ionization 
potential of the substance. An accurate theoretical expression for € involves 
the energy values and matrix elements of the molecules and is, therefore, 
in general difficult to evaluate. At r<d the interaction energy ceases to be 
given by (2), € merges into the rapidly rising exponential curve characteristic 
of the valence forces.* Here the assumption will be made that 


' e=+o forr<d (3) 


corresponding essentially to the supposition that the molecules have a di- 
ameter d and are impenetrable. Refinements, such as introducing a depend- 
ence of d on T, as would be necessary to take account of the finite slope of 
the exponential curve at r=d are possible but can be omitted in view of the 
general lack of precision of the premises on which we are compelled to make 


*H. Margenau, Zeits. f. Physik 64, 584 (1930). 
5 Reference 2, p. 256. 
6 See Fig. 1 in the paper by Eisenschitz and London. 
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the calculation. We shall therefore assume an empirical value for d, inde- 
pendent of T. Moreover we shall put 
3 a’V; (4) 
(a oe eo! _— 
4 f° 





introducing thereby a considerable error. (There is reason to suppose that 
in general € is nearer the upper limit in (2).)? Expanding the exponential in 
(1) and performing the integration using (3) and (4) we obtain 


2 Nd = 6; 6; 2 A; 3 
3 = i—0./3— — 9.4X 10-+(— —1.4X 10-+(—*) 
3 T T T 








0, 4 6; 5 
— 1.9 10 (= —2.2xX 10-+(=) 
~) t 
~ 2.2 1@*1—) —.--. 5. 5 
G § (8) 
where 
a? V; 
6; a 
d® k 
If we put 
2xNd* 
Wiles 
we have 


4r?N?2 a?V; 


OR Bo? 








By can be determined from the empirical B curve and is essentially the same as 
van der Waals’ b. The values of Bo used in the computation of B are taken 
from Beattie and Bridgeman’s® work, except for O2 and COs. For these two 
gases the experimental and the theoretical B obtained from the data listed 
by the latter authors proved to be peculiarly discrepant, more so than could 
be explained by the crudeness of the calculation. It may be observed, how- 
ever, that while for all other gases listed here Beattie and Bridgeman’s By 
is in the neighborhood of the usual b-values (as computed from critical data), 
it is considerably higher for O2 and CO, (46 and 105). Also, in the case of He, 
Ne, A, N2 their constants are derived from isothermals up to 400°C, whereas 
for the two gases in question only temperatures up to 100°C were used. If 
now it is remembered that By is the limit of B for high temperatures one could 
be led to suppose that though Beattie-Bridgeman’s Bo, together with their 
other constants, describe the behavior of O2 and CO: very well over the range 
of lower temperatures, it may not actually be the limit which is here required.® 


7 See reference (1), action of the continuous spectrum. 

8 Cf. for instance Zeits. f. Physik 62, 95 (1930). 

® At 258° Beattie and Bridgeman’s constants for CO, give a value of B which is, compared 
to Amagat’s measurements, about 50 percent too small. 
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If this supposition is erroneous the discrepancy spoken of would have to be 
admitted as unexplained. We have in the following assumed By values for Oz 
and COs: equal to the b-values from critical data, which are in agreement with 
our knowledge about the size of molecules derived from other sources (vis- 
cosity). If they are chosen the experimental and the theoretical B-curves 
agree quite well. 

Table I contains the experimental data on the basis of which the computa- 
tion of B was made. The ionization potentials of some of the gases may be 
in considerable doubt, but in view of the other uncertainty concerning the 








TABLE I. 
aX 104 By J ‘ (v olts) 0; 
: He 0.20 14 24.5 92 
Ne 0.39 20 21.5 150 
A 1.63 39 15.4 500 
He 0.81 21 16.4 440 
No 1.74 50 17 385 
O. 1.57 32 (46) 13 580 (280) 
CO, 2.9 42 (105) 10 880 (180) 











value of V to be used (Eq. (2) ) this need not disturb us particularly. It is 
interesting to observe the magnitude of 6;. The calculation of van der Waals’ 
a'° neglects higher powers of 6,/T than the first, which seems unsatisfactory, 
even if due consideration is given to the large errors in our assumptions. 
CO., though having a small dipole moment (u=0.06 X10-"*) is included in 
the list of non-polar gases. 


20 











Fig. 1. Neon. 


The curves in Figs. 1-7 show the results, the ordinates being in cc per 
mole, the abscissae in absolute degrees. Except for CO. (where the experi- 
mental B-values were taken from Amagat’s diagram) the experimental points 
were computed from Beattie and Bridgeman’s constants. The broken curves 


10 F, London, reference 2. 
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Fig. 3. Helium. 
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Fig. 7. Carbon dioxide. The broken curve was computed on the assumption B = 105. 
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in Figs. 6and7 represent the values obtained by using By =46 and 105 re- 
spectively. The agreement is seen to be good for all substances with the 
exception of He and H:2, where the large deviation can be accounted for 
qualitatively by the zero point energy of the interacting molecules. As 
London has pointed out, this is negligible for molecules of larger mass, but 
would, in the case of He and He, cause the theoretical curves to be displaced 
toward the experimental ones. 


PoLAR GASES 


It has been shown" that for sufficiently high temperatures (k T> h?/872A) 
the virial coefficient of polar molecules is given by 


N?2 ea) 
B= — | i) (eV/AT — 1)r2drdQ (6) 
8x Joa 9 


where the first integration is extended over all orientations of the two mole- 
cules and V is their classical potential energy in terms of r and the angles. 
For dipole molecules capable of polarization we find, if u is the dipole moment, 


oe aa 
V= - — 2 cos 6; cos 62 + sin 4; sin 62 cos (¢1 — ¢2)| 
r 


9 


pa 
on — cos? 6; + 3 cos? 62, + 2) neglecting small terms. (7) 
r 


The B computed from (6) will have to be added to the B resulting from the 
dispersion effect, but in doing so it must be remembered that the integration 
over r from 0 to d has already been counted, and the relevant part of (6) 
results from fa”. Expression (6) has been calculated by Falkenhagen.” He 
finds for several dipole gases, among them NH; and H.O, good agreement 
between experimental data and the results of (6). This is very surprising from 
our present point of view, for it would seem to indicate that, if the dispersion 
effect also were taken into consideration the total B would be much too large. 
This difficulty, we feel, resolves itself on closer examination. It is impossible 
to check Falkenhagen’s calculations in detail because the numerical values 
used in them are not all stated. (If we use the values of d derived from the 
By of Table II, the agreement found by him is considerably disturbed.) We 
also note that his potential energy expression contains a polarization term 
(the second term in the brackets of (7) ) twice as large as ours, which cir- 
cumstance we are inclined to attribute to a numerical error in his Ge. This 
tends to enlarge the contribution of (6) to B. 


11 H, Margenau, Zeits. f. Physik 64, 584 (1930). In Eq. (1) of this paper a bar was omitted. 
It should read, of course, 


B=2rN f (eV (2, @) kT _ 1) R2GR. 
0 


13 H, Falkenhagen, Phys. Zeits. 23, 87 (1922). 
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In evaluating (6) we remember that 


[ a = i) f f [si 6,46, sin 62d62d¢ doo, 
dQ 


and obtain, after expanding the exponential, for the dipole contribution to 
the total B, 








2r Nd* { 2 62\2 6.\3 
B’ = - —+ 2s-— + [0.33 + 0.7327] (=) + [0.272 + 0.3525] (=) 
3 lt T T T ; 
62\4 62\° 
+ [1.3 x 10-* +-0.212t] (S + 2.3 xX 10-%(=) 
T T 
+5.2 x 10-4 (=) (8) 
* T f° 
Here 
2 nN yp? 
9» = kil = - rl 
kd® = 3k Bo 
a 2xN 
" d3—- 3Bo 


There is but little sense in using more terms since (1) would probably lose 
its validity for temperatures in which these terms become significant. The 
only dipole gases for which a comparison can be made at present are NHs; 
and H:O on account of the lack of experimental data for others. Table II 
shows the data from which the B values of these two gases were computed. 














TABLE IIT. 
a X 104 Bo V5 uX10'8 6; 62 4 
NH; 2.41 34 11 1.50 880 610 0.085 
HO 1.48 34 13 1.84 470 920 0.055 








By for NHs3 was taken from Beattie and Bridgeman’s work and the experi- 
mental B computed from their constants, while for HO it was necessary to 
take the critical data b, using for comparison experimental data by M. Jacob™® 
from which Falkenhagen™ has calculated the second virial coefficient. 

In Figs. 8 and 9 the experimental curves are plotted together with values 
obtained from (a) (5), (b) (5) + (8). It is seen in both cases that the dipole 
effect is a very essential part of the interaction causing B. For NH; the agree- 
ment between B,,, and Baisp+dipote is quite good, showing that the two effects 
here considered suffice to explain the molecular interaction. For HO, how- 
ever, it is evident that they do not suffice, and that we must look, for quad- 


13 M. Jacob, Zeits. f. Ing. 1912, p. 1980. 
4 Reference 12. 
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rupole interaction to fill the space between the two lower curves of Fig. 9. 
It is plausible to assume that water molecules have a quadrupole moment, 
and one can without difficulty compute it from the last figure, for Keesom’s 














Fig. 8. Ammonia. Curve a represents values of B due to the disperson effect, curve 6 those 
due to dispersion plus dipole effect. 


theory, and the formulae he developed," have been shown to be applicable 
from the point of view of wave mechanics provided that the temperatures 
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Fig. 9. H,O. Curve a represents values of B due to the dispersion effect, curve b those due to 
dispersion plus dipole effect. 


are sufficiently high.'® The result is about 5X 10-* for the quadrupole mom- 
ent. But there appears to be no way at present of checking this value by 
independent experiments. 


16 W.H. Keesom, Phys. Zeits. 22, 129 (1921). 
16H. Margenau, reference 11. 
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ABSTRACT 


In the case of a gravitating mass of perfect fluid which has come to thermo- 
dynamic equilibrium, it has previously been shown that the proper temperature T> 
as measured by a local observer would depend in a definite manner on the gravita- 
tional potential at the point where the measurement is made. In the present article 
the conditions of thermal equilibrium are investigated in the case of a general static 
gravitational field which could correspond to a system containing solid as well as fluid 
parts. Writing the line element for the general static field in the form 

ds? = gijdxidxj + gud? i, j=1,2, 3, 


where the gs; and g44 are independent of the time ¢, it is shown that the dependence 
of proper temperature on position at thermal equilibrium is such as to make the 
quantity Tov gu a constant throughout the system. 


§ I. INTRODUCTION 


i SEVERAL previous articles the principles of relativistic thermody- 
namics have been discussed! and then applied* to determine the con- 
ditions for thermodynamic equilibrium in the presence of gravitational 
fields. In the case of a gravitating mass of perfect fluid which has come to 
thermodynamic equilibrium, it was shown in the course of the work® that 
the proper temperature as measured by a local observer would depend in a 
definite manner on the gravitational potential at the point where the meas- 
urement is made. Thus if we write the line element, for a spherical mass of 
fluid which has come to equilibrium, in the form 


ds? = — e*(dr? + r2d0? + 1? sin? Odp?) + edt? (1) 


where yw and v are functions of the coordinate 7, the dependence of the proper 
temperature 7°) on position was found to be given by the equation 


d log To 1 dv (2) 

dr 2 dr - 

This result is a very interesting one since uniform temperature through- 
out a system which has come to thermodynamic equilibrium has hitherto 
been regarded as an inescapable part of thermodynamic theory. In ac- 
cordance with this equation, however, the proper temperature is found to 


! Tolman, Proc. Nat. Acad. 14, 268 (1928); ibid. 14, 701 (1928); Phys. Rev. 35, 875 (1930); 
ibid. 35, 896 (1930). 

2 Tolman, Proc. Nat. Acad. 14, 348 (1928); ibid. 14, 353 (1928); Phys. Rev. 35, 904(1930). 

3 See reference 2, last article. 
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increase as we move inwards towards the center of the sphere, and this can 
be qualitatively interpreted by ascribing to heat the property of weight and 
regarding the increase in temperature as we move inwards as necessary in 
order to prevent the flow of heat from higher to lower gravitational levels. 
Integrating equation (2) we obtain the result 
y 


log 7) + dae const. 


or 
Toe’!? = ToV gaa = const. (3) 


and have thus obtained a quantity which has a constant value throughout 
the fluid even though the proper temperature itself does not. This result, 
however has so far only been proved for the equilibrium condition of a 
perfect fluid, and the question naturally arises whether the quantity ToV/ ga; 
would also have a constant value in the case of a more complicated system 
containing solid parts. In the present article we shall investigate temperature 
equilibrium in the case of a general static gravitational field by considering 
that the parts of the system whose temperatures are to be compared are in 
thermal contact with a small connecting tube containing radiation. Such a 
tube may be called a radiation thermometer, and by calculating the change 
in pressure as we go from one portion of the tube to another it will be easy 
to show that the radiation has the same value of Tog throughout the 
thermometer. 


§ 2. THE ENERGY-MOMENTUM TENSOR FOR BLACK-Bopy RADIATION 


In order to solve our problem by the method suggested, we shall need an 
expression for the energy-momentum tensor for black-body radiation and 
shall take this as given by the equation 


; dx, dx, 
T’ = (poo + po)— wa am Oe, (4) 
ds ds 
with 
Poo = 3Po (S) 


where poo is the proper macroscopic density of the radiation at the point of 
interest, Po its proper pressure, and the velocities dx,/ds correspond to the 
macroscopic motion of the radiation; i.e., to the motion, in the coordinate 
system which is being used, of an observer who finds on the average no net 
flow of energy in the radiation field. 

Equation (4) is well known in general relativity as being an expression 
for the energy-momentum tensor of a perfect fluid,* and the primary justifi- 
cation for adopting it as applying to black body radiation resides in its 


‘ See for example Eddington, “The Mathematical Theory of Relativity,” Cambridge1923, 
equations (54.81) and (54.82). 
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general applicability in studying the relativistic mechanics of any system 
whose local properties can be specified by the two scalars, proper density 
Poo and pressure po. And such is evidently the case for black-body radiation, 
with the additional simplification poo = 3 po. 

In addition to this primary justification, however, it will also be of 
interest to start with the usual relativistic generalization of the Maxwell- 
Lorentz electromagnetic equations and derive the above expression for the 
energy-momentum tensor by treating black-body radiation as an electro- 
magnetic phenomenon. 

In accordance with this relativistic generalization, the electromagnetic 
field at any point can be specified by a certain antisymmetric tensor F,, whose 
components are directly related in Galilean coordinates to the classical 
components of electric field strength X, Y, Z and magnetic field strength 
a, B, y. And the energy-momentum tensor is given in terms of the F,, by 
the equation® 


To = — PP ye + Lee POP gy, (6) 


Using Galilean coordinates at the point of interest and substituting for 
the components of F,, the values which they then have, it is found that the 
components of 7*” assume the values indicated by the following typical 
examples*® 


TH = — 3(X? — VY? — Z*) — 3(a? — B? — x?) 

TY = — aB — XY 

T4 = — pBZ+ yvY 

T = 3(X? + Y? + 2?) + 2(a? + B* + Y’) 
Thus in Galilean coordinates 7", T*, 7**, 7, 7, and T* become the classi- 
cal components of the stress in the field, 7, 7*4, and T** become the classical 
components of the Poynting vector, and 7“ becomes the classical density of 
energy. 

If now we consider that the electromagnetic field in question corresponds 

to black-body radiation, and take a proper system of Galilean coordinates 


in which the radiation as a whole is at rest, it is evident that we shall have 
on the average 


X?=Y?=Z* and a*® = ~*? = ¥? (8) 

since the average field strength will be independent of direction, 
XY = YZ =ZX =0 and af = py = ya =0 (9) 
since the lack of phase relations between waves will make it equally probable 


that the instantaneous values of these products will have positive or negative 
magnitude, and 


5 See Eddington, reference 4, equation (77.2). 
6 See Eddington, reference 4, equations (77.41-2-3-4). 








1794 R. C. TOLMAN AND PAUL EHRENFEST 


— pPZ+7¥ = —yX + aZ = —a¥Y + BX =0 (10) 


since there will be no net flow of energy in the proper system of coordinates 
that we are now using. Hence, introducing these results into equations 7, 
it is evident that in proper coordinates the averaged energy-momentum 
tensor for black-body radiation will have as its only surviving components 


Tis T2=TR= ~p and T* = poo (11) 
with 
Poo = 3po (12) 


where poo is the proper macroscopic density at the point of interest and the 
three surviving components of the stress are each equal to the pressure po. 

This expression for the energy-momentum tensor of black-body radiation 
holds, however, only for a set of proper coordinates x, y, 2, t, in which there 
is no. net flow of energy, and to pass to a general set of coordinates x,’, x2’, 
x3’, x,’ we shall have to substitute the above values into the general trans- 
formation equation 

Ox’, Ox,’ 


T! Sie —_—- 78 


OXq ONXg 


Doing so we at once obtain 





ge Soars 
Ox Ox Oy Oy Os Os At A 





Ox,’ Ox,’ Ox,’ dx,’ Ox,’ dx, Ox,’ dx Ax,' dx,’ 
Tl = (poo + Po) » po( . ’ (13) 
ot dl 
For the macroscopic velocity of the radiation in our new set of coordin- 
ates, however, we can evidently write 
dx,’ Ox,’ dx Ox,’ dy ax,’ dz _ Ox,’ dt dx,’ (14) 
ds Ox ds dy ds ds ds ot ds ot 


owing to the null value for the components dx/ds, dy/ds and dz/ds and the 
equality of dt and ds in proper coordinates. And from the transformation 
equation 
dx,’ dx,’ 
i at conan: damned 
- = 


OXq OXg 
we obtain 
dx,’ dx = Ax,’ dx! = Ax,’ dx,’ — Ax,' dx,’ 


ge = -— — —+— —. (15) 
Ox Ox dy Oy dz 02 ot at 














Hence, substituting (14) and (15) into (13) and dropping primes we at once 
obtain the result which was to be proved 


dx, dx, 








T” = (poo + Po) — g*" Do 


ds_ ds 
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with 
poo = 3po. 


The second justification which we have thus presented for our expression 
for the energy-momentum tensor for black-body radiation suffers from the 
fact that it makes our final expression, which contains none but macro- 
scopically measurable quantities, a consequence of the relativistic general- 
ization of the Maxwell-Lorentz equations, which assume possibilities of 
microscopic measurement in conflict with the modern ideas of quantum me- 
chanics. Nevertheless, in view of the primary justification which we presented, 
it is evident that our expression for the energy-momentum tensor for black- 
body radiation may be used with reasonable confidence at the present time, 
and may actually be regarded as more certain than the basis of the second 
justification. At some future time, the expression might have to be modified 
in the light of a successful unified field theory if that should ever be achieved. 


§ 3. THE LINE ELEMENT 


Having justified the expression which we shall use for the energy-momen- 
tum tensor of radiation, we are now ready to proceed to our problem of de- 
termining the distribution of radiation in a static gravitational field. 

For the case of a general static system which may contain solid parts we 
can not assume spherical symmetry as a necessary accompaniment of equi- 
librium, and must take the line element in the general static form 


ds? = gidxdxuj t+ gadt? i,j =1,2,3 (16, 


where g;;dx dx; is a negative quadratic form. 

We adopt the convention of using Latin indices 7, j, etc. to correspond 
with the spatial coordinates x;, x2 and x3, and shall reserve Greek indices a, 
8 etc. to correspond with all four coordinates x;, x2, x3 and ¢t. In accordance 
with the usual definition of a static system we take the potentials gis, gos and 
£31 equal to zero, and take the other potentials g;; and gs as independent of 
the time ¢, although depending in any arbitrary way desired on the spatial 
coordinates x, x2 and x3. 

For the potential g4;, we note from the form of the line element that we 
have the simple relation 


iia eel (17) 
£44 


§ 4. ENERGY-MOMENTUM TENSOR FOR BLAcK-Bopy RADIATION 
IN A STATIC FIELD 


We must now consider the form taken by the energy-momentum tensor 
for black-body radiation in the field defined by the above line element. 
Returning to our general Eq. (4) for the energy-momentum tensor 


dx, dx, 


7” = + 
(poo + po) i 








— gos (18) 
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we note that in a static system the macroscopic velocities dx,/ds will evidently 
be zero for u=1, 2,3 


d xX; 


ds 





= 0 (19) 


and taking account of Eq. (17) will reduce for the case uy =4 to 
dx dt 1 
—=—=—— = VE". (20) 
ds ds V £44 
Substituting in (18), the energy-momentum tensor degenerates into 


Tii = go 


(21) 
T*# = gitpoq 
And on lowering suffixes we have 
Tj = 8iaT'® = — giak*Po = — gpa 
Ts = 84:8**p00 = poo 
so that the only surviving components become 
Ti=T.=T;=—po Ti= po. (22) 


§ 5. APPLICATION OF THE PRINCIPLES OF MECHANICS 


We can now investigate the pressure of radiation in our thermometer 
by applying the principles of relativistic mechanics in the form of the well- 
known equation 








= 0. (23) 


Taking the case w= 1 and substituting Eqs. (21) and (22), we obtain 











—(— pv¥— 9 - —(- gipp— 9 — 


1 — Of 45 
— — (gpoov — g) =0 
Ox, 2 Ox, 2 . - : 


0 xy 


which can evidently be rewritten in the form 


___ Apo a/-—g 1 =a 5 OB a “) 








fg + fa — al = Bh +8 
Ox, is Ox, 2 Ox, Ox, 





Of 44 


1 — 
+ > (P00 + pV — gett ; = 0 (24) 


%1 
This equation can easily be simplified, however, since we are permitted to 
write in accordance with a well-known result of tensor analysis 


Og 


15 O85 4 0844 = goo S88 
Ox 


Ox . Ox, Ox, 





g (25) 


1 
g 
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and if this is substituted in Eq. (24), it is found that the second and third 
terms will cancel. Making other obvious simplifications, we then obtain for 
the dependence of pressure on position the simple expression 
Opo 4 poo + po 9 log gas 
Ox, 2 Ox, 





= 0 (26) 


and similar relations will hold with respect to the other spatial coordinates 
x, and Xe. 
For the case of radiation, moreover, we have the additional simplification 


poo = 3po 
so that equation (26) further reduces to 


0 log po 1 Pi log g44 


= 0 (27) 
Ox, Ox, 


and since similar relations hold for the other spatial coordinates, we can 
express the dependence of radiation pressure on position in our gravitational 
field by the remarkably simple equation 


po(gss)? = const. (28) 


§ 6. DEPENDENCE OF TEMPERATURE ON POSITION 


Finally, however, in the case of radiation we can connect proper pressure 
and proper temperature by the well-known result of Boltzmann 


po = —To (29) 


where a is the Stefar Boltzmann constant, and substituting this into Eq. (28) 
we at once obtain the desired conclusion for the dependence of proper tem- 
perature on position in the case of a general static gravitational field 


ToV £44 = const. (30) 
§ 7. Discussion 

In conclusion several remarks may be made with respect to the above 
result which will be of interest. 

In the first place it should be noted from the method of derivation that 
the constancy of To+/gy has been proven in the first instance sole'y for 
points inside the radiation thermometer. Nevertheless since we shall expect 
To and gs; to be continuous functions of position, we shall feel justified in 
concluding that To*/gq is also constant in the system itself where it comes 
in thermal contact with the thermometer. 

In the second place it should be noted that the derivation was carried 
out on the assumption that the system had already been provided with a 
radiation thermometer connecting the parts whose temperatures were to be 
compared. Hence in the case of a given system of interest the question arises 
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whether a thermometer can be inserted to connect the desired points without 
thereby seriously altering the system itself. Thus if we had a gravitating 
system containing solid parts it would be necessary to make a hole into the 
solid and insert a radiation thermometer if we wished to obtain information 
as to the temperature of the interior by the method that we have suggested. 
This procedure, however, would certainly affect the gravitational potentials 
Zu» Which are themselves completely determined by the distribution of the 
matter and energy in the system. Nevertheless since the equation of con- 
nection 

dx, dx, 1 

— S8rpo - = — 847" = GY” — —Gg’ + Age (31) 

ds ds 2 
is a differential one giving the distribution of matter and energy in terms of 
the g ,» and their first and second differential coefficients, it seems correct to 
assume that the insertion of a thermometer of small dimensions can be made 
without appreciably affecting the values of the g,, themselves. This 
question might bear further investigation, however, since singular cases of 
interest might be found. 

Finally, it might be emphasized that although the proper temperature 
itself, 7, varies from point to point in a gravitational system which has 
come to equilibrium, nevertheless the constancy of the combined quantity 
ToV gu provides many of the advantages of the older principle of constant 
temperature throughout as necessary for equilibrium. Indeed it would be 
possible to label To+/ga, as the temperature of a system, except for the 
undesirability of multiplying the different things that are signified by that 
word. In this connection it is also interesting to recall that Einstein himself 
was led in his early speculations on the nature of gravitation’ to distinguish 
between a quantity, called “wahre Temperatur,” which would be constant 
throughout a system in thermal equilibrium and a second quantity, called 
at the suggestion of Ehrenfest “Taschentemperatur,” which would vary 
with gravitational potential. The considerations were of only a limited 
applicability since this was at a time before Einstein’s complete development 
of the general theory of relativity; the quantities in question, however, were 
quite analogous to our present To gy and 7». 


7 Einstein, Ann. d. Physik 38, 443 (1912). 
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LETTERS TO THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may 
be secured by addressing them to this department. Closing dates for this depart- 
ment are, for the first issue of the month, the twenty-eighth of the preceding month; 
for the second issue, the thirteenth of the month. The Board of Editors does not 
hold itself responsible for the opinions expressed by the correspondents. 


The Inner Potential of a Copper Crystal 


From electron diffraction experiments planes is one-half the above value. With a 


Rupp (Ann. d. Physik 5, 453, 1930) obtains 
values for the inner potentials of several 
metals. His value of the inner potential for 
any one metal is approximately constant and 
independent of the electron wave-length. In 
some results for a copper crystal obtained by 
the writer (Phys. Rev. 34, 679, 1929) it ap- 
pears that the inner potential decreases with 
increase of wave-length or decrease of voltage. 
It is to be noted that, in his calculation for 
copper, Rupp uses the distance between ad- 
jacent (100) planes as 3.59A. This is the value 
for the edge-length of an elementary cube, 
but since the copper crystal is a face-centered 
cube, the distance between adjacent (100) 


gas-free Cu crystal Rupp obtains maxima 
at 54 and 134 volts for which he chooses the 
values 3 and 4, respectively, for n. If the 
distance 1.8A between adjacent planes is 
used, the above values of m must be divided 
by two to give the same values for the inner 
potential. This, however, requires the use of 
a half-integer. If only whole integers are 
used then equal values for the inner potential 
from the observed two maxima are not ob- 
tained. 
H. E. FARNSWORTH 
3rown University 
Providence, R. I. 
November 15, 1930 


Differences in the Absorption Spectrum of Benzene in the Liquid and Vapor State 


The absorption spectra of benzene, liquid 
and vapor, in the 1.14 and 1.6u regions have 
been recorded with an effective slit width of 
6A. Four meters of saturated vapor and two 
millimeters of liquid were used, both at room 
temperature and pressure. Several interesting 
differences between the liquid and vapor ab- 
sorptions appear when this resolution is at- 
tained. The vapor spectrum is shifted towards 
higher frequencies. The 1.14 band is shifted 
28 cm~!. Two shifts appear in the 1.6z 
region, 18 cm~! and 28 cm™, due to the fact 
that the 1.64 absorption is a combination of 
overtones of two fundamentals. The ab- 
sorption in the 1.6% region is similar in 
intensity for the two phases while the vapor 
absorbs only 4—6 percent as much as the liquid 
in the 1.1 region. 

These preliminary results yield information 


as to the distribution of energy in the vibra* 
tional states of the molecule in the two phases 
and show the effect of the proximity of the 
molecules in the liquid state. The two shifts 
in the 1.64 region are about equally divided 
This 
enables a separation of the overtones of the 
two fundamentals that overlap in this region. 
An extensive investigation of this general 


between the twenty lines observed. 


problem is now under way and a compre- 
hensive report will be published in the near 
future. 
E. D. McALISTER 
H. J. UNGER 
Division of Radiation & Organisms 
Smithsonian Institution 
Physics Department. 
University of Oregon 
November 25, 1930. 
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Relative Intensities in Hyperfine Structure Multiplets 


The extension of the Landé interval rule 
to hyperfine structure energy levels was first 
made by Bach and Goudsmit! from observa- 
tions of hyperfine patterns in bismuth. Since 
the observed intervals followed very closely 
the interval rule, the nuclear spin, J, was 
calculated to be (9,2)(h/27). This was later 
verified by them? in the Zeeman splitting 
up of each line into ten components, in a 
paper which is now considered a classic. 

The interval rule has since been shown by 
the author to hold for over a hundred hyper- 
fine patterns in praeseodymium. With the 
partially resolved patterns of three lines in 
manganese, the patterns or graphs of many 
lines were predicted for Mn I by making the 
assumption that the intensity rules as well 
as the interval rule were, in general, valid. 
These predictions fit the observed patterns 
so beautifully that the interval rule was used 
in computing all of the hyperfine levels.* 
The interval rule, as pointed out in a short 
note by Goudsmit and Fisher,* should hold 
in general and is therefore a powerful tool 
that may be used in analysing hyperfine- 
structure patterns. 

Since each tiny hyperfine multiplet in Mn I 
has not been completely resolved the intensi- 
ties of several lines were qualitatively given 
by the authors’ in diagrams or graphs. With 
intensity rules® similar to those for ordinary 


1 Bach and Goudsmit, Zeits. f. Physik 43, 
321 (1927). 

2 Bach and Goudsmit, Zeits. f. Physik 47, 
174 (1928). 

3 White and Ritschl, Phys. Rev. 35, 1146 
(1930). 

4 Goudsmit and Fisher, Paper presented at 
the Amer. Phys. Soc. November 29, 1930. 

’ For hyperfine intensity rules see Hill, 
Nat. Acad. Sci. 15, 779 (1929). 


relative intensities have been 
calculated for quite a number of the hyper- 
multiplets already published. The 
qualitative agreement with the observed 
patterns is so remarkable for all patterns 
studied that accurate measurements are now 
contemplated. A comparison of the two 
figures, here shown, with the graphs and 
microphotometer curves of the same lines® 
It is 
now clear why the off-diagonal or faint com- 
ponents are rarely observed. 


multiplets, 


fine 


strongly confirms their interpretation. 
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While very few accurate intensity measure- 
ments are known for hyperfine lines, quali- 
tative measures indicate that in general the 
intensity rules hold as well as the interval 
rule. This should be true even though the 
ordinary multiplet to which the hyperfine 
pattern belongs does not follow LS coupling 
or the intensity rules. 

, H. E. Wuite 

University of California, 

Berkeley, California 
November 22, 1930. 


The Vibrational and Rotational Analysis of the S, Bands. 


The absorption and emission spectrum of 
S; has been photographed from 45300 to \2400 
withgreater dispersion than that used byother 
investigators. The absorption spectrum was 
obtained by means of a quartz cell which con- 
tained pure sulphur and could be heated elec- 
trically up to about 700°C. A hydrogen con- 
tinuous source was used. The absorption pho- 
tographs were taken with an E, Hilger quartz 





spectrograph, and also in the first order of a 
21-foot Rowland grating. For the rotational 
analysis we have used emission plates taken 
in the third order of the Rowland grating. 
Our work confirms essentially the vibra- 
tional analysis as given by Rosen (Zeits. f. 
Physik 43, 69, 1927). However, the assign- 
ment of quantum numbers to some of the 
bands has been changed and several other 
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band heads have been measured. Numerous 
perturbations have been found mainly among 
the levels of the upper electronic state. v’ 
= 10 is the first level showing diffuseness due 
to predissociation (c.f. V. Henri, Structures 
des Molecules, Paris, 1925). This is shown 
definitely in both our emission and absorption 
plates. In absorption all bands with v’ =10 
are diffuse, while in emission no bands are 
obtained with v’=10. The 9-0 band does not 
appear in our emission plates. This can be 
attributed to the small probability in the 
Franck-Condon parabola rather than to pre- 
From v’=10 to v’=15 the dif- 
fuseness decreases and for v’=16 it increases 
greatly. Our interpretation of the predisso- 
ciation differs somewhat from that of Herz- 
berg (Zeits. f. Physik 61, 604, 1930). 

We have used for the fine structure analysis 
the bands 9-1, 7-0, 8-1, 9-2, and 7-1. These 
were photographed in emission in the third 
order of a 21-foot Rowland grating with a dis- 
persion of 0.8A per mm. The overlapping is 
such that even with this dispersion a consider- 
able number of lines are not clearly separated. 
The lines can be arranged in two P and two 
R branches, one series of each kind being con- 
siderably stronger than the other. Each of 
the stronger branches, however, split into two 
series of equal intensity for high quantum 
numbers. This structure (three P and three 
R branches of approximately equal intensity) 
is that of a *=—>*Y transition, and the bands 
are presumably analogous to the Schumann- 
Runge bands of O, (82,~—*2,,~). The rota- 
tional constants B,’ (extrapolated) and B,” are 
0.319 and 0.414 cm™. These give r,.’=0.573 
and r.”=0.503A. The upper rotational states 
belonging to v’=9 are irregular in that their 
spacing corresponds to a larger value of B 
than that for the lower rotational states. This 
is thought to be due to perturbations in the 
rotational states, since these levels are close 
Just as in the 
case of Os, alternate rotational levels, hence 


dissociation. 


to the predissociation limit. 


1801 


alternate band lines, are missing in each se- 
ries, showing that the sulphur atom S* has 
zero nuclear spin. In the lower *> state, the 
levels with even K values are missing, in the 
upper state the odd K values, as in the Schu- 
mann-Runge bands of Q). 

Our analysis does not agree with that given 
recently by Rompe (Zeits. f. Physik 65, 404, 
1930), according to which the bands are due 
to a 'II-—>' transition. Rompe’s analysis was 
based on flourescence spectra of sulphur ex- 
cited by various lines of the mercury arc and 
were photographed with relatively low disper- 
sion. He did not use single line excitation, and 
although flourescence spectra in general facili- 
tates the analysis of complex bands, the use 
of several exciting lines simultaneously, the 
proximity of the sulphur levels as well as the 
dispersion used makes it difficult to interpret 
the results. 

The analysis is complete except for a num- 
ber of weak satellite heads found mainly in 
absorption which are probably due to a differ- 
ent electronic transition. 

Gilles (C.R. 188, 1607, 1929) and Van Idde- 
kinge (Nature 51, 858, 1930) have reported a 
new system due to S; in the ultra-violet. These 
bands, however, are due to CS and not to 
Se, and the vibrational analysis of this system 
has been given by Jevons (Proc. Roy. Soc. 
A117, 351, 1927-1928). They are very in- 
tense and a small impurity of carbon in the 
discharge tube is sufficient to bring them out 
prominently. 

A detailed discussion of both the rotational 
and vibrational analyses will appear shortly. 
The authors wish to express their apprecia- 
tion to Professor Mulliken for his interest and 
assistance. 

S. MEIRING Naup£t 
ANDREW CHRISTY 


Ryerson Physical Laboratory, 
University of Chicago, 
December 12, 1930. 


Erratum 


I note in my letter concerning efficiency of 
emission of oxide-coated filaments published 
in the October 15th issue of the Physical Re- 
view that the value of A has been printed as 
0.3 ampere per cm? degree*. This is a typo- 


graphical error and should be 0.03 ampere per 
cm? degree’. B. J. THomMPpson 
Vacuum Tube Engineering Department 
General Electric Company, 
November 5, 1930. 
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Electrodeless discharge, current, pressure, and 
frequency relations—M. L. Braun—1195 

Electrodeless, Hg vapor, characteristics, H. Smith, 
W. A. Lynch, N. Hilberry—374(A) 

Geissler, hollow cathode, characteristics, C. H. 
Thomas—374(A) 

Heats of condensation of electrons on Mo, Pt 
and W-K, C. C. Van Voorhis, K. T. Compton 
1435 

Hollow cathode, excitation processes, spectra, R. 
A. Sawyer—44 

Luminosity, propagation of, J. W. Beams 


Farns- 


Bennett 


997 


ANALYTIC Sl 


BJECT INDEX 


Ne arc, metastable atom and resonance radiation, 
I. Langmuir, C. G. Found—604(L) 
New type of glow-discharge tube, H. J. 
373(A) 
Organic reactions in gaseous electrical discharge, 
EK. G. Linder—374(A), 1375 
Positive ion sheaths, effect of collisions on poten- 
tial distribution, H. P. Robertson, P. M. Morse 
375(A) 
Secondary electrons produced by metastable atom 
and positive ions, W. Uyterhoeven, M. C. Har- 
rington —709 


Reich 


Spectra, variation with discharge, characteristics, 
E. O. Hulburt—13 


Vapor jets in Hg arc, E. Kobel—1636 
Dissociation 
Nitric acid from absorption spectra, E. L. Kinsey, 
J. W. Ellis—603(L) 


Hy, velocities of H* ions formed, W. W. Lozier— 
1285 H,0 and formation of negative ions in H, 
and H.O, W. W. Lozier—1417(L) 
Dissociation, heat of (see also Molecular structure 
and constants) 


Oo, heat of formation, L. C. Copeland—1221 
Elasticity 
Longitudinal vibration of rods, D. S. Muzzey, Jr. 
935 
Plastic deformations of ductile metals, A. Nadai 
762 
Electrets 
X-ray study of structure, M. Ewing—378(A) 


Electrical conductivity and resistance 
Contact between carbon spheres, microphonic 


action, F. Gray—375(A); F. S. Goucher,—375 
(A) 

Contact resistance, wave-mechanical theory, J. 
Frenkel—1604 


Electrical discharge in gases (see Discharge of elec- 
tricity in gases) 
Electromagnetic theory 
Diffraction of a circularly symmetrical electro- 
magnetic wave by a co-axial circular disc, J. 
Bardeen—1482 
Transmission of alternating currents through a 
semi-infinite heterogeneous conducting medium, 
H. P. Evans—1579 
Unified theory, spherically symmetrical field, N. 
Rosen, M. S. Vallarta—110 
Electrons 
Deflection in magnetic field on wave mechanics 
and value of e/m, L. Page—444; criticism, C. 
Eckart—1014(L), 1514(L); reply, L. Page 
1418(L); G. E. Uhlenbeck, L. A. Young—1721; 
M. S. Plesset-—1728; E. H. Kennard—1667(L) 
e, by x-ray refraction, H. E. Stauss—1101 
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e/m from velocity measurements, C. T. Perry, 
E. L. Chaffee—904 
e/m measurements from inertia effects, S. J. Bar- 
nett—786(A) 
Heat of condensation of electrons on Mo, Pt and 
W-K, C. C. Van Voorhis, K. T. Compton—1435 
Polarization by double scattering, C. T. Chase— 
1060 
Recombination with Cs*, E. H. Kurth—374(A) 
Recombination with protons, theory, E. C. G. 
Stueckelberg, P. M. Morse—16 
Reflection coefficients and contact potentials, H. 
E. Farnsworth, V. H. Goerke—1190 
Scattering by gas molecules, M. C. Green—239, 
786(A) 
Scattering and operation of Geiger counter, C. T. 
Chase—984 
Electrons, diffraction 
Inner potential of a copper crystal, H. E. Farns- 
worth—1799(L) 
Electrons in gases 
Formation of, W. Bleakney—1303 
Ionization of helium, neon, and argon by electron 
impact, efficiency, P. T. Smith—1293 
Scattering in He, A and Hz2, angular and energy 
distribution, J. H. McMillen—1034 
Electrons, photoelectric (see Photoelectric effect) 
Electrons, scattering (see also Electrons in gases) 
Effect of small angle scattering on absorption 
coefficient, M. C. Green—239, 786(A) 
Electrons, secondary 
From metals bombarded by positive ions, C. L. 
Utterback, W. Geer—785(A) 
Metastable atoms and positive ions, produced by, 
W. Uyterhoeven, M. C. Harrington—709 
Velocity distribution of secondary electrons from 
molybdenum, T. Soller—1212 
Electron tubes 
Cathode ray oscillograph, F. Bedell, J. Kuhn— 
993 
High potential x-ray tube, C. C. Lauritsen, B. 
Cassen—988 
Energy states of atoms 
Average life of ionized helium, L. R. Maxwell— 
379(A) 
Method of calculation, C. Eckart—149(L) 
Energy states of molecules 
As effected by interaction, E. C. Kemble, F. F. 
Rieke—153(L) 
Excitation function 
He, spectra, W. C. Michels—604(L), 1362 
Filters, light 
Infrared, controllable transmission, A. H. Pfund 
—71 
Fine structure (see also Hyperfine structure) 
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He 2°P level, test of spin interaction, G. Breit— 
383 
K radiation of lighter elements, L. Y. Faust—161; 
M. Soderman—1414(L) 
K series, Ca to Ga, B. Kievit, G. A. Lindsay—648 
Ka of Fe, Co, Ni, Cu, Moand Ag, KB of Mo, J. 
Valasek—1523 
MoKa lines, double crystal spectrometer, N. S. 
Gingrich—1050 
Unmodified scattered x-rays, J. A. Bearden—791 
Fluorescence 
Polarized, of rhodamine-B and uranine, E. Mer- 
ritt, D. R. Morey—1386 
Quartz, under low voltage cathode rays, H. Peters 
—1631 
ZnS, under a, 8, y rays, decay characteristics, G. 
S. Gessner—207 
Gamma rays 
Wave-length in spectrum of Ra, RaB, RaC, and 
RaD, including intensity measurements of line 
and continuous spectra, L. T. Steadman—460 
Geiger-counters 
Sensitivity, C. T. Chase—984 
Geophysics 
Wind mixing and diffusion in the upper atmos- 
phere, S. Chapman—1014(L); reply, E. O. 
Hulburt—1264(L) 
Hall effect 
Te amalgams, P. I. Wold, J. M. Hyatt—375(A) 
Heat of dissociation (see Dissociation) 
High voltages 
Production, M. A. Tuve—1576; M. A. Tuve, L. 
R. Hafstad and O. Dahl—1261(L) 
Hyperfine structure 
Cd lines and nuclear spin, C. L. Albright—847 
Hyperfine structure and interval rule, H. E. White 
—1800(L) . 
S and P terms of Li*, G. Breit, F. W. Doer- 
mann—1732 
W Ka, F. K. Richtmyer, S. W. Barnes, K. V. 
Manning—1017(L) 
Infrared spectra :see also Spectra, molecular; 
Hg, arc, H. J. Unger—784(A) 
High resolution in near infrared, E. D. McAlister— 
784(A) 
Instruments (see Methods and instruments) 
Intensities in spectra 
Bands of diatomic molecules, theory, E. Hutchis- 
son—410 
Cameron bands, intensity changes, H. P. Knauss, 
J. C. Cotton—1099(L) 
Ni I, Ni II, Co I multiplets, L. S. Ornstein, T. 
Bouma—679 
X-ray lines in the tantalum L-spectrum, V. Hicks 
—785(A), 1273 
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X-ray satellites, intensity ratios, F. K. Richtmyer, 
L.S. Taylor—1044 
Ionization by impact 
CO; by electrons, H. D. Smyth, E. C. G. Stueckel- 
berg—472 
He, Ne, and A, efficiency by electrons 0-4500 
volts, P. T. Smith—1293 
K and Hg vapors, extra ionization potentials, A. 
L. Hughes, C. M. Van Atta—214 
Multiply charged ions, He, Ne and A, W. Bleak- 
ney —1303 
NO, and N,O, by electrons, E. C. G. Stueckelberg, 
H. D. Smyth—478 
Velocities of H* ions formed, W. W. Lozier—1285 
Ions 
Cs* in He, pressure changes, A. C. G. Mitchell 
381(A) 
Negative ions in hydrogen and water vapor, W. 
W. Lozier—1417(L) 
Recombination of H* with electron, theory, E. C. 
G. Stueckelberg, P. M. Morse—16 
Ions in gases 
Energy losses at elastic collisions, A. M. Cravath 
—248 
Recombination in A, Nz and Ho, O. Luhr—24, 
788(A); in air and O,—787(A) 
Ions, mobility 
Aged ions in air, O. Luhr, N. E. Bradbury—1394 
In dry and moist air, J. Zeleny—35 
Nat in H2, L. B. Loeb—152(L), 790(A) 
Isotopes 
Be’, band spectrum of BeH, W. W. Watson— 
1019(L) 
Cs, K. T. Bainbridge—1668(L) 
Effect on band spectrum intensities, J. L. Dun- 
ham—1553 
Nis, O1s, O17, abundances, from band spectra, S. 
M. Naudé—333 
Number of isotopes of eight metals, F. Allison, 
E. J. Murphy—1097(L) 
Kinetic theory of gases 
Second virial coefficient, experimental and theore- 
tical, H. Margenau—1782 
Light waves 
Electro-optical modification of, electric oscillations 
on Kerr cell, L. H. Stauffer—1352 
Frequency shifts in dispersing media, G. Breit, E. 
O. Salant—871 
Luminescence 
Due to a, 8, y rays, mechanism, D. H. Kabakjian 
—379(A) 
ZnS under a, 8, y rays, decay characteristics, 
G. S. Gessner—207 
Magnetic properties 
Barkhausen effect in Ni, discontinuous changes in 
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length, C. W. Heaps, A. B. Bryan—326, criti- 
cism, L. W. McKeehan—1014(L) 

Cu and Ag, susceptibilities of and structure sensi- 
tivity, F. Bitter—978 

Electrical resistance. tension and magnetization of 
Ni and permalloy, bibliography included, L. 
W. McKeehan—948 

Gyromagnetic anomaly and rotation by mag- 
netization, S. J. Barnett—789(A) 

Hall effect in ferromagnetic materials, E. M. Pugh 
—1503 

Hall effect of Te amalgams, P. I. Wold, J. M. 
Hyatt—375(A) 

Isotropy of a paramagnetic alum, C. G. Mont- 
gomery—1661 

Liz nucleus, momentum of, L. P. Granath—1018 
(L); see also G. Breit, F. W. Doermann—1262 
(L) 

Molecular field and atomic order in ferromagnetic 
crystals, L. W. McKeehan—1512(L), 1670(L) 
Pt-Co and Pd-Co alloys, F. W. Constant— 

-786(A), 1654 
Susceptibility measurements, S. Freed, C. Kasper 
-1002 
Susceptibilities of Bi single crystal, A. B. Focke 
—319 
Susceptibility of Cu crystals, C.G. Montgomery— 
498 
Susceptibility of gases, temperature dependence, 
F. Bitter—1648 
Temperature change on magnetization, W. B. 
Ellwood—1066 
Magnetostriction oscillator 
For elastic studies, D. S. Muzzey, Jr.—935 
Mechanics 
Coupled systems of linear oscillators, A. B. Lewis 
—568 
Solar system, quantization and inverse square law 
A. E. Caswell—787(A) 
Vibrations of mechanical systems representing 
organic molecules, R. C. Yates—555 
Vibrations of non-planar membrane, G. R. Sti- 
bitz—513 
Vibration of system representing benzene ring, 
R. C. Yates—563 
Mechanics, quantum 
Complex spectra, theory, E. U. Condon—1121 
Calculation of energy values, C. Eckart—149(L) 
Deflection of electrons in magnetic field and value 
of e/m, L. Page—444; criticism, C. Eckart— 
1014(L); 1514(L); reply, L. Page—1418(L); G. 
E. Uhlenbeck, L. A. Young—1i721; M. S. 
Plesset-—1728; E. H. Kennard—1667(L) 
Dielectric constant Hg, J. V. Atanasoff-——1232 
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Electrical resistance of contacts, J. Frenkel—1604 
Energy states of molecules as affected by inter- 
action, E. C. Kemble, F. F. Rieke—153(L) 
Extrapolation of atomic structure factor curves, 
D. K. Froman—1339 
Frequency shifts in dispersing media, G. Breit, 
E. O. Salant—871 
Intensities in band spectra, E. Hutchisson—410 
Isotope effect on intensities, J. L. Dunham—1553 
Magnetic moment of Li; nucleus, G. Breit, F. W. 
Doermann—1262(L) 
Motion of nucleus, effect on spectra, D. S. Hughes, 
C. Eckart—694 
Orbital valency, J. H. Bartlett, Jr.—1096(L) 
Paschen-Back effect, incomplete, J. B. Green— 
157, 379(A) 
Principle of identity and the exclusion of quantum 
states, G. N. Lewis—1144 
Recombination of electron and Ht, E. C. G. 
Stueckelberg, P. M. Morse—16 
Relativity, unified field, W. Band—1405 
Rotational motion of molecules in crystals, L. 
Pauling—430 
Screening constants, calculation of, C. Eckart — 
878 
Spin interactions as tested by fine structure of He, 
G. Breit —383 
Structure of emission line, theory, F. C. Hoyt— 
380(A), 860 
Symmetrical top, dispersion formula, Raman ef- 
fect—an erratum, M. Muskat—363 
Wave equation, solution, L. A. Young, G. E. 
Uhlenbeck—1154 
Wave functions, approximate, for any atom, J. C. 
Slater—57 
Wave functions for Be, B, C, N, O, F, Ne, C. 
Zener—51 
Mechanics, statistical 
Brownian motion, theory of, G. E. Uhlenbeck, L. 
S. Ornstein—823 
Entropy of H2, D. MacGillavry—1398, W. F. 
Giauque, H. L. Johnston—1592(L) 
Fermi-Thomas method applied to positive ions, 
E. B. Baker—630 
Metastable atoms 
Ne, in arc by resonance radiation, I. Langmuir, C. 
G. Found—604(L) 
Meteorology (see also Atmospheric electricity) 
Cosmic rays and meteorology, R. A. Millikan— 
1595 
Methods and instruments 
Calorimeter to measure small heat changes, W. 
B. Ellwood—1066 
Cathode ray oscillograph, correction, F. Bedell, 
J. Kuhn—993 


Diamagnetic susceptibility measurements, S. 
Freed, C. Kasper—1002 

False x-ray lines, J. M. Cork—665 

Filters, infrared, controllable transmission, A. H. 
Pfund—71 

Fluorescent dry plates for photographic photo- 
metry, G. R. Harrison, P. A. Leighton—779(L) 

Geiger counter, voltage effects, C. T. Chase—984 

Glow discharge tube—new type, H. J. Reich— 
373(A) 

High potential x-ray tubes, C. C. Lauritsen, B. 
Cassen—988 

High voltage tubes, M. A. Tuve, L. R. Hafstad 
and O. Dahl—1261(L) 

Infrared, high resolution in, E. D. McAlister— 
784(A) 

Low grid-current vacuum tube for current meas- 
urements, G. F. Metcalf, B. J. Thompson—1489 

Mass spectrograph, H. D. Smyth, E. C. G. Stuec- 
kelberg—472 

Mechanical models to represent dynamic proper- 
ties of molecules, C. F. Kettering, L. W. Shutts, 
D. H. Andrews—531 

Minimizing errors in reading distance on photo- 
graphic plates, D. Cooksey, C. D. Cooksey —80 

Optically excited radiation, P. Bender—1535 

Photoelectric cells, methods of enhancing sensi- 
tiveness, O. R. Olpin—251 

Photographic method of determining atomic struc- 
ture factors, D. K. Froman—1330 

Pierce acoustic interferometer, high frequencies, 
E. Klein, W. D. Hershberger—1262(L) 

Rectifier, cupric sulfide-magnesium junction, M. 
Bergstein, J. F. Rinke,C. M.Gutheil—587; 1022 

Refraction of x-rays, H. E. Stauss—1101 

Shot-circuit impedence measurements, H. N. 
Kozanowski, N. H. Williams—1314 

Specific heat variation with pressure, E. J. Work- 
man—1083 

Velocity Measurement of electrons, C. T. Perry, 
E. L. Chaffee—904 

Viscosimeter, Margules method for absolute values, 
H. R. Lille—347 

X-ray spectrometer, double crystal, J. W. M. 
Dumond, A. Hoyt—1702 


Microphone 


Microphonic action of contact resistance, F. Gray 
—375(A); F. S. Goucher—375(A) 


Mobility of ions (see Ions, mobility) 
Molecular structure and constants (see also Spectra, 


molecular) 

Atomic J values and molecular quantum numbers, 
R. S. Mulliken—1440 

Benzene ring model and Raman spectra, R. C. 
Yates—563 
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Homopolar chemical bond and Raman spectra, 
R. C. Yates—555 

K2, heat of dissociation, vibrational constants, 
W.O. Crane, A. Christy—421 

Mechanical models to represent dynamic proper- 
ties, C. F. Kettering, L. W. Shutts, D. H. 
Andrews—531 

Methane gas, new absorption band, D. M. Den- 
nison, S. B. Ingram—1451 

Oz, heat of formation, L. C. Copeland—1221 

Orbital valency, J. H. Bartlett, Jr.—-1096(L) 

Multiplets 
Intensity ratios Ni I, Ni II, Co I, L. S. Ornstein, 
T. Bouma—679 
Packing effect 
Cs, K. T. Bainbridge—1668(L) 
Paschen-Back effect 

For intermediate fields, theory, J. B. Green—157, 

379(A) 
Photoelectric cells 

Methods of enhancing sensitiveness, O. R. Olpin 

—251 
Photoelectric effect 

Cs thin films on O, L. R. Koller—1639 

Direction of emission of photoelectrons, E. O. 
Lawrence, M. A. Chaffee—1099(L) 

Effect of intense electric fields, E. O. Lawrence, 
L. B. Linford—482; see also L. B. Linford— 
1100(L) 

Inhibition by infrared light, A. R.#Ipin—376(A) 

Outgassing, action of photoelectrons, R. P. Winch 
—601(L) 

Photoionization of salt vapors, A. Terenin— 
147(L) 

Threshold and current variation with temperature 

(Hg), D. Roller—738 
Photography 

Fluorescent dry plates for photographic photo- 
metry, G. R. Harrison, P. A. Leighton—779(L) 

Unreliability of photographic emulsions for record- 
ing distances, improvements, D. Cooksey, C. D. 
Cooksey—80 

Photons 

Collisions of, A. L. Hughes, G. E. M. Jauncey— 

773; S. Vavilov—1590(L) 
Planck’s constant, / 

Method of computing not involving e, A. Lucy— 
367(L) 

Polarization, electrical (see Dielectric properties) 
Potentials, critical (see also Spectra) 

Ag I, ionization, spectra, H. A. Blair—1531 

COsz, ionization by electrons, H. D. Smyth, E. C. 
G. Stueckelberg—472 

CO II, ionization, spectra, J. H. Findlay—5 

Ionization potentials, successive, calculated, E. B. 
Baker—630 


Ionization of He, Ne and A, W. Bleakney—1303 
K and Hg vapor, extra ionization potentials, A. L. 
Hughes, C. M. Van Atta—214 
NO: and N,O, ionization by electrons, E. C. G. 
Stueckelberg, H. D. Smyth—478 
Proceedings of the American Physical Society 
Ithaca Meeting, June 19 to 21, 1930—372 
Eugene Meeting, June 19 to 21, 1930—783 
Quadrupole radiation 
Experimental evidence, R. Frerichs, J. S. Camp- 
bell—151(L); 1460 
Quantum mechanics (see Mechanics, quantum) 
Radiation 
Penetration, therapeutic use, E. Q. Adams— 
1020(L) 
Quenching of resonance radiation (Hg), M. W. 
Zemansky—919 
Radio 
Diffraction of a circularly symmetrical electro- 
magnetic wave by a co-axial circular disc, J. 
Bardeen—1482 
Radioactivity 
Range of @ particles from Th, G. H. Henderson, 
J. L. Nickerson—1344 
y-ray scattering, ThC, C. Y. Chao—1519 
Wave-length measurements of gamma rays for 
Ra, RaB, RaC, and RaD, L. T. Steadman—460 
Raman effect 
Benzene and diphenyl, anti-Stokes lines, R. W. 
Wood—1431 
Excitation efficiency, relative, of Hg lines, M. 
Werth—1096(L) 
Geometric isomers, C. F. Ffolliott —367(L) 
HBrand NI, E. O. Salant, A. Sandow—1591(L) 
Hg and He arc excitation, filters and improved 
technique, R. W. Wood—1421 
Solutions of NaNOs, V. Sterling, E. R. Laird— 
148(L) 
Structure of organic molecules, D. H. Andrews— 
544 
Trimethylethylene, D. Franklin, E. R. Laird— 
147(L) 
Recombination 
Electron and H*, theory, E. C. G. Stueckelberg, 
P. M. Morse—16 
Of ions in A, N2 and H2, O. Luhr—24, 788(A) 
Electrons with Cs*, E. H. Kurth—374(A) 
Rectifiers 
Cupric sulfide-magnesium junction, M. Bergstein, 
J. F. Rinke, C. M. Gutheil—587; errataa—1022 
Relativity 
Mechanical treatment of temperature distribution 
in case of radiation W. Anderson—365(L.); re- 
ply, R. C. Tolman—365(L) 
Non-statical solution of Einstein’s law of gra vita- 
tion, H. P. Soh—1515(L) 
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Projective relativity,O. Veblen, B. Hoffmann—810 
Relativistic wave mechanics, of deflected elec- 
trons, M.S. Plesset—1728 
Temperature equilibrium in a static gravitational 
field, R. C. Tolman, P. Ehrenfest-—1791 
Unified physical field, W. Band—1405 
Unified theory, spherically symmetrical field, N. 
Rosen, M.S. Vallarta—110 
Resistance of metals (see Electrical conductivity and 
resistance) 
Tension coefficient of, H. Rolnick—506 
Resonance radiation 
Width of line as function of geometry of tube and 
vapor density, M. W. Zemansky—219, 380(A), 
1591(L) 
Screening constants 
Theory and calculation of, C. Eckart—878 
Shielding constants 
For any atom, approximate, J. C. Slater—57 
Shot effect 
Abnormal effect of ions, J. S. Donal, Jr.—1172 
Oxide coated filaments, H. N. Kozanowski, N. H. 
Williams—1314 
Skylight 
Energy distribution in ultraviolet, B. O’Brien 
381(A) 
Spark 
Time lag, relation to photoelectrons from cathode, 
J. A. Tiedeman—376(A) 
Specific heats 
Gases, variation with pressure, calculated values, 
L. G. Hoxton—1091 
Gases, variation with pressure, experimental val- 
ues, E. J. Workman—1083 
Organic crystals, calculated from 
models, A. B. Lewis—568 
Spectra, absorption 
Benzene in the liquid and vapor state, E. D. 
McAlister, H. J. Unger—1799(L) 
Nitric acid dissociation, by absorption spectra, 
E. L. Kinsey, J. W. Ellis—603(L) 
Spectra, atomic 
Ag Il and Pd II, extension, H. A. Blair—173 
Ag I, classification and new lines, H.A. Blair—i531 
Atomic J values and molecular quantum numbers, 
R. S. Mulliken—1440 
Auroral green line, Zeeman effect, quadrupole 
radiation, R. Frerichs, J. S. Campbell—151(L), 
1460 
Co II, classification, J. H. Findlay—S5 
Cu arc, correction to published list, A. G. Shen- 
stone—603(L) 
Diaphantine equation, hydrogen spectrum, T. H. 
Gronwall—1671(L) 
Energy states, calculation, C. Eckart—149(L) 


mechanical 


1817 


Excitation function of helium, W. C. Michels— 
604 (L); 1362 

Excited in hollow cathode discharge, R. A. Saw- 
yer—44 

Fe lines, classification, H. D. Babcock—784(A) 

Forbidden lines and quadrupole terms, L. D. 
Huff, W. V. Houston—842 

He, new lines in arc and spark, P. G. Kruger—855 

Hg lines, fluorescence, polarization of, A. C. G. 
Mitchell—1589(L) 

Hg 2537, ratio of emission to absorption line width, 
M. W. Zemansky—219, 380(A), 1591(L) 

Hyperfine structure Cd lines, C. L. Albright—847 

Hyperfine structure of S and P terms of two elec- 
tron atoms, G. Breit, F. W. Doermann—1732 

Kr and Xe, wave-lengths, C. J. Humphreys— 
380(A) 

Li I and Li II, calculated effect of motion of 
nucleus, D. S. Hughes, C. Eckart—694 

Nitrogen atomic, resonance series, J. J. Hopfield— 
789(A) 

Ni I, Ni II, Co I, intensity measurements, multi- 
plet intensity ratios, L. S. Ornstein, T. Bouma 


—679 

O, singlet system, origin of auroral line, R. Fre- 
richs—398 

Oxygen, atomic, resonance series, J. J. Hopfield 
—789(A) 


Pb III and Pb IV, classification, S. Smith—1 

Pd I analysis, A. G. Shenstone—669 

Polarization of Hg resonance line, S. Mrozowski 
1591(L) 

Quenching of Cd lines by gases, P. Bender—1535 

Rh, term values, A. G. Shenstone, J. J. Livingood 
—380(A) 

Structure of emission line, theory, F. C. Hoyt— 
380(A), 860 

Theory of complex spectra, E. U. Condon—1121 

Wave equation, application to hydrogenic atoms, 
L. A. Young, G. E. Uhlenbeck—1154 

Wave-length list for extreme ultraviolet, J. M. 
MaclInnes, J. C. Boyce—368(L) 

Spectra, continuous 

From gases lighted by strong discharge, E. O. 

Hulburt—13 
Spectra, general 

Excitation function of helium, W. C. Michels— 
604(4); 1362 

Fluorescence of quartz, H. Peters—1631 

Fluorescence, polarization of, E. Merritt, D. R. 
Morey—1386 

From gases lighted by strong discharge, E. O. 
Hulburt—13 

General Electric type S-1, spectrum, D. Dooley— 
1476 
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He, continuous spectrum in condensed discharge 
500-1100A, J. J. Hopfield —784(A) 
Hg vapor absorption bands, origin of, S. Mrozow- 
ski—1168 
Optical excitation, P. Bender—-1543 
Oscillating arc in spectrochemical analysis, E. Z. 
Stowell, W.S. Huxford—1348 
Oxygen, auroral and nebular lines, I. S. Bowen 
600(L) 
Quenching of mercury resonance radiation by 
nitrogen and carbon monoxide, J. Kaplan-—-788(A) 
Radiation from a high potential x-ray tube, C. C. 
Lauritsen—1680 
Tungsten mercury arc, peculiarities, W. E. For- 
sythe, M. A. Easley —150(L) 
Tungsten mercury arc, spectral energy, B. T. 
Barnes—1468 
Spectra, molecular 
Absorption bands of Ke, W. O. Crane, A. Christy 
—421 
Atomic J values and molecular quantum numbers, 
R.S. Mulliken—1440 
Band in Be®H spectra, W. W. Watson—1019(L) 
Bands in the absorption spectrum of mercury, J. 
G. Winans—1020(L) 
CaF and SrF, interpretation of, A. Harvey, F. A. 
Jenkins—1413(L) 
Cameron bands in CO intensities, H. P. Knauss, 
J. C. Cotton——1099(L) 
CO third positive group, J. Kaplan— 7? 34(A) 
Electronic states in halogen bands, R. S. Mulliken 
—699 
Energy states, effect of interaction, E. C. Kemble, 
F. F. Rieke—153(L) 
Halogen bands, interpretation, R. S. Mulliken— 
364(L) 
H2, B-A resonance bands, quantum assignments, 
H.H. Hyman—187 
Intensities in bands, J. Kaplan—778(L) 
Isotope effect on intensities, J. L. Dunham—1553 
Liquids, absorption in infrared, effect of pressure, 
J. R. Collins—305 
Methane gas, new absorption band, theory of 
overtones, D. M. Dennison, S. B. Ingram—1451 
Nz, absorption bands, J. J. Hopfield—789(A) 
N2 and Os», abundances of Nyg, Orig and Oy7, S. M. 
Naudé—333 
Ne band spectra, by canal ray excitation, H. D. 
Smyth, E. G. F. Arnott—1023 
Notation, report on, R. S. Mulliken—611 
Oz bands, new, J. J. Hopfield—789(A) 
Optical excitation of CdH and ZnH bands, P. 
Bender—1543 
Organic liquids, absorption in infrared, high res- 
olution, R. B. Barnes—296 
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S:, vibrational and rotational analysis—S. M 
Naudé, A. Christ y—1801(L) 
Symmetrical diatomic molecules, intensities, 


theory, E. Hutchisson—410 
ZnH and CdH bands, Zeeman effect, W. W. Wat- 
son—1134 
Spectrograph, spectrometer 
X-ray, double crystal, J. W. M. 
Hoyt—1702 
Spectroscopic technique 
Unreliability of photographic emulsions for re- 
cording distances, improvements, D. Cooksey, 
C. D. Cooksey—80 
Sun (see Astrophysics) 
Surface tension 
Of saturated vapors, maximum, J. 
shefsky—377 (A) 
Retention of liquidsbet ween homogeneous spheres, 
W. O. Smith, P. D. Foote, P. F. Busang—524 
Susceptibility, magnetic (see Magnetic properties) 
Thermal expansion 
M-M-M alloy, P. Hidnert, W. T 
787(A) 
Thermionic emission of electrons 
Efficiency of emission from oxide coated filaments, 
B. J. Thompson—1415(L); erratum—1802(L) 
Oxide coated filaments, shot effect, H. N. Koza- 
nowski, N. H. Williams—1314 
Oxide coated Ni-Ba alloy, N. C. Beese-—-1309 
Shot effect, abnormal, of ions, J. S. Donal Jr. 
1172 
Work function and constant A for thoriated tung- 
sten, W. B. Nottingham—376(A) 


Dumond, A. 


L. Shere- 


. Sweeney 


Thermodynamics 
Cp Variation, calculated from p-v-t data, L. G. 
Hoxton—1091 
Entropy, elastic strain, W. S. Kimball—377(A) 
Entropy of crystalline hydrogen, L. Pauling—430 
Entropy of H.2, D. MacGillavry—1398; W. F. 


Giauque, H. L. Johnston—1592(L) 
“quations of state and thermodynamic functions 


for a substance with variable specific heat, W. 
P. Boynton—787(A) 

Equilibrium, ammonia synthesis reaction. I. Mass 
action function, L. J. Gillespie, J. A. Beattie 
743. Il. Heat of reaction entropy, L. J. Gilles- 
pie, J. A. Beattie—1008 

Fundamental correspondence to geometry, J. Q. 
Stewart—377(A) 

Gibbs-Dalton law of partial pressures, L. J. 
Gillespie—121 

Mechanical treatment of temperature distribution 
of radiation, W. Anderson—365(L); reply, R. C. 
Tolman—365(L) 

Real gases and mixtures, J. A. Beattie—132 

















Surface tension of saturated vapors, J. L. Shere- 
shefsky—377 (A) 
Temperature equilibrium in a static gravitational 
field, R. C. Tolman, P. Ehrenfest—1791 
Thermoelectric effect 
In crystals grown in magnetic fields, A. Goetz, M. 
F. Hasler—1752 
Tribothermoluminescence 
In crystals, glass and frits, R. E. Nyswander, B. 
E. Cohn—1257 
Ultrasonic waves 
Velocity and absorption in O2, W. H. Pielemeier— 
1005, 1667(L); J. C. Hubbard—1668(L) 
Velocity in water vapor, G. E. Thompson—77 
Vacuum tubes 
High voltage, M. A. Tuve, L. 
Dahl—1261(L) 
Low grid-current tube, type F. P.-54, G. F. Met- 
calf, B. J. Thompson—1489 
Viscosity 
Margules method, absolute values, H. R. Lille— 
347 
Vision theory of 
Electromagnetic theory of sight and color vision, 
S. R. Cook—790(A) 
Wave functions 
For any atoms, approximate, J. C. Slater—57 
For Be, B, C, N, O, F, Ne, C. Zener—S51 
Work function (see Thermionic emission) 
X-rays, absorption 
By Li, K. C. Mazumder—457 
Elements Ca to Ga, fine structure, B. Kievit, G. A. 
Lindsay—648 
In gases, W. W. Colvert—1614 
Penetration, therapeutic use, E. Q. Adams— 
1020(L) 
Soft rays in gases, C Ka in air, E. Dershem, M. 
Schein—378(A) 
X-rays, diffraction, reflection and scattering 
Compton effect, position and structure of modified 
line, F. L. Nutting—1267 
Compton shift, accurate measure, N. S. Gingrich 
—364(L) 
Electric fields and diffraction pattern of a liquid, 
G. W. Stewart—1413(L) 
Fine structure lines, with double crystal spectro- 
meter, N.S. Gingrich—1050 
Fine structure of unmodified line, J. A. Bearden 
—791 
Liquids, effect of temperature, E. W. Skinner— 
1625 
Multiple scattering in the Compton effect, J. W. 
M. Dumond—1685 


R. Hafstad, O. 
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Scattering coefficient as a function of wave-length 
and atomic number, E. N. Coade—778(L), 1109 
Scattering power of Ni, Cu, Fe, for various wave- 
lengths, R. W. G. Wyckoff—1116 
Thin metallic films, reflection from, J. E. Hender- 
son, E. B. Jordan—785(A) 
Width and structure of modified Ka, Kao, P. A. 
Ross, J. C. Clark—378(A) 
X-rays, emission 
Fine structure, Ka of Fe, Co, Ni, Cu, Mo and Ag, 
KB of Mo, J. Valasek—1523 
Hyperfine structure W Ka, F. K. Richtmyer, S. , 
W. Barnes, K. V. Manning—1017(L) 
Kas of Cu, origin of, and Richtmyer double group 
hypothesis, J. W. M. Dumond, A. Hoyt—799; 
J. W. Dumond—1015(L) 
Ta, L-spectrum relative intensities of lines, V. 
Hicks—785(A), 1273 
X-rays, polarization 
Continuous x-rays from single electron impacts, 
B. Dasannacharya—1675 
X-rays, reflection and refraction 
Glancing angle of calcite for Ag Ka, C. D. Cook- 
sey, D. Cooksey —85 
X-rays, scattering 
Breadth of modified line, J. W. M. Dumond— 
146(L) 
X-rays, spectra 
Fine structure in K-radiation of lighter elements, 
L. Y. Faust—161; M. Soderman—1414(L) 
X-rays, spectra and spectroscopy 
False lines in x-ray grating spectra, J. M. Cork 
—665 
Intensity ratios of satellites, F. K. Richtmyer, L. 
S. Taylor—1044 
Radiation from a high potential x-ray tube, C. C. 
Lauritsen—1680 
X-rays, tubes 
Currents through walls, D. P. LeGalley, W. R. 
Ham, M. W. White—379(A) 
Zeeman effect 
Auroral green line. R. Frerichs, J. S. Campbell— 
1460 
Incomplete Paschen-Back effect, calculation, J. B. 
Green—157, 379(A) 
In Co II spectrum, J. H. Findlay—5 
In quadrupole radiation, R. Frerichs, J. S. Camp- 
bell—151(L) 
Misplacements and g values, H. N. Russell— 
1590(L) 
ZnH and CdH bands, W. W. Watson—1134 
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